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I 


SUMMARY 


The scope of this program was to include a thorough critique of the 
JANNAF sub-critical propellant injection/combustion process analysis computer 
models and application of the models to correlation of well documented hot 
fire engine data bases. These programs are the Distributed Energy Release 
(DER) model for' conventional liquid propellant injectors and the Coaxial 
Injection Combustion Model (CICM) for gaseous annulus/liquid core coaxial 
injectors. The critique would identify model 'inconsistencies while the 
computer analyses would provide quantitative data on predictive accuracy. 

The program was comprised of three tasks; Task I - Computer Program Review 
and Operation, Task II - Analysis and Data Correlations, and Task III - 
Documentation. ' 

There were three objectives of Task I. (1) Critique of the DER and 
CICM Computer Programs, (2) Correction of coding errors, updating of inade- 
quate formulations, and addition of diagnostic printout statements, and 
(3) Identification of inconsistencies between the analysis computer programs 
and the JANNAF prediction procedures documented in CPIA 246. The results of 
the DER and CICM reviewsare comprehensively reported in Appendices A and B, 
respectively. Complete summaries of the corresponding conclusions and 
recommendations of -the reviews are contained in Section III, Computer Pro- 
gram Review and Operation. There were two major conclusions resulting from 
the DER review. First, the intended predictive accuracy of the JANNAF rigorous 
performance evaluation procedure (to within 1 percent for predicted specific 
impulse) is, in general, currently out of the question for a priori performance 
prediction with DER. Secondly, the DER analysis originally planned to be 
conducted during program Task II should rather be concerned with improvement 
of a DER technical shortcoming. The primary conclusion of the CICM review 
was that the applicability and accuracy of the model is currently limited by 
the absence of an intra-element coaxial gas/liquid mixing model. This limita- 
tion not only makes the mixing loss calculation dependent on correct appli- 
cation of empirical cold flow mass distribution data, but hinders the develop- 
ment of general program coaxial jet atomization and drop size constants that 
control the program vaporization calculation. 


- 1 - 



I Summary (cont.) 

There were originally three primary objectives of Task II. (1) Provide 
information on the present prediction capabilities of the JANNAF DER and CICM 
injection-combustion computer analysis techniques, (2) Identify conditions 
where rel iable -predictions can be obtained, and (3) Identify areas requiring 
further improvement and research. The CICM analysis task was completed as 
.originally planned. The results of the CICM analysis are reported in Section 
IV, CICM Analysis and Data Correlations. The CICM analysis was performed by 
establishing the existing M-l H^/Og engine data base, executing a nominal 
operating point CICM analysis, correlating the CICM prediction with the test 
data, conducting two off-nominal test point analyses to determine the influence 
of velocity ratio changes on injector performance, and identifying prediction 
ranges and required model improvements. The CICM analysis results verified 
the accuracy of the CICM vaporization model for the case where injector intra- 
element mixing losses are negligible. 

The objective of the DER Phase of Task II was altered based on the 
recommendations of the Task I DER computer model review. Improvement of the " 
LISP subprogram ZOM plane mass distribution and mixing methodology was 
selected as the new Task II DER goal. This task was conducted in four parts. 

(1) An a priori ZOM plane prediction model was formulated that accounts for 
combustion gas acceleration effects on inter-spray fan mixing, (2) A subscale 
test data base was developed for analysis and the ZOM model was used to predict 
mixing performance for each test, (3) The model predictions were correlated 
with the hot fire test results, and (4) Recommendations for continuation of 
model development were formulated. The primary discovery of this initial 
ZOM model development work was that a physically mechanistic near-zone model 
that will predict the ZOM mixing plane location must account for both gas 
acceleration and reactive stream ("blowapart") forces on droplet spray fan 
formation and mixing. 

Task III of the program resulted in eleven monthly status letters 
and this comprehensive final report containing explicit recommendations for 
improvement of the JANNAF performance prediction computer programs. The 
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I Summary- (cont. ) 

English system of units has been exclusively employed inthis report since 
SI units have yet to be adapted to the JANNAF system .of computer programs. 
The program COR has concurred with and approved this choice. 
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II 


INTRODUCTION 


The rCRPG (now JANNAF) Performance Standardization Working Group 
was formed in 1965 for the purpose of improving and recommending methodology 
for the analytical and experimental evaluation of the performance of liquid 
propellant rocket engines. In 1968, the working group published a Performance 
Evaluation Manual (Ref. 1) which described the procedures and computer 
programs recommended for the prediction, correlation, and extrapolation of 
the performance of liquid propellant thrust chambers. The scope of this 
first effort was limited to assembling, into a compatible overall system, 
the best relevant analytical and experimental techniques existing through- 
out the industry at that time. During this effort, it was concluded that 
the energy, release phenomenon could not be adequately described or predicted 
by existing analytical techniques. As a result, an interim empirical proce- 
dure was. adopted. 

Since this first attempt at achieving a standard performance evalua- 
tion model, a semi-empi'rical , but mechanistic, computer model has been 
developed for the analysis of the liquid injector-combustion chamber energy 
release process. This model, termed the Distributed Energy Release (DER) 
model (Ref. 2) has reached the stage of development where it is being in- 
corporated into the Improved JANNAF Performance Evaluation Methodology (Ref. 

3). DER is composed of two major programs which link the atomization, vapor- 
ization and mixing processes within the combustion chamber. The first is the 
Liquid Injector Spray Patterns (LISP) program which calculates propellant 
mass and mixture ratio distributions at a specified chamber cross-sectional 
plane (ZOM) downstream of the injector face. The second is the Stream Tube • 
Combustion (STC) program which calculates the propellant vaporization, 
reaction and acceleration from the LISP specified collection plane to the 
combustion chamber throat plane. Additionally, a third JANNAF recommended 
program has been developed for the specialized case of injector elements 
containing central circular orifice liquid propel lant 'injection surrounded 
by annular gaseous injection. The Coaxial Injection Combustion Model (CICM) 
(Ref. 4) is designed to replace the DER LISP subprogram for this injector 
type. 
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II 


Introduction (cont.) 


While these programs provide analytical methods for evaluation of the 
energy release process, the program developers have identified analysis 
parameters which are critical to the accuracy of the resulting performance 
predictions. These include specification of propellant mass median droplet 
diameters and the LISP Spray distribution correlation coefficients, which 
have been established over limited ranges of element type and design condi- 
tions. Additional studies using DER have shown that the specification of 
the LISP-STC interface plane (ZOM) is also critical to the end performance 
prediction. 

The objective of this program was to develop quantitative data on 
the present prediction capabilities of- the JANNAF sub-critical propellant 
injection/combustion process analysis programs (LISP, STC, and CICM). The 
desired program end product was identification of conditions for which 
reliable predictions could be conducted and. areas which need further improve- 
ment and research. 

Future attainment of a broader overall objective was continued with 
conductance of the Injection Processes Program. The JANNAF Performance 
Standardization Working Group has the purpose of improving methodology 
for analytical design modeling of rocket engines. The current and future 
economics of rocket development do, and will certainly, make it imperative 
that cost saving analytical methods replace more expensive hardware develop- 
ment and test programs. Of course, such tools are only cost effective if 
they- model the applicable physical processes realistically and accurately. 

The Injection Processes program and other related efforts have provided - 
information on the state of JANNAF model development through application 
to real rocket engine systems. During this program the CICM computer program 
was used to correlate performance data obtained with the M-l 1 million Ibf 
hydrogen/oxygen engine. The DER computer program has been successfully 
applied to design analysis of the Orbital Maneuvering System (OMS) engine 
for Space Shuttle, the Improved Transtage Injector Program (ITIP) currently 
being conducted by the USAF, and an advanced development monomethyl hydrazine/ 
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II 


Introduction (cont.) 


fluorine-oxygen engine tested by the NASA. Each of these efforts has resulted 
in constructive criticism of the computer models that, when applied, results 
in further advancement of the state-of-the-art of rocket engine analytical 
design. The final end product of programs that support the JANNAF predictive 
methodology will someday be a capability to eliminate major hardware develop- 
ment technology programs through verified standardized analysis techniques. 

A superior development procedure would be constituted of initial JANNAF model 
analysis, fabrication and test of the full scale engine, re-analysis, full 
scale hardware modification, and final engine verification test. The In- 
jection Processes. Program has made this seemingly optimistic goal a bit more 
achievable through a comprehensive evaluation of the DER and CICM models. 
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COMPUTER PROGRAM REVIEW AND OPERATION 


There were three primary objectives of the first program task. 

(1) Critique of the JANNAF DER and CICM programs, 

(2) Correction of codir-g errors, updating of inadequate formulations, 
and addition of diagnostic printout statements, and 

(3) Identification of inconsistencies between the analysis 

computer programs -and the JANNAF prediction procedures described in CPIA 
246 (Ref. 3). ' 

The complete results of the DER and CICM reviews are contained in Appendices 
A and B, respectively, of this report. The computer programs are introduced 
and their functions in the JANNAF performance prediction procedure briefly 
described in the following paragraph. A complete summary of the findings 
and corresponding recommendations of the computer model reviews follows 
the program descriptions. 

A flow chart showing the DER and CICM programs and their relationship 
to the JANNAF Two-Dimensional Kinetic (TDK) Computer Program (Ref. 5) is 
illustrated in Figure 1, taken from Ref. 3. DER is composed of LISP and STC, 
two major programs- that link atomization, vaporization, and mixing processes 
within the combustion chamber. The Liquid Injector Spray Patterns (LISP) 
program calculates propellant mass and mixture ratio distribution at a speci- 
fied chamber cross-sectional plane (termed ZOM) downstream of the injector 
face. LISP was developed for conventional (i.e., circular orifice) liquid/ 
liquid injection 'elements. The Stream Tube Combustion (STC) program calculates 
propellant vaporization, reaction, and acceleration from ZOM to the combustion 
chamber throat plane. STC can provide direct computer input data for the TDK 
program that continues the multiple stream tube analysis through the supersonic 
expansion process. CICM replaces the LISP program for the analysis of gas/ 
liquid coaxial elements. CICM is a highly specialized program that has currently 
only been applied- to the analysis of injection elements with a central liquid 
^2 circular core surrounded by a gaseous t-^ or l-^/C^ combustion gas mixture 
annulus. 
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Ill Computer Program Review and Operation (cont.) 

A. DER. Computer Model Review Recommendations and Conclusions 

Four subtasks were accomplished during the DER review. 

(1) Identification and Correction of Coding Errors, 

(2) Addition of Diagnostic Comment Cards and Print-Out 
Statements, 

(3) Identification of Inadequate Formulations and Model 
Technical Formulations, and 

(4) ' Review of the JANNAF Performance Prediction Procedures 

(CPIA 246) with Regard to Use of DER. 

The review is applicable strictly the DER subcritical K-Prime version des- 
cribed in Ref. 2. The corresponding user's manual referred to in this 
report is Ref. 6. 

The third subtask listed above was emphasized during the review 
for two reasons. The initial results of the review indicated that DER still 
requires major technical improvements and therefore subtasks (1) and (2) 
were considered to be of less current interest. Secondly, SDER, a new 
"standardized" version of DER (Contract FO 461 1 -75-C-0055),was developed 
concurrently with completion of this program. It was intended that the 
improved DER model be influenced by the findings surranarized in this report; 
therefore the discovery of DER technical formulation shortcomings was con- 
sidered to be of prime importance. 

A major conclusion of the DER review was that the DER analysis 
originally planned to be conducted during program Task II should rather be 
concerned with improving a DER technical shortcoming. It seemed inappropriate 
to conduct the analysis with a computer model that possessed vaporization 
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Ill Computer Program Review and Operation (cont.) 

and mixing models containing several questionable solution formulations, 
as summarized in the following paragraphs concerning review recommendations. 
Improvement of the LISP ZOM plane mass distribution methodology was selected 
as the new Task II DER analysis goal. The current status of the mixing model 
improvement work is described in Section V of this report. Key recommenda- 
tions and conclusions, resulting from the DER review results detailed in 
Appendix A, are listed in the following four paragraphs corresponding to 
the previously described review subtasks. 

1 . Identification and Correction of Coding Errors 
a. LISP Subprogram 

(T) An unsymmetrical pie section input problem 
was identified for the LISP program. It should be eliminated by adjusting 
the collected p.ie section mass flowrate to 0/360 of the total injected flow 
of each propel.! ant. 


(2) Inconsistencies between published DER drop 
size equations and those actually existent in the DER code must be resolved. 

(3) The DER code should be changed to eliminate 
a mass flux calculational error for triplet elements caused by an improper 
rotation of the ZOM collection plane around the normal x axis. 

(4) The ZOM mass distributions should consider 
the influence of .baffle height. 

b. STC Subprogram 

(1) The STC program limits the number of radial 
and circumferential mesh lines to twenty; this limitation should be noted in 
the DER user's manual, or preferably removed. 
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Computer Program- Review and Operation (cont.) 


2. Addition of Diagnostic Comment Cards and Printout 
Statements 


The recommended statement additions and improvements 
are presented in Section B of Appendix A. 

3. - Identification of Inadequate Formulations and Model 

Technical Shortcomings 

a- Drop Size Prediction 

(1) The inconsistencies cited, between referenced 
drop size correlations and those appearing in the DER code, must be resolved. 

(2) It is recommended that the DER drop size 
equations be comprehensively reviewed with respect to available atomization 
correlations and their impact on DER performance prediction accuracy. A task 
performed during the SDER development program was to be concerned with such 

a review, although the results have not been published. 

(3) Interim to release of SDER, all DER drop 
sizes should be user input and justified. 

b. ZOM Plane Selection 

(1) The ZOM point source flow assumption should 
be tested empirically. That is, it should be determined if the LISP spray dis- 
tribution coefficients are a function of the cold flow collection plane dis- 
tance. 


(2) The ZOM mass distribution methodology should 
account for combustion effects such as gas acceleration and reactive stream 
separation forces. A proposed model approach is detailed in Section V of this 
report. 
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Computer Program Review and Operation (cont.) 


. (3) The LISP spray coefficient matrix should 

be expanded if the ZOM technique is retained- in DER. 

c. DER Vaporization Sensitivity Study 

(1) The implications of the work of Bracco (Ref. 
7) with respect to DER vaporization modeling should be evaluated. 

(’2) The DER K-Prime vaporization model insen- 
sitivity to chamber pressure should be investigated. The argument suggested 
in Appendix A to be the source of this error should be evaluated. 

(3) The DER integration technique droplet 
downstream station velocity error should be eliminated. Additionally, the 
Euler predictor-corrector technique should be evaluated through a study 
using different calculational step sizes and number of corrective iterations. 
The possibility of developing a more efficient integration technique should 
be investigated. 


(4) The results of this study and the work 

of Bracco- both indicate the importance of the droplet drag coefficient (Cp) 
assumption. The drag coefficient literature shbuld be reviewed, and the 
selected DER drag coefficient formulation justified. 

(5) The DER vaporization model should account 

for droplet heatup. 

(6) The DER user manual and CPIA 246 should 
include an expanded section on droplet size distribution input selection. 

d. Near-Zone Combustion and Monopropellant Flame 
Considerations 

(1) It is recommended that DER incorporate 
a monopropellant flame model for reasons cited in Section C.4. of Appendix A. 
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Ill Computer Program Review and Operation (cont. ) 

e. Combustion Gas Acceleration and Reactive Stream 

Separation (RSS) Effects on Cold Flow Mass 

Distribution 

(1) It is recommended that a RSS model be 

considered for DER. 

(2) The initial development of an a priori 

. ZOM plane selection methodology (See Section V) should be brought to fruition. 

f. Turbulent Mixing Model 

(1) The characterization of turbulent mixing 
effects in DER would comprise a large step toward providing DER with the 
desired a priori prediction capability. It is recommended that such a model 
be considered for DER. 

g. Development of an A Prior i DER Mixing Model 

(1) It is recommended that the current LISP ZOM 
model be improved by incorporating the influences of combustion gas accelera- 
tion, reactive stream separation, and turbulent mixing. As previously men- 
tioned, an a priori ZOM calculational technique is also required. This topic 
is expanded in Section C.7. of Appendix A. 

4 • Inconsistencies Between JANNAF Procedures and DER 

Computer Program- Operations 


The primary conclusion is that the intended predictive 
accuracy of the JANNAF (DER) rigorous procedure (to within 1 percent for 
predicted specific impulse) is currently out of the question for a priori 
performance prediction. This directly relates to the program decision to 
forego the originally planned Task II DER analysis and concentrate, instead, 
. on improvement of the ZOM plane mass distribution methodology. 
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Ill Computer Program Review and Operation (cont.) 


B. CICM Computer Model Review Recommendations and Conclusions 
The CICM review was accomplished in three subtasks. 

(1) Identification of Operational Problems Including a Code 
Review and Inclusion of Diagnostic Print-Out Statements, 

(2) Identification of Inadequate Formulations and Model 
Technical Shortcomings, and 

(3) Review of the JANNAF Performance Prediction Procedure 
(CPIA 246) with Regard to the Use of CICM and Identification of Inconsistenci 


The review is applicable to the CICM version described in Ref. 4, which also 
contains the user's manual referenced continually in this report. 


The review was initiated by executing the program documented 
sample case and attempting to interface the program output with the STC 
subprogram of DER, as recommended in CPIA 246 for gas/liquid coaxial injector 
rigorous performance analysis. It was determined that the current CICM 
interface routine, DERINI, was incomplete and punched several improperly 
formated cards for input to the STC subcritical K-Prime version. First 
priority, during the review, was given to development of a new CICM/STC 
interface procedure because of the need for an accurate and cost-effective 
method of interfacing CICM and STC during the program Task II CICM analysis. 
The resulting new procedure is detailed in Section C.3. of Appendix B. The 
key recommendations and conclusions resulting from the CICM review results 
detailed in Appendix B are listed in the following three paragraphs corres- 
ponding to the previously described review subtasks. 

1 ■ Coding Errors and Diagnostic Statements 


It is recommended that the CICM calculational problem 
that results in periodic "dropping" of drop size groups from the calculation 
be investigated. 
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Ill Computer Program Review and Operation (cont.) 

2* Identification of Inadequate Formulations 
and Model Technical Shortcomings 

The identification of inadequate CICM formulations 
and technical shortcomings was considered to be the next most important review 
task after improvement of the CICM interface procedure. CICM is a relatively 
new JANNAF program that has not been used extensively, except by the developers 
of the model. Therefore, it was considered important that basic model assump- 
tions and analysis techniques be critically evaluated. The recommendations 
and conclusions resulting from the CICM technical formulations review are 
summarized below. 


a. A review of the CICM stripping rate correlation 
should be conducted. The derivation of the current, or any proposed alternate 
correlation, should be substantiated and be made open to critical review. 

b. A review of the CICM drop size correlation should 
be conducted. Such a study could also investigate the sensitivity of coaxial 
injector performance to the predicted jet mass median drop size. This would 
allow deterni nation of the performance prediction uncertainty due to the 
availability of many different drop size correlation equations. 

c. The drop size distribution tabulated at the end 
of a CICM run i,s only the summation of several constant mass median diameter 
groups; each group being calculated over a particular axial step. This 
resultant distribution is quite different than a drop size group calculated 
with distributions typically used to model rocket combustor sprays (e.g., 
Nukiyama-Tanasawa, Logarithmic-Normal, etc.). It is recommended that the 
significance of this CICM model simplification be evaluated. 

d. It 'is strongly recommended that the CICM technique 
for accounting for intra-element mixing be improved. If the use of single 
element cold flow data to specify the intra-element mass distribution is ' 
continued, a standard measurement technique should be developed. A standard 
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Ill Computer Program Review and Operation (cont.) 


methodology for interpreting and inputting the data to CICM is also required. 
Preferably,- an intra-element mixing model should be developed for CICM. 
Applicable models have been derived from experiment for gas/gas coaxial 
element mixing. The first step in adapting such models would be to determine 
the feasibility of applying a gas/gas mixing model to the solution of gas/ 

1 iquid mixing. 

e. All JANNAF engine analyses should record estimated 
manifold maldistribution performance losses, to build up a reference data 
base. 


3. Inconsistencies Between JANNAF Procedures 

and Program Operations 

The new CICM/STC interface procedure was written during 
this review subtask. The recommendations and conclusions resulting from the 
review of CICM's role in the JANNAF performance procedures are listed below. 

* 

a. The original provision of the CICM/STC interface 
was for the supercritical DER program version. The new CICM/STC interface 
procedure described in Section C.3. of Appendix B should be used for sub- 
critical propellant analysis. This procedure should also be adopted for use 
in the new "standardized" DER program currently being developed. 

b. The CICM and STC programs should be interfaced 
at a chamber axial plane where all the calculated oxidizer drop size groups 
have been heated to the chamber "wet bulb" temperature. 

c. A standard JANNAF procedure or technique should 
be developed to predict single coaxial element intra-element mass distribu- 
tion. 


d. A procedure should be developed for allowing -for 
the effect of diffusion mixing on face plane measured manifold mass distri- 
butions. 
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Ill .Computer Program Review and Operation (cont.), 

e. An accurate CICM mass distribution analytical 
model or empirical approach is required to allow OANNAF standard atomi- 
zation coefficients (C^ and B^) to be backed out from coaxial injector hot ' 
fire data. 


- 17 - • 



IV 


CICM ANALYSIS AND DATA CORRELATIONS 


The original objectives of Task II were: (1) Provide information on 

the present prediction capabilities of the JANNAF DER and CICM' injection- 

combustion computer programs; (2) Identify conditions where reliable predictions 

can be obtained; and (3) Identify areas requiring further improvement and 

research. The CICM phase was completed as originally planned, while the DER 

phase of the task-was rescoped (see Section V). The CICM model was applied to 

correlation of characteristic exhaust velocity efficiency (n^*) for three 

■tests conducted with the M-l pressure fed 600,000 Ibf (at 550 psia chamber 

pressure) hydrogen/oxygen engine. The CICM analysis was limited to tests with 

subcritical liquid oxygen inlet conditions. Excellent agreement was obtained 

between op* and n r * from the JANNAF simplified prediction methodology 
U TEST. . L PRED 

for two of the three tests analyzed. The results of the analysis have verified 
the accuracy of the CICM model for the case where injector intra-element mixing 
losses are negligible. 

A. M-l Engine Experimental Data Base 

The data base selected for the analysis and correlation of the 
CICM computer program was that of the M-l thrust chamber developed by ALRC under 
NASA Contracts NAS 3-2555 (Ref. 8) and NAS 3-11214 (Ref. 9). The M-l engine 
was designed to utilize liquid oxygen/liquid hydrogen propellants and deliver 
1,500,000 of thrust when operating at its nominal design conditions of 1000 
psia chamber pressure and 5.49 mixture ratio. During development, the thrust 
chamber was tested with LO^/GHg propellants with a low area ratio ablative 
combustion chamber over a range of chamber pressure (550-1050 psia), mixture 
ratio (4-6), and hydrogen inlet temperature (80-130°R). The CICM data base 
met all the pre-defined program requirements for the following eight reasons: 

1. Conventional injector element applicable to CIC-M (gas/ 
liquid coaxial); 

2. Capable of direct modeling with CICM/DER; 

3. Subcritical propellant conditions (P c = 550 psia); 

4. Propellants of future interest (0 2 /H 2 ); 
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CICM Analysis and Data Correlations (cont.) 


5. Low area ratio test configuration (e = 2:1); 

, 6. - Simple wal.l boundary conditions (no mass addition, 

minimal fufe-1 film, cooling of 1/2 percent of the 
total flow rate)-; 

7. ' - Test.data at -nominal and off-nominal operating conditions 

(0/F, hydrogen density variations); 

8. ' Element to element mass distribution cold flow, data. 

Detailed descriptions of all the-M-1 test hardware, facilities, 
and' data measurement techniques are contained within the JANNAF-Simpl ified 
Performance Prediction narrative of Appendix C. The S/N 012 injector analyzed 
during the study is pictured in Figure 2. The injector contained 3,248 elements 
with gaseous hydrogen being injected annularly around the oxidizer. A row of 
360 orifices drilled through the porous > igimesh face were located around the 
injector periphery and provided the chamber wall film cooling. Approximately 
3.7 percent of the total fuel flow rate was used for chamber wall film cooling. 
Total .fuel element flow rate was 89.. 8 percent of the thrust chamber fuel flow- 

rate with a baffle fuel film cooling flow percentage of 3.9 percent. The re- 

maining 2.6'percent of the fuel flowed through the rigimesh injector face. The 
coaxial element consisted of two basic components which were threaded together. 

An oxidizer tube was recessed within the fuel sleeve producing a fuel annulus 
between the two -parts. The oxidizer tube was flared at a fifteen degree half 
angle- and was. recessed 0.231 inches from the injector face. Elements were 
arrayed in 33 concentric rows. The low area ratio combustion chamber used for 
testing with the M-l injector was comprised of an outer steel shell and an 
inner ablative liner (tape wrapped silica-reinforced phenolic). The assembled 
combustion chamber (See Figure C-4 of Appendix C) consists of an upper fuel 
torous and a lower conical combustion chamber. 

The test data that was reduced during the task data evaluation 

effort is' tabulated in Table I. Nomenclature for Table I is shown in Figure C-l 

of Appendix C. The three tests that were selected for CICM analysis are detailed 
in Table II. Test 009 was at the nominal operating point. Test 010 was analyzed 
to investigate the influence of mixture ratio on performance. Test 016 was 
analyzed to correlate the effect of injection velocity ratio change due to 
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TABLE I 


SEA LEVEL SUB-CRITICAL TEST DATA 


t 


Test No. 

Summary T<me 

Duration 

Throat Area 

Chamber 

0/F 

Thrust 

Isp 

Meas. 

( s .ec) 

305.0 

U 

Wp 

“fell 

w_. 

W T 

(#/sec) 


Start 

(sec) 

Stop 

(sec) 

(sec) 

Pre 

(InZ) 

Post 

(in 2 ) 

Pressure 

(psia) 

Meas. 

. IIM.. 

492840 

(#/sec) 

(#/sec) 

FT 

(#/sec) 

007 

44.2 

44.7 

44.72 

707.370 

711.860 

582.8 

4.87 

1340.7 

249.4 

26.0 

275.3 

<616.0 

009 

44.3 

44.8 

44.81 

728.269 

735.994 

556.6 

5.46 

495409 

300.5 

1393.4 

220.2 

35.0 

255.2 

1648.6 

010 

46.8 

47.3 

47.33 

735.994 

736.308 

572.0 

4.04 

510096 

310.4 

1317.6 

296.4 

29.6 

325.9 

1643.6 

014 

46.8 

47.3 

47.38 

706.495 

722.048 

541.1 

5.30 

481765 

303.5 

1335.4 

213.4 

38.6 

251.9 

1587.3 

016 

45.0 

45.5 

45.56 

722.048 

727.902 

567.9 

5.53 

501 304 

301.7 

1407.1 

209.0 ( 

45.6 

254.6 

1661.7 

017 

46.3 

46.8 

46.89 

727.902 

728.368 

571.0 

4.76 

506116 

307.5 

1360.1 

245.7 

40.1 

285.7 

1645.8 

019 

44.3 

44.8 

44.89 

733.644 

736.391 

576.0 

5.15 

516590 

304.4 

1421.3 

236.1 

39.9 

276.1 

1697.3 

020 

46.5 

46.5 

46.5 

736.391 

748.222 

569.4 

5.07 

510642 

298.7 

1428.1 

240.0 

41.6 

281.5 

1710.0 


Test 

PFT 

IP s .|a) 

PFFM-2 

(psia) 

PF Mt v-2 

iM) 

PFTCV-1 

jpsTa) 

PFTCV-2 

is ilil 

PFJ-3A 

(psia) 

TFFM 

.m. 

TFTCV-2 

TFJ 

(°R) 

POT POFM P0TCV-1 
(ps1a)(ps1a)(ps1a) 

P0TCV-2 

(psia) 

P0J-2A 

(psia) 

TOFM 

(°R) 

171 

TOTCA-2 

(°R) 

TOJ PC4B-1 PC48-2 
(°R) (psia) (psia) 

007 

805 

731 

722 

719 

703 

624 

44 

102 

84 

749 

729 

729 

724 

680 

186 

173 

482.4 

482.6 

009 

808 

748 

740 

741 

720 

619 

44 

117 

97 

750 

732 

729 

737 

674 

168 

181 

169 

464.4 

463.8 

010 

878 

773 

761 

763 

742 

638 

45 

89 

82 

749 

737 

734 

737 

685 

173 

177 

174 

477.5 

476.9 

014 

832 

778 

753 

763 

746 

523 

45 

116 

110 

730 

720 

717 

705 

662 

173 

180 

174 

451.5 

450.1 

016 

872 

823 

805 

808 

788 

646 

44 

127 

122 

769 

750 

746 

734 

686 

173 

181 

174 

474.0 

472.3 

017 

897 

831 

804 

812 

787 

652 

45 

108 

106 

759 

686 

742 

732 

685 

170 

181 

171 

476.7 

475.4 

019 

899 

83C 

811 

816 

792 

658 

45 

117 

no 

788 

769 

762 

740 

700 

171 

179 

172 

480.4 

478.3 

020 

900 

832 

814 

316 

793 

656 

44 

115 

107 

787 

769 

762 

753 

706 

169 

180 

170 

475.3 

473.7 
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TABLE II 

M-1 TESTS SELECTED FOR CICM ANALYSIS 


TEST 

W 

0 

(Ibrn/sec) 

W F 

(Ibm/sec) 

' T o 

(°R) 

A 

(°R) 

0/F 

P c 

(psia) 

V V o 

AV 

(ft/sec) 

P F 

(lbm/ft 3 ) 

n c* 

009 

1393 

• 255.2 

169 

97 

5.46 

524 ' 4 

18.2 

310 

1 .45 

.959 

010 

1318 

325.9 

174 

82 

4.04 

538 

16.2 

264 

2.16 

.964 

016 

1407 

254.6 

174 

122 

5.53 

534 

25.8 

456 

1.0 

.980 


009 Nominal Conditions 

010 Effect of Fuel Gas Density at Constant A V 

016 Effect of AV 



IV CICM Analysis and Data Correlations (cont.) 

hydrogen density variation. 

B. M-l Coaxial Injector Analysis with JANNAF Simplified 

Prediction Procedure 

The procedures and results of the CICM analysis of the M-l 
engine tests are summarized in the following three subsections, that describe 
in turn: (1) calculation of test characteristic exhaust velocity efficiency; 

(2) prediction of C* efficiency with the JANNAF simplified performance evalua- 
tion methodology; and (3) determination of test measured C* uncertainties. 

The JANNAF simplified prediction procedures described in CPIA 246 were utilized 
to economize and speed the analysis. 

Examination of the DER and CICM review results previously 
presented in Section III can, admittedly, lead to the conclusion that the M-l 
performance analysis described below has been conducted with inadequate models. 

An important consideration was the fact that the M-l thrust chamber design is 
very similar to the J2-S design used to calibrate key CICM jet stripping rate 
and drop size constants. (See Ref. 6 and J2-S sample case in CPIA 246). Also, 
both the M-l and J2-S engines posses extremely long chambers that eliminate 
significant intra-element mixing losses. Therefore, the M-l predictions were 
not invalidated by assuming uniform intra-element mass distribution, as described 
in a following paragraph. Additionally, using the STC subprogram of DER down- 
stream of CICM was -not considered an analysis weakness because STC utilizes 
similar key vaporization model analytical techniques to those of CICM (e.g., 
both models use the same droplet drag coefficient model). It should be remem- 
bered that a primary objective of the analysis was to verify that an independent 
user of the CICM/STC JANNAF analysis methodology could obtain an accurate 
performance prediction for a gas/liquid coaxial injector. 

1 . Calculation of Test C* Efficiency 

Test C* was calculated from the equation shown below, 
taken from Section 2.1.2 of CPIA 245. 


L TEST 


Pc eff A T 


TEST 


M 

t test 
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IV Cl CM Analysis and Data Correlations (cont.) 

Pc eff is the effect ive throat stagnation pressure, calculated from available 
chamber static pressure measurements. Two static pressure measurements were 
taken; at the Pc^ and Pc^ locations shown in Figure C-2 of Appendix C. The 
chamber combustion total pressure loss resulted from the CICM/STC computer 
run executed during the C* prediction analysis described in the next section. 
The CICM/STC calculated chamber static pressure profile correlated extremely 
well with the measured static pressures, as explained in Section IV.C.l. 

This correlation verified the CICM/STC calculated combustion (Rayleigh Line) 
total pressure loss. The test summary periods for analysis were selected to 
occur just prior to test FS2 so that the post-test ablative chamber throat 
diameter measurement would result in an accurate test throat area value. 


Test C* efficiency is simply the ratio of the test C* 
to the theoretical ODE C* value at the test propellant inlet, mixture ratio, 
and chamber pressure conditions. 


'C* 


TEST 


C* 

L TEST 

r* 

L ODE 


( 2 ) 


C* ODE was calculated with JANNAF TDK computer program (Ref. 5) _ at the test 
.conditions indicated in Table III. The resulting test C* efficiencies are 
also shown in Table III. 


2 . JANNAF Test C* Prediction 

The JANNAF simplified performance prediction methodology 
described in Section 3 of-CPIA 246 was utilized. Appendix C of this report 
contains a narrative of the application of the procedure to analysis of the 
selected M-l testsand sample input for all the JANNAF computer programs exe- 
cuted. The predictive equation for C* is expressed in terms of efficiencies 
for the significant chamber loss processes. 


'C* 


Pred 


- n 


C* 


HL 


x n 


C* 


TD 


x n 


C* 


KIN 


x n 


C*, 


BL 


x n 


C* 


MIX 


x n 


C* 


VAP 


(3) 
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TABLE III 


TEST 

0/F 

PC eff 

(psia) 

T o 

(°R) 

009 

5.46 

' 514 

169 

010 

4.04 

532 

174 

016 

5.53 

534 

174 


t 

IV) 

cn 


TEST 

CONDITIONS FOR n G * 

CALCULATION 

TEST 

H 

If 

H f 


(°R) 

0 

(cal/q-mole) 

(cal/q-mole) 

97 

-3027 

-1827 

82 

-2991 

-1918 

122 

-2991 

-1733 


r* 

u ODE 

b TEST 

n c * 

(ft/sec) 

• (ft/sec) 

TEST 

7694 

7376 

.959 

7960 

7674 

.964 

7685 

7529 

.980 
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CICM Analysis and Data Correlations (cont.) 


The purpose of the M-l test data analysis was to verify the capability of the 

CICM model to calculate the n r * (mixing) and n r * (vaporization) effi- 

_ MIX- L VAP 

ciencies for a GI^/LO^ coaxial injector. The meaning of and the technique used 
to evaluate each of the efficiency terms are explained in the following six 
paragraphs. 


a. Heat Loss Efficiency (n r * ) 

L HL 

The chamber heat loss efficiency was assumed to be 
I • 0 for each test. This assumption was made for two reasons.- (1) The thrust 
chamber wall was composed of an ablative silica-reinforced (tape-wrapped) 
phenolic that resulted in an effective adiabatic wall condition; and (2) Chamber 
heat loss to the injector face would be directly transferred to the propellants 
because of the plenum manifolds on the injector face backside. 


b. Two-Dimensional Flow Efficiency (o r * ) 

L TD 

The two-dimensional C* flow efficiency accounts for 
the reduction of -the throat potential flow area due to inlet effects. The 
equation used is simply the inverse of the inviscid flow discharge coefficient. 


M, 


'C* 


TD 


ODE 


TDE 


J 

C ' 

U INV 


( 4 ) 


The JANNAF ODE and TDE programs contained in TDK calculated the M-l chamber 
tiq* value of 1.002 (C^ = 0.998). This high throat C^ value occurs because 
of ffie large M-l chamber throat inlet radius ratio value of 2.132. 


c. Reaction Kinetic Efficiency ( n r * ) 

L KIN 

The reaction kinetic C* efficiency was calculated 

with the ODK option of the TDK program. For all mixture ratios from 1.0 to 

12.0 n r * was calculated to be 1.0 for the M-l engine. This occurs because 
L KIN 
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Cl CM Analysis and Data Correlations (cont.) 


of the high operating chamber pressure and thrust level of the engine (550 psia 
and 500,000 Ibf, respectively). 

d. Boundary Layer Efficiency (n r * ) 

L BL 

The C* boundary layer efficiency accounts for 
the displacement boundary layer effect on the throat potential flow area. 


>C* 


BL 


A 


T 


Ay - 2ll -Ry S*y 


(5) 


The TDK program was run at the Test 009 nominal 0/F to establish edge conditions 
for a boundary analysis with the JANNAF BLIMP computer program (Ref .10) . Wall 
temperature and calculated ablative chamber regression rates documented in 
Ref. 9 were used to establish input for BLIMP. BLIMP was executed by using 
the assigned wall temperature and assigned blowing rate input options, and 
edge gas properties for a mixture ratio of 2.5:1. This mixture ratio is the 
nominal Test 009 wall mixture ratio, based on M-l injector manifold mass dis- 
tribution results described in the next paragraph. The BLIMP calculated throat 
displacement thickness was -5. x 10 6 ft which resulted -in n c * of 1.000.’ Since 
the boundary layer effect on 0* was found to be small, this v§^Tue was assumed 
to be correct for all three tests analyzed. 


e. Mixing Efficiency (n r * ) 

L MIX 

The purpose of the M-l data analysis is to verify 
the capability of the JANNAF CICM computer program to predict energy release 
efficiencies for GH 2 /LO 2 coaxial injectors. The C* energy release efficiency 
is composed of a mixing and vaporization term. 


is shown below. 


'C* 


ERL 


•'V* x 
L MIX L VAP 


The C* simplified mixing efficiency definition 
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( 7 ) 


CICM does not calculate intra-element (shear) or inter-element (diffusion) 

mixing, however, the program has the capability to accept multiple zones of 

varying mixture ratio and to calculate the corresponding effect on the LO^ 

atomization and vaporization rates. Since CICM simply solves the equation 

shown above for n r * , this calculation was evaluated externally from the 
• MIX 

CICM program to allow inexpensive parametric evaluation of the M-l injector 
mass distribution data. 


The M-l injector manifold radial mixture ratio 
distribution is shown in Figure 3. The three levels of mixture ratio are 
due to a segmenting of the fuel manifold at the location of two injector 
baffle rings. Because of symmetric inlet conditions, circumferential dis- 
tributions were calculated to be within + 2 percent of nominal, and thus 

were ignored for purposes of the n r * calculation. 

c MIX 

Intra-element maldistribution data was not 
available for the M-l design configuration, therefore no intra-e lenient mixing 
loss was calculated for the injector. The mixing efficiency term accounts only 
for manifold induced element-to-element mass maldistribution. The H 2 /0 2 gas/gas 
empirically based mixing model developed in Ref. 11 was used to estimate the 
intra-element mixing .efficiency for the M-l injector. The model indicated that 
intra-element mixing losses would be insignificant because of the long (29.75 
inch) M-l chamber design! 


A simple computer program was written to sum 
streamtube performance and to evaluate the injector manifold induced mixing 
loss; by solving the following equation. 


C* 



C* 


ODE 

INJ MR 
MULTI ZONE 


x 


ODE 

ZONE MR 


(8) 



Local Row Element Mixture Ratio 
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Figure 4 indicates the results of the n c * evalua- 
tion. Calculations were made ranging from 1 to 36 streamtubes (33 iWjictor 
rows plus two baffle ring and one outer film cooling row)- to determine the 

influence of stream tube mass assignment on the n^* calculation. The calcu- 

L MIX 


lated efficiency is seen to be extremely sensitive to the selected number of 

streamtubes for flow division. The value decreases as the number of 

u MIX 


streamtubes is increased as would be expected. This sensitivity points out a 
general weakness of. the- JANNAF performance prediction methodology, that is, 
there are no standardized techniques for streamtube mass assignment in any 
of the JANNAF performance programs (i.e., CICM and DER). Since, as shown in 
Figure 3, the M-l manifold design resulted in three distinct chamber flow 
field mixture ratio zones, a three zone n r * calculation was performed. 


This result is indicated by -the dashed line in Figure 4 . The calculated value 

was equal to the case where a streamtube was assigned to each injector row. 

This. Oq* calculation technique was selected for analysis because it was 
MIX 


consistent with the physical injection zones created by the injector baffle 
design. The calculated n c * ranged from 0.976 for tests 009 and 016 to 
0.980 for the lowmixture r^Ho test number 010. 


,f. Vaporization Efficiency (iw ) 

L VAP - 

The JANNAF CICM and STC computer programs were 
utilized to calculate the injector LO^ vaporization efficiency. As explained 
in Appendix B, the .recommended program interface technique, which was utilized 
during the analysis, is to run CICM until all L0^ droplets have approached 
the chamber wet-bulb temperature. The CICM analysis was conducted by input ing 
required M-l injector/chamber geometry and selecting the program user's manual 
recommended atomization rate (C^) and vaporization rate (B^) constants shown 
in Table IV. The test vaporization calculations are summarized in Table IV. 

CICM was run to a chamber axial location of 4.10 inches (wet bulb plane deter- 
mined through one trial CICM run) from the injector face plane for all three 
tests. STC completed the calculation to the chamber throat plane axial location 
of 29.75 inches. .One zone analyses (at the test mixture ratio) were executed 
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TABLE IV 


CICM/STC VAPORIZATION CALCULATION SUMMARY 


TEST 

PROGRAM 

• . ZONES 

0/F 

009 

CICM/STC 

1 

5.46 

010 

CICM/STC 

1 

4.04 

016 

CICM/STC 

1 

5.53 

009 

CICM only 

1 

5.46 


CA . 

BA 

«AP 0X . 

' T)C *vap 

0.08 

120 

.973 

.982 

0.08 

120 

.992 

.994 

0.08 

120 

.997 

.997 . 

0.08 

120 

-v. 98 

^.99 
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for all three tests to calculate n^* - Multiple zone analyses were not 

conducted for two reasons. First, iX^tial correlation of the test 009 C* 

prediction with the test value showed excellent agreement utilizing a one 

zone o r * value. Secondly, approximately 75 percent of the injector mass 
u VAP - 

flow is contained in the outer zone (rows 16-33, See Figure 3). All of these 
rows have mixture ratio values only slightly lower than the nominal ' injector 
core mixture ratio. 

In addition to the three CICM/STC runs for each 

test, a CICM only run was conducted for test 009 to note any difference 

between a CICM/STC calculation and a complete CICM chamber calculation. The 

CICM run stopped at an axial station of 24 inches in the '29.75 inch M-l 

chamber because of a continuity check error caused by improper input of the 

chamber throat area. For this reason, the corresponding efficiency values 

shown in Table IV were deduced through extrapolation. A complete discussion 

of the CICM and STC vaporization calculation results is included in the section 

on data correlation arid analysis to follow. The CICM/STC n r * calculations 

L VAP 

were utilized in the C* efficiency predictions summarized in the next subsection. 

■ g. C* Efficiency' Prediction ( n f * ) 

L PRED 

The calculated test C* efficiencies are tabularized 
below in Table V. A discussion on correlation of the predicted and test values 
follows the next -section on test measurement uncertainties. 


TABLE V 

TEST n c * PREDICTION SUMMARY 


TEST 

'Y* 

U HL 

’Y* 
u TD 

*Y* 

KIN 

n r * 

L BL 

■Y* 

U MIX 

r* 

u VAP 

r W 

PRED 

Y* 

u TEST 

009 

1.000 

1.002 

1.000 

1.000 

0.976 

0.982 

0.960 

0.959 

010 

1.000 

1.002 

1.000 

1.000 

0.980 

0.994 

0.977 

0.964 

016 

1.000 

1.002 

1.000 

1.000 

0.976 

0.997 

0.976 

0.980 
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Test Measurement C* Uncertainties . 


The correlation of the test and predicted depend on 
the uncertainty of both values. The net correlation uncertainty is defined by 
CPIA 245 (Ref. 12) as: 




TEST 


+ 


PRED 


B TEST + B PRED 


(9) 


The precision (S) and bias values (B) depend on a knowledge of measurement 
and prediction calibrations and trends. To correlate the M-l prediction and 
test values the. following simplifications were made, because of lack of data. 


' S PRED °» B TEST “ °» B PRED ” °’ 

These assumptions indicate that the only uncertainty that can be accurately 
evaluated for the M-l analysis is the precision of the test data C* measurement. 
The following C* measurement 2 a data uncertainties were known. 

Total Weight Flow ’ +0.8% 

Chamber' Pressure + 0.4% 

- Ablative Throat Area + 0.7% 

The resultant uncertainty in test measured C* is + 1.1%. Therefore, even 
by assuming zero uncertainty in the C^ prediction and no measurement or pre- 
diction bias the agreement between measured and predicted C* (See Table V) 
is well within the accuracy of the test data, except for test 010. This 
result is discussed 'in the next section. 


C. Data Correlation and Analysis 

The results of the M-l test data correlation will be discussed 
in two parts: (1) a discussion on the results of the CICM/STC and CICM computer 
model combustion chamber energy release predictions; and (2) results of the 
correlation of the JANNAF simplified prediction procedure C* efficiencies with 
the test values. 


reproducibility OF THE 
P&aE IS POOR' 
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1 . ' Vaporization Model Results 

The CICM/STC calculated chamber pressure profiles for 

the three tests analyzed are shown in Figure 5. The analytically calculated 

profiles pass closely to the test measured static pressure values, indicating 

that the chamber energy release characteristic is being realistically modeled 

with CICM. These good correlations verified the use of the CICM/STC calculated 

chamber total pressure loss for the determination of the P value for each 

c eff 

test, as previously described in Section IV.B.l. 

As previously mentioned, a CICM only run was executed 
for test 009 to determine if the use of the simpler STC vaporization model of 
DER was compromising the accuracy of the vaporization calculation. The L0 2 
vaporization profiles for each calculational method is shown in Figure 6. 

The two calculations agreed within one to two percent over the entire 
chamber length. The CICM only calculation was extrapolated beyond the 24- 
inch axial station because of an input throat area error described in the next 
paragraph. 


The test 009 chamber pressure profiles calculated by 
CICM/STC and CICM only are compared in Figure 7. As displayed, the pressure 
profile agreement is excellent. The slight differences are attributable to the 
incorrect throat area input to CICM for the CICM only calculation. This input- 
error resulted in a continuity check error as the throat plane was approached. 

2. Correlation of Predicted and Test C* Efficiencies 


The predicted and test C* efficiencies summarized in Table 
V are graphically compared in Figure 8. Agreement was excellent for tests 009 
and 016, while there was a 1.4 percent difference (compared to a test measurement 
uncertainty of + 1.1 percent) between prediction and test for test 010. 

The test conditions are compared in Table II. The primary 
operating difference between test 016 and the nominal test 009 is an increase 
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in the injection velocity difference of from 310 to 456 ft/sec. The increase 
occurs because of the fuel density decrease associated with increasing the 
fuel inlet temperature from 97°R to 122°R. The CICM equations accurately 
predict the performance increase due to the smaller drop sizes produced by a 
higher velocity difference between the gaseous H 2 annulus and the liquid 0^ 
core- This inverse relationship is evident from the CICM mass median drop 
size correlation equation shown below. 


°j - b a 


P-5 (a -/p - ) 

P u 2 

g r 


1/2 


2/3 


( 10 ) 


The JANNAF/CICM Oq* prediction for test 010 was 1.4 
percent higher than the test value. As protrayed in Figure 8, the test 
performance for test 010 is only slightly higher than the nominal test 009 ' 
value. Referring again to Table II, it can be seen that a test 010 increase 
in fuel flowrate is offset by a higher fuel density that results in a net 
decrease in the gas to liquid jet relative gas velocity. This effect should 
lower predicted performance. However, the higher inlet density increases 
predicted performance as can be seen from equation (10). • The mass median drop 
size is inversely proportional to the fuel gas density (p g ) raised to the 2/3 
power. As described in Section B.2 of Appendix B, this CICM correlation 
dependency on the gaseous annulus density is much more severe than predicted 
by the other empirically based circular jet drop size models that has correlated 
a gas density influence. The model of Ingebo (Ref. 13) shows drop size to-be 
inversely proportional to gas. density raised to the 3/10 power. It is therefore 
suggested that CICM overpredicts the performance of test 010 because the gas 
density term is too. significant in the equation (10) drop size relationship. 


The following two observations, that resulted from the 
CICM analysis, are reiterated here to help clarify the results of the M-l data 
correlation work. (1) The M-l thrust chamber design is very similar to the J2-S 
design used to calibrate key CICM jet stripping rate and drop size constants. 
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(See Ref. 4 and J2-S sample case in CPIA 246). This is a definite reason for 
the success of the M-l performance predictions. (2) Both the M-l and J2-S ■ 
engines possess extremely long chambers that .el iminate large intra-element 
mixing losses. Therefore, the M-l predictions were not invalidated by 
assuming uniform intra-element mass distribution. 

D. Conclusions and' Recommendations 

1 . Conclusions 

The following conclusions have resulted from the JANNAF/ 

CICM analysis of the M-l thrust chamber. H ‘ 

a. The CICM model has been verified for high performing 
thrust chambers with negligible intra-element mixing losses. 

b. The CICM mass median drop size dependency on the 
gaseous annulus density is overly significant. It must be noted that changing 
the equation would most likely result in the requirement of recorrelating 

the key drop size constant, B^. 

c. The primary weakness of the CICM model is the simplified 
methodology for calculation of intra-element and inter-element (manifold induced) 
mixing losses. 


2. ' Recommendations 


The following recommendations are made based on the above 
conclusions regarding the M-l analysis. 

a. An intra-element mixing model should be developed 

for CICM. 
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b. CICM should be applied to correlation of test data 
obtained with a short chamber coaxial injector thrust chamber with a finite 
intra- element mixing loss. 

c. Reformulation and verification of the CICM mass 
median drop size correlation equation should be considered. 
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PER MASS DISTRIBUTION MODEL IMPROVEMENT 


The original objective of Task II was to provide information on 
the present prediction capabilities of the JANNAF DER and CICM computer 
programs through correlation of well documented hot fire data bases. DER 
was to be used to analyze a 6000 Ibf like doublet pair injector developed 
on the OMS engine program while CICM was to be applied to the 500,000 Ibf 
M-l engine gas/1 iquid coaxial injector. The CICM analysis was completed as . 
originally planned and is documented in Section IV of this report. 

After a- careful evaluation of the Task I DER Computer Program 
Review, it was concluded that the DER subcritical K-Prime program contains 
inadequacies in the analytical formulations that could produce invalid data 
when applied to the OMS thrust chamber analysis. It was decided that the 
originally considered funds for this task should rather be used to remove 
detected shortcomings in the model. 

Improvement of the LISP Z0M plane mass distribution methodology was 
selected as the new Task II analysis goal for three reasons. First, the 
"standardized" DER' (SDER) development program' (Contract F0 461 1 -75-C-0055) , 
conducted concurrently with this program, has concentrated on improve- 
ment of the DER vaporization modeling, but not on mass distribution and mixing 
modeling. Secondly, as discussed in Appendix A, the Z0M plane location is known 
to be a key DER input parameter which significantly influences the calculated 
chamber mixing performance efficiency. Lastly, recent empirical investigations 
have led to formulation of a model for calculation of the Z0M plane location 
on an a^ priori basis. 

The current development status of the new Z0M mass distribution 
model is summarized in the following four paragraphs that concern, respectively, 
(1) an explanation of the hypothesized model, (2) presentation of the subscale 
like doublet pair injector data base used to correlate the predictions of the 
formulated model, (3) results of data analysis and model correlation effort, 
and (4) conclusions and recommendations of this initial model development work. 
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V DER Mass Distribution Model Improvement (cont.) 

' * . ✓ 

A. Model Approach 

During a recent development effort on the Space Shuttle OMS 
engine program subscale injectors were tested to model combustion stability 
response (Ref. 14). The test combustion chamber was densely instrumented 
with static pressure transducers to allow calculation of the local combustion 
gas flowrate and velocity through the use of isentropic flow relationships. 
Bracco (Ref. 15) has also utilized this technique and developed a method for 
accurately interpreting such measurements. The availability of the OMS test 
data has resulted in empirically based mass vaporization profiles that eli- 
minate the uncertainty associated with calculating chamber gas profiles with 
DER or other available vaporization models. The uniquely accurate OMS data 
allowed calculations' of the influence of near-zone combustion gas formation 
and acceleration on liquid spray fan profiles. The results of initial cal- 
culations indicated that these effects are significant, and that further 
investigation and formulation of an analytical model was warranted. 

That .the initial model development effort described in the 
following paragraphs of this section utilized empirical energy release rate 
data as the primary model input does not imply that such data will always 
be required. The test data was used instead of analytical predictions made 
with DER because accurate vaporization profiles near the injector face were 
required. DER does not account for monopropellant burning of hydrazine 
based fuels (the OMS subscale test propellant combination was-NTO/MMH) that 
is known to significantly effect near zone energy release rates. (Monopro- 
pellant flame effects are discussed in Section C.4 of Appendix A). If the 
proposed model is ever adopted as a standard analytical procedure in DER it 
is probable that th.e DER vaporization models would have to account for mono- 
propellant burning to result in accurate mixing loss predictions. 

The originally proposed calculational technique is graphically 
portrayed in Figure 9 . The top plot in Figure S displays an empirically 
determined near zone (0-2 inches from the injector face plane) mass vaporiza- 
tion profile. Static pressure measurements included the five axial locations 
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DER Mass Distribution Model Improvement (cont.) 


shown; 0.0, 0.3, 0.6, 1.0 and 2.0 inches from the face. Isentropic flow 
relationships were used to determine the local gas flowrate, resulting in 
the plot of percent mass vaporized versus axial distance. The equations 
used to develop gas flowrate (i.e., mass vaporization) profiles from chamber 
static pressure measurements are detailed in Appendix D , taken from Ref. (15). 

The local gas flowrates were then used to calculate a chamber 
combustion gas axial velocity profile. Knowing the gas velocity profile 
allowed calculation of droplet velocity profiles through use of the standard drag 
equation and an. assumed droplet drag coefficient model. These results are shown 
in the middle plot of the figure. A mass median droplet with a constant dia- 
meter of .002 inches was assumed to have an initial velocity vector as shown. 

The droplet axial velocity increases as the combustion gas axial velocity 
increases, because of axial aerodynamic drag. The droplet radial velocity 
decreases because the combustion gas was assumed to have a radial velocity 
component of zero. 

The bottom plot on the figure shows the effect of combustion 
gas acceleration on the trajectory of a propellant droplet assumed to be 
on the outer spray fan streamline. Cold flow correlation techniques (e.g. 
the DER Z0M mass distribution method) assume a constant droplet velocity 
resulting, for the given initial droplet conditions, in the 30° spray fan 
half angle shown. If gas acceleration effects are accounted for the droplet 
trajectory, or spray fan profile, changes significantly. One of the corrected 
trajectories shown in the figure assumes the droplet is accelerated in the axial 
direction only. The other includes the effect of radial deceleration. 

The results shown in the figure indicate that, for the case 
considered, spray fan radial spreading becomes insignificant at distances 
beyond 1.8 inches of the injector face. This result implies that little 
interelement mixing would occur downstream, thus pinpointing the area for 
selection of the correct value of the DER cold flow mixing plane, Z0M. The 
initially proposed Z0M determination technique, indicated in the figure, was 
to project the corrected spray fan radial dimension back to the cold flow case. 
The hot fire spray fan mass distribution was assumed to be correctly charac- 
terized by the cold flow mass distribution at the calculated Z0M plane location. 
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A four part task was conducted to develop the proposed ZOM 
calculation technique. 

(1) Model Formulation 

The purpose of this task was to formulate the proposed 
model for calculation of a predicted hot fire ZOM plane location. The model 
was coded for the. digital computer to allow rapid reduction of the test data 
to be correlated in the data analysis subtask. ' 

(2) Data Analysis 

A test data reduction program was written to calculate 
test C* efficiencies and chamber axial gas velocity profiles. The ZOM pre- 
diction model used the gas velocity profile for each test to calculate 
the combustion corrected spray fan radial dimension and project back to the 
corresponding cold flow radial location to calculate the ZOM plane location. 

(3) Performance Data Correlation 

The DER LISP subprogram was used to predict C* mixing 

efficiency (n^* ) .as a function of the ZOM plane location. An empirically 

determined miX value was backed out for each test knowing the measured C* 

^ mix 

efficiency and analytically calculating the test vaporization efficiency. 

An empirical ZOM value was calculated for each test from the versus ZOM 

relationship calculated by LISP. Test determined ZOM values and trends were 

\ 

compared to those calculated by the analytical model. 

(4) Results and Recommendations 

The results of the initial model development effort 
were evaluated and conclusions reached. Recommendations for continuation 
of model development were formulated. 
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B. QMS Subscale Injector Experimental Data Base 

The OMS subscale injector test program documented in Ref. 14 
provides a uniquely accurate and comprehensive data base for correlation 
of predictions of the new ZOM model. Sixty-eight multi-element combustion 
tests with intensive chamber pressure profile instrumentation were used to 
infer axially distributed combustion profiles for the various injector designs. 
The OMS engine utilizes NTO/MMH propellants at a nominal chamber pressure of 
125 psia. Mixture ratio, chamber pressure, and propellant temperature varia- 
tions were tested 'to gain quantitative data on the combustion response influ- 
ences" of these engine operating variables. 

The .combustion' chamber design utilized during the testing is 
sketched in Figure 10. Pressure measurements were made- at' planes located 0., 
0.3, 0.6, 1..0, 2.0, 3.5 and 5.4 inches from the injector face plane. The 
Chamber was 8.0 -i riches in length, resulting in measured test C* efficiencies 
of 80 to 90 percent of theoretical. . The relatively low test C* efficiency 
for the coarse subscale injectors resulted in data that provided excellent 
insight into the effect of test variables on injector/chamber performance. 

Two conventional circular orifice like doublet pair (quadlet) 
and four platelet injectors were tested. A quadlet injector design was selected 
for analysis because the DER LISP subroutine contains empirical spray distri- 
bution coefficients for only conventional circular orifice element types. 

The six element, 135 Ibf thrust, quadlet injector is pictured in Figure 11. 

The fuel doublet is. positioned nearest the wall and the oxidizer- doublet is 
located inboard. A sketch of the quadlet element design is detailed in Figure- 
12. The quadlet tests selected for the -ZOM model development effort' are 
summarized in Table VI. 
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TABLE VI 

SUBSCALE QUADLET TEST SUMMARY 


Test 

0/F 

(psia) 

0 

(°F ) 

(°F) 

Oq* ( % ) 

175 

2.05 

152.5 

69 

77 

89.3 

176 

1.87 

120.7 

69 

77 

88.9 

177 

1.60 

120.8 

69 

75 

87.4 

CO 

1.59 

97.6 

71 

75 

88.9 

179 

1.69 

99.3 

72 

75 

88.8 

180 

1.71 

79.9 

73 

75 

90.2 

181 

1.66 

141.0 

74 

76 

87.0 

182 

1 .70 

142.1 

75 

75 

86.4 

183 . 

1.64 

121 .5 

73 

190 

86.5 

184 ‘ 

1.67 

123.9 

69 

184 

87.2 

185 

1'.72 

' 124.2 

69 

217 

86.6 

186 

1.72 

123.1 

141 

215 

85.1 

187 

1.68 

141.1 

137 

283 

86.3 

188 

1.73 

146.3 

130 

271 

84.3 


Statistical characterization of C* efficiency and calculation 
of empirically determined combustion gas velocity profiles for these tests 
is detailed within the following section concerning model data analysis and 
correlation. 

C. Model Data Analysis and Correlation 

1 • Quadlet Injector Test Data Reduction 

A computer program was coded to reduce the quadlet injector 
tests selected for analysis and summarized in Table VI. The primary test 
variables input to the program are injector flow areas, chamber throat area, 
propellant flowrates, temperatures, and manifold pressures and the measured 
chamber static pressures. 
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A subroutine was included in the program that contained 
parametric NTO/MMH combustion gas properties as a function of chamber pressure, 
mixture ratio, and propellant temperatures. The one dimensional equilibrium 
(ODE) properties calculated with the routine included characteristic exhaust 
velocity (C*), molecular weight, stagnation temperature, dynamic viscosity, 
and the ratio of specific heats (y). The ODE C* value was used to define test 
C* efficiency through -comparison to the -test calculated value. The remaining 
gas properties were used to compute throat effective chamber pressure and the 
test combustion gas velocity profile from the chamber axial static pressure 
measurements. . - 

A sample output case of the test data reduction program 
is displayed in Figure 13. The gas velocity profile printed as a function 
of 0.1 inch axial chamber increments was generated by applying a 2nd order 
'curve fit to the measured static' pressure data. The primary program outputs 
used. as input to the ZOM calculational model described in the next paragraph 
are the gas velocity profile and the calculated propellant injection velocities. 

2. ZOM Prediction Mode! Formulation 


The ZOM prediction model approach introduced previously was 
coded for the computer to allow rapid reduction and correlation of the sub- ' 
scale quadlet injector tests. The- function of the computer model is to inte- 
grate the basic equation for droplet acceleration based on input droplet size, 
injection velocity, ..spray fan half angle (i.e., the initial droplet trajectory) 
and the computed chamber gas velocity profile. The droplet acceleration equation 
is shown below. 


dV 

dt 





(v o - v 

D 


01 ) 


The equation was converted to allow integration with respect to the axial 
chamber distance, x. 
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dV 

dx 



( 12 ) 


The computer program utilized a special subroutine formulation of the Adams- 
Bashforth integration method. The Adams-Bashforth method is a extremely 
efficient predictor-corrector variable step size integration technique. 


The Ingebo (Ref. 16) drag coefficient correlation was 
built into the computer model coding. 

C D = 27 Re ,,' 0 - 84 ( 13 ) 

The influence of the drag coefficient assumption on the predictions of the 
ZOM model was not investigated during this initial development effort. 


The model begins execution at a designated axial plane. 

A spray droplet of mass median diameter "D is introduced at the initial plane 
with an input radial and ax;ial velocity component. The droplet acceleration 
equation is integrated and the droplet trajectory calculated versus chamber 
axial distance. The calculation is terminated at the axial plane at which 
the droplet axial velocity vector is within 0.1 percent of the total droplet 
velocity vector. That is. 



0.999. 


At this point droplet radial velocity forces that would induce inter-spray fan 
mixing are negligible. The final droplet trajectory point radial dimension 
is used to calculate the predicted ZOM value assuming a cold flow linear spray 
fan half angle consistent with the droplet initial radial and axial velocity 
components. This calculational process is explained in equation form below. 
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0 = tan 


-1 


ri 


(14) 


ai 


ZOM = tan fix r. 


where: 


ri 


ai 


(15) 



initial drop radial velocity vector 

initial drop axial velocity vector 

final droplet radial location corresponding to point 
where axial droplet velocity forces are predominant 


A sample case output of the ZOM prediction model is shown 
in Figure 14. .The droplet location can be traced through the calculated axial 
and radial' locations, X (1) and.X (2), respectively. .The calculated local 
axial arid radial velocity components at these locations are V (1) and V (2) 
respectively. The ZOM value tabulated at the final - cal culational point is 
the model predicted cold flow spray fan ZOM value for mixing efficiency pre- 
diction with the LISP sub-program of DER. 


3. Model Analysis and Data Correlation Results 


-a. Statistical Evaluation of Quadlet Injector 
Test Data 

The tests selected for analysis were subjected to 
a statistical evaluation to allow characterization of injector performance as 
a function of engine operating variables. A computer model was utilized that 
combines least squares curve fits with standard multiple regression and co- 
variance techniques. The primary test variables that were evaluated during 
the test program were chamber pressure and propellant temperature. 
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FIGURE 14. ZOM MODEL SAMPLE OUTPUT 
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DER Mass Distribution Model Improvement (cont.) 


The results of the statistical analysis are plotted 
in Figure 15. The analysis indicated that chamber pressure and fuel tempera- 
ture variances significantly influence injector C* efficiency. The statistical 
analysis resulted in the curve fit equation written below. 

n r * = 95.36 - .05279 P - .009468 T. (16) 

L TEST C T 

As shown in the figure, the equation results in decreasing C* efficiency as 
chamber pressure and' fuel temperature increase. 

The statistical analysis results indicated real 
injector operating variable influences on performance that could, hopefully, 
be modeled with the ZOM prediction model. Also, the analysis indicated that 
the quadlet injector tests comprise a high quality, repeatable data base void 
of significant measurement error or bias influences. 

b. Model Analysis 

The initial model analysis work concentrated on 
the influence of chamber pressure on test performance and' evaluation of the 
model's capability to calculate the correct absolute magnitude of the ZOM 
plane location. 


The data statistical analysis indicated a significant 
test performance efficiency sensitivity to chamber pressure. Examination of 
the test combustion gas velocity profiles calculated with the test data reduc- 
tion program gave the initial indication that the model would accurately pre- 
dict the chamber pressure influence trend. Figure 16 shows the empirically 
based gas velocity profile for a low Pc test (# 180) and a high Pc test (# 182). 
Both tests were conducted with ambient temperature propellants. As shown, the 
low Pc test resulted in a C* efficiency nearly 4 percent higher than the high 
Pc test. Interestingly, as displayed in the figure, the lower performing 
high Pc test actually possessed a significantly faster rate of near injector 
zone energy release, as reflected by the higher calculated combustion gas 
velocity. In other words, the test that exhibited high performance near the 


- 58 - 




□ 75 - 77° F Tests 
H 184- 21 5°F Tests 
S 271 - 283°F Tests 


Low and High Pc Benchmark Tests (#180£182) 
Data Points are Indicated 





Gas Velocity (ft/sec) 



FIGURE 16. CHAMBER PRESSURE INFLUENCE ON GAS VELOCITY PROFILE 


V DER Mass Distribution Model Improvement (cont.) 

injector face (the area of spray fan formation and mixing) possessed lower 
overall performance. This result appears to mechanistically agree with the 
formulated ZOM prediction model for the following reason. The higher axial 
combustion gas velocity near the face results in more rapid axial acceleration 
of spray fan droplets, thus flattening the droplets trajectory. The more 
rapid attainment of an axially directed spray fan results in a lower calculated 
value of ZOM. A lower ZOM value results in a reduction in predicted mixing 
"C* efficiency with the LISP computer model . ' ' 

► 

The initial ZOM model prediction results, for the 
cal culational technique that will be termed the baseline model, are shown in 
Figure 17. In the baseline case radial velocity deceleration is calculated 
by assuming a combustion gas velocity component of zero in the radial ..direc- 
tion, thus the droplet radial velocity component is reduced as the calculation 
proceeds axially down the chamber. The initial quadlet spray half angle was 
selected to be 40 degrees based on cold flow spray fan photographs and mass 
distribution measurements. The calculated ZOM value for each test is plotted 
versus test chamber pressure. To allow clear interpretation of the model 
predictive trend only the ambient propellant temperature test point predictions 
are .plotted. The calculated trend is opposite from that expected; that is, 
the lower performing high pressure tests have high calculated ZOM values. 

Model predictions were repeated for the same tests 
with varying calculational assumptions to ascertain the reason for the in- 
correctly calculated trend of ZOM versus chamber pressure. 'The results are 
displayed in Figure 18. The test data points were eliminated for clarity. 

: he first calculational change (Case 2 in the figure) eliminated radial decel- 
eration by assuming a constant droplet radial velocity equal to the injection 
radial component. The predicted trend of ZOM versus chamber pressure is the 
same but the absolute ZOM value is increased. ZOM increases because the con- 
stant radial velocity assumption results in a greater time to flatten the droplet 
trajectory because velocity is only changing in the axial direction. The 
second calculational change (Case 3 in the figure) was to initiate the calcu- 
lation at an axial distance of 0.4 inches from the injector face plane (using 
the empirically calculated gas velocity consistent with this location). It was 
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V DER Mass Distribution Model Improvement (cont.) 

reasoned that since jet impingement and breakup require a finite time to occur 
the droplet acceleration calculation should begin at an axial plane consistent 
with initial development of atomized droplets. This cal cul ational method 
did not significantly affect the trend or absolute magnitude of the predicted 
ZOM value. The third cal cul ational change was based on the following observa- 
tion. Although the tests at higher chamber pressure possessed higher. near 
zone gas velocities that should flatten the droplet trajectory more rapidly 
they also result in higher initial injection- vel ocities. The high initial 
radial velocity component is decelerated at a much slower rate than the axial 
component is accelerated. This results in the initial radial component pre- 
dominating the. calculation of the local droplet velocity vector as the droplet 
is marched downstream. Therefore, to verify this observation the third cal- 
culational change .(Case 4 in the figure)was to assume an initial qas and droplet 
velocity for each test of 100 ft/sec. The droplet acceleration calculation 
was initiated when the empirically calculated gas velocity exceeded 100 ft/sec. 

The resultant ZOM trend is opposite to the previous cases because the influence of 
droplet injection velocity has been eliminated. This ZOM trend is consistent 
•with the test data trend of decreasing C* efficiency with increasing chamber 
pressure. This method was varied only slightly in the final case 5 calculation 
by accounting for an influence of jet injection velocity on the atomized mass 
median drop diameter. 

These initial model ZOM predictions were evaluated 
by "backing out" test ZOM values based on actual measured test C* efficiency. 

The test correlated ZOM trend is shown as curve 6 in Figure 18 . This curve 
was developed through use of Figures 19and20. Figure 19 displays calculated 
test vaporization efficiency versus chamber pressure. The calculation was 
made with a "two-flame" modified version of the Priem L-General model (Ref. 17). 
Figure 20 shows the DER LISP subprogram predicted relationship between the ZOM 
plane location and mixing efficiency. This sensitivity curve of n^*. • versus 
ZOM was based on quadlet mass distribution coefficients developed in-^use at 
ALRC. The test "backed out" ZOM value was calculated knowing the measured test 
C* efficiency and the predicted test vaporization efficiency. 
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DER Mass Distribution Model Improvement (cont.) 


n C* . n c* / n r* 

mix u test u vap 


( 17 ) 


ZOM = f (n r * ), from Figure 20 (18) 

u mix - 

The absolute magnitude of the test ZOM value can be affected by error in- 
fluences of the vaporization calculation, the LISP calculation, and the test 
C* efficiency measurement. However, the trend of test ZOM versus chamber 
pressure accurately reflects the actual test results of increasing performance 
with increasing test chamber pressure. Returning to Figure 18, it is encour- 
aging that the model predicts ZOM values that have absolute values near those 
determined from' the test data. However, it is apparent that only the case 4 
and 5 ZOM calculational methods produce a trend approaching that deduced 
through "backing out" ZOM from the test data. 

The correlations shown in Figure 18 resulted in the 
observation that droplet injection momentum forces dominate the ZOM calculation 
in a way that overshadows the influence of higher combustion gas velocity forces 
on droplet trajectories. Opposingly, the test data trend clearly reflects a 
test variable influence that affects measured performance to a greater degree 
than injection velocity. For this reason, the possibility that Reactive Stream 
Separation (RSS or "blowapart) forces affected quadlet injector performance 
was investigated.' A discussion on RSS is included in Section C.5 of Appendix i 
A. A recently completed subscale injector test investigation (Ref. 18) indicates 
that RSS can be accurately modeled and predicted in terms of injector/chamber 
design and operating variables. The application of the Ref. 18 quadlet RSS 
model to the task test data base is shown in Figure 21. The model predicts 
that the majority of the quadlet test data is in the separated operating mode. 

The previously presented ZOM model prediction results indicate that strong 
reactive forces, such as produced by RSS, are required to result in the measured 
quadlet test data trends. Encouragingly, for modeling purposes, the test data 
indicates that RSS is a continuous process (note the linear test Hq* trend versus 
chamber pressure in Figure 15) that does not result in step function changes in 
injector performance. The same conclusion was reached in Ref. 18 . after reduction 
and correlation of several hundred tests conducted with many different injector 
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V DER Mass Distribution Model Improvement (cont. ) 

types and designs. Additionally, work has already been initiated towards the 
analytical modeling of the RSS phenomenon (Ref. 19). 

The influence of fuel temperature on the ZOM prediction 
trends was also investigated to gain further data in support of the conclusions 
reached from the chamber pressure correlation effort. Figure 22 displays the 
influence of fuel temperature on test C* efficiency for six tests conducted at 
a chamber pressure of 130 psia. Test C* shows a significant decreasing trend 
as the fuel temperature is increased. The relationship between the empirically 
based combustion gas velocity profile for a low and high fuel temperature test 
is shown in Figure -23. Again, the higher performing test (Test #176) has a 
lower rate of energy release in the injector near zone. The ZOM model predic- 
tions for the six' tests are presented in Figure 24. As before, an incorrect 
ZOM trend was produced with the baseline model. As fuel temperature increases 
the fuel density decreases resulting in increased injection velocity that again 
dominates the influence of increased axial acceleration forces. The fuel temp- 
erature correlations supports the previous results of the chamber pressure corre- 
lation effort. 

•D. Conclusions and Recommendations 
1. • Conclusions 


The following conclusions have- been reached from the initial 
a priori ZOM prediction model development effort. 

a. The OMS subscale test program (Ref .14 } has resulted 
in an excellent data base for the investigation of near-zone combustion and mixing 
phenomenon. 


b. The formulated ZOM prediction model should be tested 
with a data set that is void of significant "blowapart" forces. 

c. The gas acceleration effects ZOM model calculates ZOM 
values on the level of those required to accurately predict injector mixing 
performance. Therefore, the model most probably accurately accounts for near 
zone injection and gas acceleration momentum forces. 
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FIGURE 22. 
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FIGURE 23. FUEL TEMPERATURE INFLUENCE ON CHAMBER GAS VELOCITY PROFILE 
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DER Mass Distribution Model Improvement (cont.) 

d. Combustion reactive forces due to the mechanism termed 
"blowapart" strongly^ alter droplet inertial forces. 

e. A physically mechanistic near-zone model that will 
predict the ZOM location must account for both gas acceleration and reactive 
stream forces on droplet spray fan formation and mixing. 

2. Recommendations 

The following recommendations are made based on the above 
conclusions reached from the ZOM prediction model correlation task. 

a. The gas acceleration effects model should be further 
tested through application to a data base void of significant "blowapart" forces. 
Subscale (IK Ibf) quadlet injector data, similar to the OMS data, was developed 
on the current Improved Transtage Injector Program. This, data is at low chamber 
pressure and injection velocities and therefore is well suited for such an 
evaluation. 


b. The ZOM prediction model development effort should 
be continued with emphasis on the analytical modeling .of reactive stream forces. 
The work initiated in Ref. 19 should be evaluated for application to the ZOM 
model . • - 
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VI CONCLUSIONS 


The major conclusions from this program were: 

1. The JANNAF Performance Evaluation Methodology is being advanced 
with regard to accuracy'and applicability through conductance on this and other 
related technology programs. Such programs must continue based on resultant 
recommendations to end in valuable, standardized analytical prediction proce- 
dures. 


2. The intended predictive accuracy of the JANNAF rigorous prediction 
procedure(to within 1 percent for predicted specific impulse) is, in general, 
out of the question for a priori performance ’prediction. 

3. The generality -of the CICM program is limited due to absence 
of an intra-element mixing model. If the use of’ single element cold flow 
data to specify the intra-element mass distribution is continued, a standard 
measurement technique should be developed. A standard methodology for inter- , 
preting and inputing the data to CICM is also required. Preferably, an 
intra-el ement mixing model should be developed for CICM. 

4. The CICM analysis results verified the CICM vaporization model 
for the case where. injector intra-element mixing losses are negligible. 

5. The new ZOM mixing model development initiated during the program 
should be continued with the emphasis on the analytical modeling of reactive 
stream forces. A physically mechanistic near-zone model that will predict the 
ZOM location must account for both gas acceleration and reactive stream forces 
on droplet spray fan formation and mixing. 
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This appendix details the results of the DER computer program review. 

The review was accomplished in four subtasks. (1) Identification and Correction 
of Coding Errors, (2) Addition of Diagnostic Comment Cards and Print-Out 
Statements, (3) Identification of Inadequate Formulations and Model Technical 
Shortcomings, and (4) Review of the JANNAF Performance Prediction Procedures 
(CPIA 246) with regard to use of DER. A complete summary of the reconmendations 

resulting from the' review are included in Section III. A. of this report. 

A. Identification and Correction of Coding Errors 

A new ‘"standardized" DER program is currently undergoing final 
development. General release is planned for the fall of 1976. This effort 
will result in a code considerably changed from the subcritical K-Prime version 
reviewed during. Task I of the Injection Processes (IP) Program. For this reason 

no attempt was made to verify every formulation in the DER code. The results 

of the coding review are presented below. The majority of the comments do hot 
concern errors, as such, but points that should be brought to the attention 
of the DER user to generate increased understanding of- program limitations. 

1 . LISP Unsymmetrical Pie Section Input Problem 

The following must be true to- result in an accurate 
total propellant flowrate integration calculation at the LISP collection 
plane (ZOM). 

(a) For an injector slice of 8 degrees the slice must 
contain exactly (0/360 x 100) percent of the total number of injector elements. 
This requirement is sometimes difficult to achieve for fine patterns. If the 
above stipulation is met LISP will execute properly. However, the total flow- 
rates used in STC will be in error unless the following is also true: 

(b) 9 must be a integer divisor of 360 degrees. That 
is., a 0 value of 40 degrees will work, but a 0 value of 39 degrees will cause 
an error in the STC total flowrate. 
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In the case of an unsymmetrical injector it is 
sometimes impossible to satisfy both points (a) and (b). An improved technique 
would be to adjust the collected pie section mass flowrate to 0/360 of the total 
injected flow of each propellant. 

2. STC Mesh Point Dimensional Limits . 

LISP will execute properly if the total number of mesh 
points (NRML, constant’ radius lines x NTHML , constant 0 lines) is equal to 
400 or less. Any combination of NRML and NTHML will work. However, dimensional 
arrays in STC' require that NTHML < 20 and NRWALL <_ 20 (number of NRML to wall). 
Otherwise, the STAPE and SCRMBL routines will compute inaccurate streamtube 
flowrates-. This STC limitation is not noted in the DER user's manual. It should 
be noted in the user's manual, or preferably removed to allow any NTHML-NRWALL 
combination in STC. 

3. Drop Size Equation Inconsistencies 

Examination of the drop size routine DSIZE indicated that 

two drop size equations differ from the equations given 'in DER documentation. 

The equations were inconsistent for (1) the center orifice of a Triplet or 

Pentad (4-on-l )' el ement, and (2) the contraction ratio adjustment factor for 

secondary atomization of like doublet elements. The differences should be 
✓ 

resolved, but basic questions concerning the validity of DER drop size pre- 
diction equations are a more important issue. The drop size equations are 
thoroughly evaluated in Section C.l. of this appendix. 

4. ' Triplet and Pentad Collection Plane Rotation Error 

An error exists in the LISP calculation of the Z0M mesh 
point mass fluxes for triplet and pentad elements. Inline triplet and pentad 
elements are symmetrical about both the face plane X and Y axis. The LISP 
subroutine SC0EF calculates a rotation of the Z0M collection plane around the 
normal X axis based on the relative fuel and oxidizer element momentum. For 
a regular symmetrical triplet or pentad the resultant spray fan will always 
be normal to the chamber longitudinal (Z) axis. The skewing of the collection 



plane calculated by LISP results in incorrectly computed ZOM nodal point mass 
fluxes. This error can be eliminated by setting the variable ALFMOM equal 
to 0 in the triplet section of the SCOEF code. 

5! . LISP Infinite Baffle Height Assumption 

When calculating mass distributions for injectors with 
baffles, the LISP subroutine BNDY assumes infinite baffle height. This 
technique results- in large accumulation of mass at the baffle boundary. The 
ZOM mass distribution should be a function of baffle height. This limitation 
should be noted in the DER user's manual. 

B. ' Addition of Diagnostic Comment Cards and Printout Statements 

Error message requirements identified from Section A of this 
appendix and previous DER analyses are listed below. They are confined to 
the main LISP and STC routines. 

1 . LISP Error Messages 

The main inconvenience for the user of LISP is that input 
errors that are detected do not have accompanying error messages that speci- 
fically pinpoint the problem. 

(a) An error message is required to explain inconsistency 
between the program - NTHML, NTHL, and NTHR inputs. 

(b) An error message i-s required to identify input that 
assigns excessively high values to the NMESH, NEL, and NLSPEC program variables. 
NMESH is -the total number of LISP nodal points, i.e., the product of the 
inputs NRML and NTHML. NEL is the total number of injector elements and NLSPEC 
the number of .different element specifications input. 

(c) An error message is required to specifically state 
the program failure associated with improper input of the Type 8 injector spray 
coefficients. 
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2 . 


STC Error Messages 


, (a) The limitation on the STC radial and angular mesh 

lines has been alluded to (NTHML and NRWALL <_20). If this limitation is not 
removed an error message should be included in STC to identify the problem. 

C. Identification of Inadequate Formulations and Model Technical 

Shortcomings 

Several inadequate formulations and shortcomings have been' identi- 
fied which limit; the current predictive capability of DER. The purpose of 
this subtask of the DER review was to identify model problem areas and, if 
possible, to propose alternate approaches for improvement. Incorporation of 
the improvements is, for the most part, beyond the scope of this program, but 
their identification will provide a basis for future DER work. The model criti- 
que is summarized in the following six separate sub-sections concerned with, 
respectively; drop size predictions, ZOM plane selection, a DER vaporization 
model sensitivity study, near-zone combustion and monopropellant flame consid- 
erations, combustion gas acceleration and reactive stream separation effects 
on cold flow mass distributions, and the need for a turbulent mixing model in 
the STC subprogram. A seventh and final section contains a proposed approach 
for combining ZOM plane, combustion gas acceleration, reactive stream separation, 
and turbulent mixing considerations inot a physically realistic mass distri- 
bution and mixing model for DER. 

1 . Drop Size Prediction 

The inaccuracy of the LISP drop size predictions has been 
a major DER shortcoming identified by two studies that used DER to analyze 
engine performance data (Refs. 20 and 21). The DER drop size correlations 
were examined for coding accuracy and for reference DER equation predictions 
were compared to those made with the empirically based drop size model of 
Priem (Ref. 17). The drop size prediction equations for LISP element types 
1-5 (unlike doublet, like doublet, like-doublet pair, triplet, pentad (4-on-l), 
respectively) were included in the study. Several DER publications indicate- 
that Ref. 22 includes a section showing all current DERdrop size equations, 
but examination of this report uncovered no such write-up. 
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a. Unlike Doublet Drop Size Equations 

The DER unlike drop size equations are presented in 
Ref. 23. They are shown below. 

Unli-ke Doublet (larger Diameter Orifice - ) 

.38 



Unlike Doublet (Smaller Diameter Orifice) 




2.29 


O’ 27 D 


.023 


0 ££_ 


U - 74 U 
D D,opp. 
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D^ is the mass median drop diameter (inches) determined from molten wax 
atomization experiments (Ref. 24). D^ is corrected in theDER'code to account 
for secondary (aerodynamic) break-up with the following general equation 
(developed in Ref. 25). 



_T_ 

J A, 


- + B 


(A-3) 


For unlike doublets DER uses values of .8 and 250 for and B, respectively. 


In Figure A-l drop size predictions made with the 
DER like doublet equations are shown. The predictions shown were made for 
equal .orifice diameters (D-| a e = D $ma -|-j) for orifices from .01 inch to .1 inch 
in diameter (a typical orifice design range). There are four apparent anomalies 
with the equations. The first peculiarity is that for the same jet size and 
velocity the large orifice and small orifice D^ equations predict drop sizes 
one order magnitude different. Secondly, whereas the secondary break-up 
■correction reduces the predicted small orifice drop size significantly, the 
correction actually increases the large orifice drop size prediction. The 
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FIGURE. A-l. UNLIKE DOUBLET DROP SIZES 
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third anomaly is that the large orifice equation correlated from the hot wax 
experiments shows that-jet diameter has no influence on the mass median drop 
size. Included in the figure, for reference, are predicted drop size trends 
made with the Priem model for liquid heptane drops. The Priem data shows a 
significant effect of jet diameter on drop size. Last, the unlike- doublet 
equations do not allow for propellant property effects on the predicted drop 
size. Priem ‘"s -model results in drop size being proportional to the propellant 


properties grouping .shown. 


D 




(A-4) 


Typical' atomization models found- in the literature, such as those of Weiss 
and Worsham (Ref. 26) ; Ingebo and Foster (Ref. 27), and Nukiyama and Tanasawa 
(Ref. 28) also allow for the influence of liquid properties on the atomized 
drop diameter. 


b. Triplet and Pentad Drop Size. Equations 

The drop size equations for triplets and pentads, 
taken from- Ref. 23, are shown below. 

Triplet and 4-on-l (Center Orifice) (Ref. 23) 
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Triplet dnd 4-on-l (Outer Impinging Streams) 
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Examination of the DER code revealed a different equation for the center orifice. 
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Triplet and 4-on-l (Center Orifice) (PER Code) 


, / 

For triplets and pentads DER uses and B values of 0.03 and 310, respectively, 
in the secondary breakup equation (Eq. A-3). The variation in and B con- 
stants for unlike doublets and triplets imply that equal drop sizes produced 
with different injector element types result in different secondary breakup • 
characteristics. There would appear to be no physical basis for such an effect. 

. ' Predictions made with the DER triplet equations, 

for a range in orifice diameter from .01 to .1 inches, are shown in Figure A-2. 
The (hot wax) equations show realistic trends, but- these results are 
obliterated by the secondary breakup correction. As an example, the outer 
orifice D^ predictions range from about .006 to .013 inches. When inserted 
in Eq. A-3 the predicted drop size range is from .00317 inches to .00320 inches. 
It is apparent from this result that the DER triplet equation is effectively 
a constant (1/310), since the first term in the Eq. A-3 -denominator will always 
be small compared to the constant second term. Additionally, as was the case 
for the uni ike. doubl et, the triplet equations do not allow for the effect of 
propellant properties on the predicted mass median drop size. 

c. Like Doublet (single or pair) Drop Size Equations 
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Ref. 2 shows the following equation to account for 
the combined effects of hydraulic and secondary breakup for like doublet 
elements. ' • 


D = 1.524 




(A— 8) 


The formula for the contraction ratio function (e c ) varies between Ref. 2 and 
the DER code. 
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, x 5 (c-1) 

f ( e ) = 9 ( Re f. 2) 


(A-9) 
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(DER Code) 
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Predicted drop sizes for these equations are shown in Figure A-3. As shown, 
the Priem like doublet correlation shows considerably more sensitivity to jet 
orifice diameter than the DER correlations. The variable C pr .is a propellant 
properties term, that is only included for the secondary breakup part of the 
drop size equation. 


2. ZOM Plane Selection 


ZOM is the interface plane (measured from the injector 
face plane) between the LISP and STC subprograms that comprise .DER. ZOM 
is the plane at which the LISP mass distribution is calculated. The mixing 
limited performance loss is directly dependent on the ZOM plane calculated 
mass distribution because STC does not account for turbulent mixing or com- 
bustion effects. Two analytical studies (Ref. 20 .and 21) . have determined that 
DER predictions are quite sensitive to selection of a value for ZOM. ' Ref. 2 
provides only the following two guidelines- to selection of ZOM: (1) the collec- 
tion plane should be far enough downstream to account for substantial spray 
spreading and wall impingement, and (2) because LISP does not account for inter- 
element spray interaction, spray mixing and impingement effects can be over- 
predicted if ZOM is too far downstream. It is apparent that an a priori method 
for selection of ZOM does not currently exist. 


LISP calculates the mass distribution from an injector 
element with the foil owing general equation. 



(f) ( A— 1 1 ) 

-b (z) 2 


1 is the distance from the element impingement point (H), thus ZOM is the sum 
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FIGURE A-3. LIKE DOUBLET DROP SIZES 
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of Z and H. The mass distribution coefficients (C-,, C-^, C 3 , C 4 , Cg, C g , a 5 b) 
are evaluated empirically from single element, cold flow data generated with 
propellant simulants. The equation assumes that the element spray can be charac 
terized as point source flow and that the spray coefficients for the equation 
are constant, independent of the Z distance. This assumption has not been 
verified experimentally. 

The LISP mass distribution equation results in a linear 
half angle spray fan spreading characterization as shown in Figure A-4. 
Inter-spray mixing increases as.ZOM is increased because of spray fan overlap. 
LISP does not account for any spray fan interaction thus adjacent spray fan 
mass distributions are simply superimposed on one another. Since no consis- 
tently accurate a prior i method for selection of the ZOM value exists, attain- 
ment of an accurate value usually depends on an iterative process utilizing’ 
available hot-fire performance data. 


DER . 
COLD FLOW 



HOT FIRE 



SPRAY FAN RADIAL 
FORCES INSIGNIFICANT 


FIGURE A-4. COMBUSTION EFFECTS ON COLD FLOW SPRAY FAN PROFILE 
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Finite mixing losses are experienced in hot firings because 
combustion gas acceleration and reactive stream separation (RSS) forces (if 
any) combine to impede inter spray fan mixing. It therefore seems reasonable 
that selection of the ZOM plane should account for the influence of combustion 
on the chamber spray fan mass distribution. As combustion gas is formed and 
accelerated significant axial droplet drag forces are generated. These forces 
result in an effective bending of the cold flow spray fan, as shown in Figure 
A-4. Eventually, the axial spray fan drag forces dominate any radial droplet 
velocity forces and inter spray fan mixing stops. During Task II of the IP 
program a methodology for accounting for combustion effects and for a priori 
selection of the ZOM plane location was developed. The current status of this 
model is detailed in Section V of this report. 


Recent work (Ref. 18) indicates that, in addition to normal 
combustion gas acceleration effects. Reactive Stream Separation (RSS or "blow- 
apart") can also significantly affect chamber mixing performance. RSS has been 
successfully correlated as a function of the injector element design and 
operating point. ' Section B.5 of this appendix expands the RSS topic and suggests 
ways of incorporating RSS' modeling in the LISP mass distribution formulation. 

A proposed overall development plan for general improvement and update of the 
LISP ZOM plane mixing technique is included in Section C.7of this appendix.' 

A final limitation of the ZOM plane methodology is the 
relatively narrow parametric range .over which spray distribution coefficients 
have been catelogued in LISP. Currently, if the design point to be analyzed 
is not within the spray coefficient range for the element type in question 
the spray coefficients for the nearest available design point are selected. 

No technique currently exists for extrapolation-of the spray coefficients. 

The current parametric range for the LISP coefficients for the five primary 
liquid/liquid element types is shown below in Table A-I. Expansion of the 
spray coefficient matrix is required if the ZOM technique is retained in DER. 
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TABLE A-I 


LISP SPRAY COEFFICIENTS PARAMETRIC RANGE 


(1) 

Element 

Orif. Dia. 
(in.) 

Imping. 
Angle 
(deg) . 

Momentum 

Ratio 

(f/o) 

Unlike Doublet 

0.020-0.079 

45-70 

9.42-1.0 ' 

(2) 

Like Doublet 

0.020-0.079 ■ 

45-70 

1.0 

(3) 

Like Doublet Pair 

0.020-0.028 

60 Like, 
40 Uni ike 

Not Correlated 

(4) 

Triplet ’ 

0.085-0.067 Outer 
0.043-0.067 Inner 

70 

0.3-8. 0 

(5) 

4-on-l (Pentad) 

0.21-0.47 Inner - 
0.1-0.22 Outer 

60 

0.2-1.25 ’ 


3. PER Vaporization Sensitivity Study 

The previous subsections have detailed the three most 
critical DER analysis input parameters'; the mass median drop diaemter, the LISP 
.spray distribution coefficients, and the LISP/STC interface plane, ZOM. The 
specification of these LISP inputs controls, in large part, the accuracy of 
the DER calculation. The most important function of the Streamtube Combustion 
(STC) subprogram is to compute propellant vaporization to the chamber throat 
plane, after correct input is established and STC has segregated the LISP 
calculated mass distribution into a finite number of axisymmetric streamtubes. 
The third subtask of the DER review was a vaporization sensitivity study 
conducted to determine the influence of engine design' and operating variables 
on the DER vaporization calculation. DER predictions were compared to similar 
calculations made with the simplified Priem L-General model (Ref. 17) for 
reference. The Priem L-General model is an empirical correlation of an analy- 
tical vaporization model that accounts for droplet heating. The L-General 
model accounts for the effect of chamber length, contraction ratio, chamber 
pressure, injection velocity, drop size, initial propellant temperature and 
propellant properties on vaporization rate. 
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The vaporization sensitivity study was conducted by running 
single stream calculations with STC. The NTO/MMH propellant combination was 
selected because of experience obtained during the Space Shuttle OMS Engine 
program. Oxidizer and fuel mass median drop diameter, initial velocity, and 
total flowrates were input to STC. A log normal drop size distribution 
(a = 2.3), with five size groups for each propellant, was used throughout the 
study, except during the phase of the study which evaluated the effect of the 
specified drop'distribution on the vaporization efficiency prediction. The 
independent design and program input variables evaluated during the study 
included chamber length, mass median drop diameter, chamber pressure, droplet initial 
(injection) .velocity ,ehamber contraction ratio, propellant inlet temperature, 
and the propellant drop size distribution. The study nominal calculation 
point and the parametric range of the independent variables is detailed in 
Table A-II. 


TABLE A-II 

DER VAPORIZATION SENSITIVITY STUDY 
VARIABLE RANGES 


Variable 

Nominal 

Value 

Range 

Chanber Length, L 1 - (in) 

6 

6-14 

Mass Median Drop Diameter, D (in) . 

.002 

.001 -.004 

Chamber Pressure,- Pc (psia) 

120 

• 60-240 

Injection Velocity, V. (ft/sec) 

J 

65 

65-200 

Contraction Ratio, 

1.9 

1.9-5 

Propellant Temp., Tp (°F) 

70 

40-130 


In addition to the DER vaporization sensitivity study, 
the recent comprehensive combustion model evaluation conducted by Bracco (Ref. 
7) was reviewed. This study resulted in three conclusions that are relevant 
to the DER review. - First, the Priem vaporization model, used as the reference 
technique in the sensitivity study, was judged to be capable of accurate corre- 
lation of empirical ethanol mass vaporization profiles. Secondly, the Priem 
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technique was most accurate if the droplet drag coefficient was assigned a 
low but finite value (0 < < 24/Re). A major conclusion by Bracco was that 

low drag coefficients give better results than higher drag equations, such as 
the Rabin equations (See Section C.3.d. of this appendix) used in DER. In 
support of lower drag coefficients he cites Eisenklam's suggestion (Ref. 29) 
that burning droplets actually have lower drag coefficients than solid spheres. 
The final significant conclusion of Braccos work is that a drop size distri- 
bution function is not necessary to accurately reproduce the steady combustion 
profile, although it does tend to improve the results. 

The results of the vaporization study are presented below 
for each design- variable shown in Table A-II. 

- a. Chamber Length 

Propellant vaporization was calculated at three 
chamber lengths (6, f 10, and 14 inches) with DER and the Priem L-General model.' 
Figure A-5 shows that both models predict similar trends of propellant vapori- 
zation versus chamber length. The calculated absolute levels differ somewhat 
for the MMH vaporization characteristic. This difference is believed to be 
related to the fact that DER does not allow for the finite time required to 
heat the liquid propellant to the chamber "wet bulb" condition. This omission 
is described more fully in the section dealing with propellant inlet temperature 
considerations. The trend agreement versus chamber length for the Priem and DER 
models suggests that both model the gas dynamic, droplet ballistic, and steady 
state heat and mass transfer processes similarly. This result is clouded some- 
what by the chamber pressure sensitivity result presented in the next paragraph. 

b. Chamber Pressure 

Propellant vaporization was calculated for chamber 
pressures ranging from approximately 60 to 240 psia. The results of this phase 
of the sensitivity study are shown in Figure A-6. The figure indicates a 
significant difference exists between the predicted effect of chamber pressure 
on propellant vaporization for the DER and Priem models. The DER model indi- 
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FIGURE A-5. CHAMBER LENGTH EFFECT ON VAPORIZATION 
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% Vaporized At Throat Plane 



FIGURE A- 6; CHAMBER PRESSURE EFFECT ON VAPORIZATION 
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cates no sensitivity to chamber pressure, while the Preim model shows a signi- 
ficant influence*. The DER and Priem vaporization formulations were briefly 
investigated to determine the reason for the chamber pressure influence differ- 
ence. The DER K-Prime model is based on the early work of El Wakil and others 
(Ref. 30). An admitted weakness of this model was the need for a correction 
factor accounting for unidirectional, as opposed -to equimolal, droplet vapor 
diffusion in the mass transfer equations. The Priem model accomplishes the 
transformation through the following equation. 


where: 


a 



unidirectional 



) 


equimolal 


x 




(A-12) 


The mass transfer rate equation without the unidirectional diffusion correction 
is written as: 


w 


K 


a,s 


(A-13) 


Including the a term yields the following, which 
results in a term showing a directly proportional relationship between the mass 
transfer rate and .the local static pressure. 


W = A K p In 
s s 



(A-14) 


It is suggested that possibly the DER K-Prime model 
does not accurately account for the effect of chamber pressure on the vapori- 
zation rate because of this omission. Complete resolution of this question was 
beyond the scope of the current work. 


*Priem correlated an effective chamber length 


p 0.66 

c 


^ L gen^ 


as being proportional to 
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c. Mass Median Droplet Diameter 

Propellant vaporization was calculated for mass 
median driplet diameters of .001 , .002, .003, and .00-1 inches. Relatively 
small droplets were selected because state-of-the-art injector designs are 
attaining performance consistent with such drop sizes. The results of the 
drop diameter study are plotted in Figure A-7. The trends for the DER and 
Priem models are nearly identical. The figure indicates that both formula- 
tions account for the influence of drop diameter on mass transfer, heat 
transfer, and droplet ballistics. 

d. Droplet Initial (Injection) Velocity 

Propellant vaporization was calculated for droplet 
initial velocities of 65, 100, and 200 ft/sec. The results are shown in 
Figure A-8. -The Priem model indicates a much greater sensitivity to initial 
velocity than does DER. Two differences in the droplet ballistic equations 
of the two models have been discovered. 

First, the DER integration technique (a simple 
step-by-step Euler approach) often predicts a droplet downstream station velo- 
city greater than the downstream gas velocity. When this occurs DER sets the 
downstream station droplet velocity equal to the upstream station droplet 
velocity, resulting in an unrealistically long droplet chamber residence ..time 
and increased propellant vaporization. This result is indicated in Figure 
A-8 by the high DER vaporization efficiencies predicted for high droplet 
initial velocities. 

The second droplet ballistics difference between 
the two models is in the formulation of the droplet drag coefficient. The 
Priem model employs the empirical correlation developed by Ingebo (Ref. 16). 

C d = 27 Re D " 0 - 84 (A-15) 

DER uses a variation of the Ingebo result, developed by Rabin (Ref. 31) above 
Reynolds numbers of 80. 
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FIGURE A -7 .MASS MEDIAN DROP DIAMETER EFFECT ON VAPORIZATION 




C d = 27 Re D 


(A-16) 


-0.84 


Re^ < 80 


C d = 0.271 Re D ' 217 80 £ Re Q < 10 4 (A-17) 

C d = 2 Re D > 10 4 (A-18) 

Considerable discussion occurs in the literature over the validity of the 
two models. As introduced previously, Bracco (Ref. 7) has concluded that 
high C D ‘s (such as the Rabin equations) give erroneous mass vaporization 
profile results. Also, an excellent synopsis of relatively current thought 
is contained in Ref. 32. 


The effect of substituting the Ingebo correlation 
for the Rabin correlation on DER vaporization predictions is shown in Figures 
A-9 and A-10. In Figure A-9 the predicted vaporization efficiencies are 
plotted versus injection velocity for both drag coefficient correlations. The 
Ingebo equation increases the absolute vaporization efficiency level and results 
in a slope more nearer the Priem model result. In Figure A-10 DER predictions 
for both models are plotted versus chamber length. The slope of the predictions 
are nearly equivalent, while the Ingebo equation results in a significantly 
higher rate of propellant vaporization. 

e. Chamber Contraction Ratio 

Propellant vaporization was calculated for chamber 
contraction ratios of 1.9, 3, and 5:1 for a constant chamber length of 10 
inches. The results are shown in Figure A-ll. The differences in the slooe 
of the predicted DER and Priem model results are likely to be attributable 
to the two droplet ballistic model inconsistencies commented on in the previous 
paragraphs. 


f. Propellant Temperature 

The DER K-Prime model does not allow for the finite 
time required for a droplet to heat from its injection temperature to the "wet 
bulb" state. DER has no prescribed method for accounting, for the effect of the 
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FIGURE A-9. INJECTION VELOCITY EFFECT ON VAPORIZATION FOR DIFFERENT 
DROPLET DRAG COEFFICIENTS 
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FIGURE A-ll. CONTRACTION RATIO EFFECT ON VAPORIZATION 
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propellant inlet temperature on vaporization efficiency because of this model 
simplification. Priem used his time based mass and heat transfer equations 
to compare the time to reach the "wet bulb" condition to the time to vaporize 
99 percent of the mass for several different propellants. The results are shown 
below. 


TABLE A- III 

PROPELLANT HEAT-UP TIME CHARACTERISTICS 

Length to Wet Bulb 
Propellant Length to Vaporize 


(Pc = 300 psia) 


Heptane 1/5 
-Hydrazine 1/12 
Ammonia 1/16 
Oxygen ' 1/10 
Fluorine 1/10 


It is apparent that this initial unsteady state can be significant when 
accounting for a complete droplet time history. . The time to reach the "wet 
bulb" condition is primarily dependent on the droplet diameter and the initial 
propellant temperature. 


Figure A-12 shows the Priem model predicted effect 
of propellant temperature on vaporization for the study baseline calculation 
point. An attempt to allow for this effect with DER was made by adjusting the 
MMH input latent heat of vaporization. This result, shown in the figure, is 
considered to be unsatisfactory. The most physically correct way to solve 
this shortcoming of DER would be to adopt a time dependent vaporization model. 


g. ■ Drop Size Distribution 

The importance of the droplet size distribution 
on propellant vaporization has long been recognized. The DER user should 
realize that the DER builtin drop size distribution may not be physically 
accurate for his particular injector design*. Figure A- 1 3 shows various drop 
size distributions found in the literature. • The builtin DER drop size dis- 

*The DER user can override the builtin distribution through input. 
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FIGURE A-l 2. PROPELLANT TEMPERATURE EFFECT ON VAPORIZATION 
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FIGURE A-13. DROPLET SIZE DISTRIBUTIONS 



tribution is based on the Rocketd yne molten wax experimental results (Ref. 
24). It is quite similar to the well-known Nukiyama-Tanasawa empirically 
determined cold flow distribution which is also shown in the figure. Priem 
correlated his hot fire data with log-normal distribution equations. 


dR 

dD 


- 1/2 
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Priem determined that doublets and triplet hot test sprays were best described 
with standard deviations of 2.3 and 3.6, respectively. These two distributions 
are also shown' in Figure A-13. The effect of the drop size distribution on 
propellant vaporization is shown in Figure A-14. The significance of the dis- 
tribution on the predicted level of vaporization is evident. The DER and' 
Nukiyama-Tanasawa distributions are cold flow (i.e., gas static) distributions. 
The implication of the Priem correlations is that a dynamic (i.e., accelerating) 
gas environment affects the atomization process thus resulting in a different 
hot test distribution. 


4. Near-Zone Combustion and Monopropellant Flame Considerations 

Monopropellant decomposition burning for hydrazine based 
fuels has been .verified by several investigators. Decomposition burning 
results in higher energy release than the bipropellant reaction in the injector 
face near zone. A recommendation has been made recently (Ref. 33) not to in- 
clude a decomposition flame model in the "standardized" DER program. Two reasons 
for this recommendation were cited: (1) "Two Flame" effects are only important 
close to the injector face and do not significantly affect the vaporized pro- 
pellant mass fraction at the chamber throat plane, and (2) the combustion cham- 
ber near-zone flow field can not be well defined analytically. 

Monopropellant burning does not significantly affect 
performance for most thrust chamber designs. However, a valid DER performance 
model would extend itself naturally to combustion stability and chamber com- 
patibility analytical modeling. Accurate stability and compatibility pre- 
dictions hinge on a realistic representation of the near zone flow field, where 
monopropell ant effects are significant. Also, the new ZOM mass distribution 
model described in Section V of this report depends' on an accurate vaporization 
rate calculation near the injector face. 
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Recent work on the Space Shuttle OMS. engine program (Ref. 

14) indicates that the near combustion zone can be well modeled. This con- 
clusion was reached through correlation of analytical predictions made with 
the Priem model with empirically measured energy release rate data. Figure 
A-15 shows actual -and predicted energy release data for a. subscale platelet 
V-doublet injector. The analytical prediction was made with a "two flame" 
transformation to the simplified Priem L-general model. .The OMS engine 
utilizes the Nitrogen Tetroxide (N^O^j/Monomethyl Hydrazine (MMH) propellant 
combination, thus justifying a two flame correction to the MMH vaporization 
calculation. Excellent agreement between prediction and test was obtained 
for the injector near zone, from 0-2 inches from the injector face, as' shown 
in the plot. The correlation continued to be valid to the throat plane, 8 
inches from the injector face plane (not shown).- The predicted vaporized 
propellant (gas) mixture ratio profile for the two flame model prediction is 
also, shown in the figure, along with the mixture ratio profile made with the 
Priem model without the "two flame" correction. The difference is quite signi- 
ficant, indicating the near zone is not modeled .correctly unless decomposition 
burning is accounted for. The local gas composition’ is directly related to 
the local gas mixture ratio, indicating the importance, of the "two flame" 
correction to chamber stability and compatibility modeling. 

The semi -empirical technique developed to convert measured 
chamber axial static pressure profiles to injector energy release rate profiles 
is graphically illustrated in Figure A-16. The measured static pressure at any 
chamber axial location is used to calculate the local gas velocity and flowrate 
through isentropic -relationships. The cumulative energy release to the same 
plane is calculated knowing the percentage of the total propellant burned. The 
energy release rate is predicted by taking the first derivative of the cumulative 
energy release profile. It has been determined that the resulting injector 
energy release rate characteristic can be directly related' to the injector 
combustion stability characteristic. Additionally, data obtained in this manner 
for OMS multi-element subscale like-doublet pair injectors are being used to 
develop and verify the new ZOM mass distribution model described in Section V 
of this report. 
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FIGURE A-16. SEMI-EMPIRICAL NEAR ZONE COMBUSTION MODEL 


The OMS test results indicate a "two flame" model is 
required to accurately model the injector near zone for injectors employing 
hydrazine based fuels. Examples of available models are Refs. 34 through 37. 

5. Combustion Gas Acceration and Reactive Stream 

Separation Effects on Cold Flow Mass Distribution 

The LISP mass distribution calculation technique does 
not allow for the influence of combustion on the elemental cold flow mixing 
characteristics. The formation and acceleration of combustion gas affects 
spray distribution for all liquid propellant injectors. Additionally, dependent- 
on the injector operating point, the phenomenon termed Reactive Stream Separation 
(RSS or "blowapart") can also alter the hot fire case mass distribution from 
one measured under cold flow conditions. 

A -model has been proposed to account for the influence of 
combustion gas acceleration on the calculated chamber mass distribution. Also, 
inherent in the proposed model is a technique for a prior i estimation-of the 
ZOM plane location for mass distribution characterization. The model is des- . 
cribed, along with a report on the current status of a model .verification 
effort, in Section V of this report. 

The effect of RSS on injector performance can be significant. 
Figure A-17 indicates the influence of RSS on the energy release efficiency (ERE) 
for platelet "splash plate" injectors as a function of engine chamber pressure. 
For one particular injector tested RSS decreased injector ERE approximately 
10 percent for a chamber pressure range of from 50 to 110 psia. 

A recently completed investigation (Ref. 18) indicates 
that RSS phenomenon can be accurately modeled and predicted in terms of injector/ 
chamber design and operating variables. Single element unlike doublet, F-O-F 
triplet, and platelet injectors were tested. The mode of operation of a parti- 
cular test (i.e., mixed, mixed-separated , or separated) was determined through 
filming of the combustion with a high speed motion picture camera. Results 
were correlated with the test design and operating point to result in a mechan- 
istic RSS model. 
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An example RSS correlation is shown in Figure A-18. The 
occurrence of RSS is plotted for a single unlike doublet element as a function 
of chamber pressure and a parameter containing Weber Number, Reynolds Number, 
and a propellant temperature influence term. The Weber number relates the jet 
aerodynamic drag force to the liquid surface tension force. The results indi- 
cate that the occurrence of RSS can be accurately predicted if the injector 
operating point is well characterized. 

Inclusion of RSS results; such as shown in Figure A-18, 
in LISP could provide a significant improvement- in the DER performance pre- 
diction capabilities under conditions where RSS occurs. There appears to be 
two types of RSS models that could be incorporated into LISP. (1) A simple 
"warning" model that' could tell the user that the occurrence of RSS is 
predicted for his' input injector design and operating point, and' (2) a model 
that would predict the occurrence of RSS and would adjust the LISPZOM plane 
mass distribution based on actual hot test 'mass distribution results. 

. The first model would serve only to provide the designer 
with more information. It would not provide a quantitative estimate of the 
effect of RSS on injector performance. 

The second model approach represents the most quantitatively 
accurate approach to RSS performance modeling. Testing would have to be per- 
formed .that would result in measured hot fire mixed and separated mass distri- 
butions. LISP mass distribution coefficients would be developed to account 
for the influence of RSS. Such measurements have been performed for a gas/gas 
swirl coaxial element at ambient chamber pressure (Ref. 38). A double-walled, 
hot hydrogen cooled probe was used to withdraw the gas sample. The gas sample 
mass composition was measured by a mass spectrometer. Measured mixed and 
separated mass distribution coefficients could be correlated as a function of 
the injector design and operating point, to result in an experimental ly verified 
model to be readily inserted into DER. 

An additional RSS model approach exists that is adaptable 
to the JANNAF simplified performance prediction methodology, but not necessarily 
compatible with the DER program computational techniques. This model would 
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Figure A-18. Correlation of. RSS Test Data 





predict the occurrence of RSS and adjust the performance prediction through an 
empirically based performance correlation technique. The model would require 
development of a correlation between a RSS prediction parameter (such as the 
abscissa of Figure A-18) and actual injector performance. A possible technique 
would be to relate mixed and separated mass and mixture ratio distributions 
through stream tube performance relationships. 

6. Turbulent Mixing Model 

The LISP cold flow mass distribution calculated at the ZOM 
plane should account for the influences of combustion' gas acceleration and RSS, 
as previously suggested. Another significant omission in DER mass distribution 
modeling is the characterization of turbulent mixing effects on performance 
and chamber compatibility (remembering that DER is the computational base for 
the Injector/Chamber Compatibility (ICC) model). Turbulent mixing effects down- 
stream of ZOM could range from minimal (for uniform patterns with a large number 
of injector elements) to substantial (coarse patterns or film cooled chambers). 
The DER streamtube modeling does not calculate any inter-streamtube mass ex- 
change downstream of the ZOM plane. The characterization of turbulent mixing 
effects would provide a large step in the direction toward providing DER, with 
the desired a priori prediction capability. The effect of streamtube mixture 
ratio changes due to turbulent mixing on propellant vaporization efficiency can 
be shown to be a second order influence for most chamber designs. Therefore, 
the turbulent mixing and vaporization calculations could be separated, resulting 
in a relatively simplified computational approach. Conversely, a simultaneous 
mixing with vaporization model would result in a more accurate solution with an 
inherent increase in programming complexity and computer run time. 

Two test programs have been conducted that resulted 
in the development of semi -empirical models for the prediction .of mixing limited 
injector performance. The programs investigated gas/gas intraelement (Ref. 11) 
and film cool in'g/injector core mixing (Ref. 38).- The physically mechanistic 
analytical modeling developed on both programs is naturally extendible to the 
DER program. 
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O'Hara, et. al., (Ref. 39) have recently completed 
work that resulted in a quantification of the intensity of turbulence and 
Lagrangian correlation for turbulent mixing in rocket combustion chambers. 

The analysis was based on gas sample measurements taken from a oxygen/heptane 
300 psia chamber pressure small rocket engine. It is concluded that the Lagrangian 
correlation could be used in rocket engine diffusion calculations. Another 
significant discovery of the work was that turbulence intensity was very high 
near the injector face due to the rapid rate of combustion and presence of 
liquid spray in this area. 

7. • Development of an A Priori PER Mixing Model 

'• ' Two primary weaknesses in the DER ZOM plane mixing "technique 
were determined during the DER review. (1) No methodology exists for calculating 
the correct ZOM plane value on an a priori basis, and (2) RSS and turbulent 
mixing effects can significantly alter the correct mass distribution for chamber 
throat plane mixing loss calculations. 

Recent research indicates that RSS is a near-zone phenomenon. 
That is, combustion gas mass distribution is affected by RSS within an inch or 
two of the injector face plane. Combustion gas acceleration effects on spray 
distribution are also predominant in the near zone. (See Section. V of this 
report). Conversely, turbulent mixing effects can extend to the chamber throat 
plane. For these reasons modeling of gas acceleration effects, RSS, and turbu- 
lent mixing can be easily adapted to the current DER computational methodology, as 
graphically suggested in Figure A-19. 

New DER mixing methodology could be executed in the 
following three steps. 

(1) ZOM is calculated with the gas acceleration effects 
model presented in Section V of this report. The 
ZOM mass distribution would be caTculated based on 
empirical cold flow spray coefficients (identically 
to the current model ) . 
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Figure A-l 9. Proposed DER Mixing Model Approach 







(2) The ZOM cold flow distribution would be adjusted 
for RSS effects. 

(3) The turbulent mixing model would adjust the new 
ZOM mass distribution in axial increments to the 
chamber throat plane. 

D. Inconsistencies Between JANNAF Procedures and DER Computer 

Program Operations 

The JANNAF Performance Prediction Procedures described in CPIA 
246 were reviewed' with regard to use of the DER computer program. , The primary 
functional purpose of DER in the rigorous JANNAF procedure is to provide STC 
output which can be directly input to the Two Dimensional Kinetic (TDK) Reference 
Program (Ref. 5). TDK analyzes the supersonic nozzle expansion process in the 
JANNAF methodology. During Task II of this program the STC/TDK interface problem 
was to be objectively evaluated. Rescoping of Task II eliminated STC and TDK 
analyses that would have determined if the interface is easily accomplished. 

Review of CPIA 246 indicated several inconsistencies between 
the JANNAF methodology and the results of the DER review. Each point is 
elaborated on to the extent that task scope allowed. It is hoped that the 
significant questions will be resolved with future DER review and applications 
work. 


1. The DER Review results indicate that the intended predictive 
accuracy of the JANNAF rigorous procedure to within 1 percent for predicted 
specific impulse) is out of the question for a priori performance prediction. 

ZOM plane, mass median drop size, drop size distribution, droplet drag coeffi- 
cient, and combustion effects on mixing considerations can each affect the 
prediction on the order of 1 percent. 

2. Section 2.8.1 of CPIA 246 suggests that the DER subcritical 
K-Prime vaporization model is valid for chamber pressures 20 percent below 
propellant critical pressures in which the droplet heating time to the "wet bulb 
temperature is negligible and combustion gas solubility is low. For other cases 
the DER supercritical version is recommended for use. The problem is that the 
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JANNAF methodology includes no formulation for predicting the significance 
of the droplet unsteady temperature state or gas solubility effects. 


3. Section 2.8.3 of CPIA 246 recommends that the most physically 
realistic technique for selection of ZOM is to run LISP to a ZOM plane value 
at which the spray patterns from different elements start to overlap. It is 
almost impossible to recognize this point in a typical LISP output unless one 
possesses an intimate knowledge of the injector element spray characteristics. 
Also, this technique ignores interelement mixing effects on injector perfor- 
mance. 
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This appendix details the results of the critical review of the 
JANNAF CICM computer model. Three subtasks were accomplished during the 
review. (1) Identification of Operational Problems Including a Code Review 
and Inclusion of • Diagnostic Print-Out Statements, (2) Identification of 
Inadequate Formulations and Model Technical Shortcomings, (3) Review of the 
JANNAF Performance Prediction Procedure (CPIA 246) with Regard to the Use 
of CIGM and Identification of Inconsistencies Between the Procedure and 
Program Operations . A complete summary of recommendations resulting from 
the CICM review is included in Section III.B. of this report. 

A. Coding Errors and Diagnostic Statements 

The review of the CICM computer code consisted of verification 
of the key model equations'. The equations checked in the code were: 

(1) jet stripping rate, (2) mass median drop size, (3) droplet drag coeffi- 
cient, (4) droplet drag force and acceleration, (5) droplet heating, and 
(6) droplet vaporization. These formulations were all coded correctly. 

One possible code (or formulation) error was discovered during 
the review. Examination of the documented sample case output revealed that 
periodically, in the chamber vaporization calculation, drop size groups that 
should' not have been completely "vaporized are dropped from the calculation. 

As an example, refer to pages 117-122 of the Appendix C sample case output 
in the user's manual ( Ref. 4)-. At the 0.75 inch axial station there are 
22 drop size groups.. At the 1.00 inch axial station the number 2 drop size 
group is missing, though the smaller drops in group number 1 still remain. 

At the 0.75 inch station the group 1 D was 75 microns and the group 2 IT was 
88 microns. It is physically incorrect that the group 2 drops would vaporize 
more quickly than group 1 drops. This error occurs again at the 1.250 axial 
station when the number 3 drop size group vanishes, but the group I still 
exist. The source of this computation error was not discovered during the review. 
It should be added that, for the following reason, this possible error was 
not expected to affect the CICM analysis of the M-l engine reported in report 
Section IV. The drop size groups in question were consistently vaporized com- 
pletely long before the chamber throat plane calculation station was reached. 
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B. Identification of Inadequate Formulations arid Model Technical •• 

Shortcomings 

The purpose of this subtask of the CICM review was to identify 
model technical problem areas and, if possible, to propose approaches for 
improvement. Incorporation of the improvements was, for the most part, 
beyond the scope of the program, but their identification provides a basis for 
future CICM work. The CICM user's manual states that the model controlling 
processes are": (l)the local stripping rate of the liquid jet-, M ft , (2) the 
local mean drop' size produced when M^ is stripped from the jet, IT, (3) 'the 
droplet heating and vaporization rates, (4) the assumed droplet drag coefficient 
formulation, and (5) for the chamber flow, the rate of mixing of the external 
“rigimesh" face flow. Careful review of the program input requirements and 
model analytical assumptions indicated that two additional important controlling 
CICM parameters should be defined, (6) the input specification of the intra- 
element fuel and oxidizer mass and mixture ratio distribution, and (7) the 
input specification of separate flow analysis zones to allow for manifold mass 
and mixture ratio maldistribution. The-CICM program does not calculate 'mixing 
and requires that mass distribution input be user justified. Currently no 
standard guidelines exist in the CICM user's manual or the JANNAF Performance 
Prediction Manual. (CPIA 246, Ref. 3) for measurement or input specification of 
.these propellant mass distributions. The results of the CICM formulations review 
are presented below in seven sub-sections that deal with the model controlling 
processes defined. above. 


1 . . Jet Stripping Rate Correlation 


defined below. 


The circular stripping rate correlation used by CICM is 
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The CICM and JANNAF open literature does not include a derivation of the 
stripping rate correlation. A cursory examination of the literature on- the 
atomization of liquid jets injected concurrently into gas streams yielded no 
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directly applicable correlation for jet disintegration rate. The scope of 
the CICM review did not allow for a comprehensive literature review on liquid 
jet atomization. Qualitatively, the equation appears to be correctly formu- 
lated. As the jet to gas relative velocity increases the stripping rate in- 
creases due to aerodynamic drag. The stripping rate will also increase as the jet 
density, viscosity, and surface area (irD. Az) and the gas density increase. 

u 

It is also correct that the stripping rate should decrease as the jet liquid 
surface tension increases. For the atomization constant (C ft ) to be a universal 
constant the respective- terms in the atomization equations must be raised to 
the correct power. Also, no physical variables that have a significant in- 
fluence on the stripping rate can be omitted from the equation and still result 
in development of a universally applicable value of C^.' Variables that could 
possibly fall in this category are the absolute liquid jet velocity, the absolute 
gas stream velocity, and the gas stream viscosity!. 

The stripping rate equation calculates the time lag between 
jet initial contact with a concurrent gas stream and final jet disintegration. 

For a coaxial injector, the initial contact can occur in the recessed portion 
of the element cup or at the injector face plane. Typical gas/liquid coaxial 
injector designs require relatively long chamber lengths to reduce mixing and 
vaporization performance losses. Therefore, the atomization time lag is usually 
small compared to the total chamber residence time. As an example, the M-l 
engine design analyzed during Task II has a conical chamber length (face plane 
to throat plane) of nearly 30 inches. Based on documented previous CICM runs 
it was expected that the element oxidizer jet would be completely atomized from 
2 to 4 inches of the injector face. It was apparent that the drop sizes 
calculated to be shed from the oxidizer jet will have a far more significant 
effect on M-l predicted engine performance than the rate of jet atomization. 

2. Mass Median Drop Size Correlation 


defined below. 


The mass median drop size correlation used by CICM is 
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Similarly to the jet stripping rate correlation described previously, the 
open literature does not include explanation of development of the CICM drop 
size correlation equution. However, a number of investigations are documented 
that concentrated on measurement of drop sizes generated from injection of 
liquid jets into non-accelerating concurrent gas streams. Table B-I summarizes 
the results of three of these investigations and compares their correlation 
results to the CICM equation. The table shows the power exponent correlated 
by each study for nine different independent variables.- The Ingebo study 
(Ref. 13) developed the most mathematically stringent correlation by assuming 
that the measured maximum drop size was controlled by six non-dimensional para- 
meters that characterized liquid hydraulics, gas dynamics, gas acceleration, 
and liquid hydrostatic and surface tension forces. The CICM correlation accounts 
for five of the nine variables modeled by Ingebo. The most significant CICM 
omission identified by Ingebo appears to be that the absolute gas velocity 
(not just the velocity differential) has a significant influence on the atomized 
drop size. Comparison of the exponents of the variables modeled by both Ingebo 
and CICM shows that the exponent sign .agrees for all variables except the liquid 
jet viscosity. The Ingebo result is inconsistent with other atomization models 
that indicate drop size increases with increasing liquid viscosity. For example, 
Priem (Ref. 17) uses a propellant properties grouping that results in a liquid 
jet viscosity exponent of + 0.25. The other variable exponents shown agree 
in sign but consistently disagree on the absolute magnitude. It is apparent 
that the cited jet drop size equations differ because of the influence of 
measurement technique, measurement error, and the method of data correlation. 

The CICM code was examined to determine the drop size dis- 
tribution relationship used by CICM as an addition to the drop size equation 
review. During the DER review (Appendix A) it was shown that the assumed drop 
size distribution, about the mass median diameter, significantly influences the 
predicted total mass vaporization rate. Review of the CICM code indicated 
that subroutine ATOM calculates the portion of the liquid jet that is atomized 
over one axial computational increment. ATOM calculates a droplet spray group 
based on the total jet mass shed during the axial step. All the drops in the 
group are assigned an initial' diameter equal to the mass median diameter cal- 
culated with the previously introduced drop size correlation equation. Thus, 
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TABLE B-I. COMPARISON OF LIQUID JET 
BREAKUP CORRELATIONS IN NONACCELERATING GAS STREAMS 
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in reality, CICM does not calculate a real droplet size distribution. CICM 
assumes that all the mass shed during a finite time period, defined by the 
axial step distance, is shed with a constant diameter defined by the F equa- 
tion. The influence of this calculational assumption on the total liquid 
vaporization rate can not be estimated simply. It is apparent,' though, that 
the drop size distribution tabulated at the end of a CICM run is only the summa- 
tion of several constant mass median diameter groups, each group being cal- 
culated over a particular axial step. This resultant distribution is quite 
different than a drop size group calculated with distributions typicallyused 
to model rocket' combustor sprays. 

3. Droplet Heating and Vaporization Rate Formulations 

CICM contains an advanced droplet heat-up and vaporization 
model that is described in detail in the program user's manual; The CICM 
formulation is far superior to the DER K-Prime model in that the droplet 
temperature transient and continuous vaporization through subcritical and 
supercritical propellant states are allowed for. The CICM/STC interface 
procedure review reported in Section C.2 of this appendix resulted in a recom- 
mendation that, in the JANNAF performance prediction methodology, CICM should 
compute to the chamber plane at which all the calculated drop size groups have 
reached the chamber "wet bulb" temperature. For oxygen this temperature trans- 
ient typically takes place over a time period equal to about 10 percent of the 
-total time required to vaporize 99 percent of the propellant (Ref. 17). After 
CICM has calculated the transient the STC program calculates droplet steady 
state burning from the interface plane to the chamber throat plane. Therefore, 
STC is responsible for calculating droplet vaporization rates over approximately 
90 percent of the total droplet chamber residence time. 

The two most important functions of CICM in the JANNAF 
methodology, based on the information in the previous paragraph, are to; (1) 
calculate liquid drop sizes resulting from aerodynamic stripping of the liquid 
jet; the bulk of which will be vaporized in STC, and (2) calculate the droplet 
temperature time transients. The CICM drop size correlation was discussed in 
the previous sub-section. The CICM droplet heat-up formulations were checked 
by comparing, for reference, CICM calculated heating rates for oxygen to 
heat-up rates calculated in Ref. 17. The results of the comparison are shown 
in Figure B-l . Initially, heat-up rates were compared, as a function of 
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Mass Median Drop Diameter, D (microns) 


FIGURE B-l. • 


COMPARISON OF OXYGEN HEATING RATE CALCULATIONS 



mass median drop diameter, for the solid line conditions shown in the figure. 
The agreement- between the two models improved when chamber pressure and velo- 
city differential effects on the local drop, heat transfer rate were accounted 
for in the C I CM calculation. Since no attempt was made to ensure that both 
model predictions were made with exactly the same gas and liquid properties 
the agreement can be considered to be excellent. The, results of the droplet 
heating rate comparison verify the CICM droplet heat-up model. 

4. Droplet Drag Coefficient Correlation • 


Cl CM employs the droplet' drag correlation equations developed 
By Rabin (Ref. 31). ‘ 


. Re D < 80 (B-3) 

80 £ Re D < 10 4 (B-4) 

Re D > 10 4 (B-5) 

me influence of the assumed droplet drag correlation on the vaporization 
rate was examined during the DER review. This investigation indicated that 
the droplet drag correlation significantly affects the final performance 
prediction made by STC. The review recommended that the Rabin drag coeffi- 
cient correlation be reviewed and compared to other available correlations. 

5. Chamber Mixing of Face ."Rigimesh" Flow 

The CICM' program calculates mixing of any face rigimesh - 
flow by assuming the rate of mixing to be a linear function between the face 
and an input downstream distance. The mixed flow is spread uniformly over 
the cross-sectional area of the element flow field and becomes part of the pro- 
pellant (usually fuel) to be reacted. The CICM user's manual states that 
calculations have been performed for rigimesh mixing that indicate that rapid 
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acceleration reduces the rigitnesh flow area to only approximately 3 percent of 
its inj.ection area. This rarefaction occurs on the order of only 2 inches 
from the injector face plane. The rigimesh area reduces to an annulus trapped 
between coaxial element flows; the average thickness of the annulus was cal- 
culated to typically be on the order of .01 inch. It is therefore argued 
that turbulence sweeps the flow into adjacent element flow fields. This cal- 
culational technique would seem to be satisfactory for ordinary amounts of 
rigimesh flow (on the order of five percent of the total fuel flow) because 
the axial variation of the expansion area of the combusting flow field of 
adjacent elements will not be affected significantly. As an example of a 
typical case in point, the M-l injector design analyzed during Task II of the 
program has 3 - percent of the total hydrogen fuel flowing through the riginiesh 
portion of the injector face. 

6. Intra-Element Mass Distribution Specification 

CICM allows for the effect of intra-el ement mass and 
■mixture ratio distribution through user input specification. For each zone 
(i.e., single element) analyzed by CICM, the user is instructed to input radial 
zonal oxidizer and fuel mass fractions based on single element cold flow data. 
An example of such input is shown in Figure B-2 , taken from the CICM J-2S 
sample case in CPIA 246. There are several problems associated with accounting 
for intra-element mass non-uniformities in this manner. 

(1) There is no available standard technique for 
measuring single element cold flow gas/liquid coaxial mass distribution. 

(2) The JANNAF methodology does not specify the axial 
plane (i.e., collection plane) at which the intra-element mass distribution 
should be specified. Face plane measurements are most easily accomplished 
but will be significantly altered by the high AV shear mixing inherent to 
coaxial element designs. 

(3) The test cases used to back out the recommended 
atomization and drop size input constants to CICM assumed that the thrust 
chamber in question had uniform throat plane mixtue ratio distributions. For 
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most real coaxial injectors there will be a finite mixing loss because the 
coaxial element is a relatively slow mixing element. It is apparent that 
the correct values for the and coefficients will be directly dependent 
on the assumed single element mixture ratio, distribution. Unless a standard 
method for measuring or calculating single element mixture ratio distributions 
is developed it is extremely doubtful that universal values for the and 
constants can be verified. 

(4) Similarly to the DER program for 1 iquid/1 iqu'id in- 
jectors, CICM does not allow for the influence of combustion on the single 
element mass and mixture ratio distribution. 

Currently, it appears that, without a standard coaxial 
element mixing model or approach, standardization of the parameters that 
influence the propellant vaporization rate will be difficult. That is, two' 
processes affect coaxial injector performance (mixing and vaporization) and each 
process must be physically modeled to a comparable degree to result in a model 
that can calculate an accurate superimposed solution.' At this stage CICM has 
been verified for engines that apparently have only one effective performance 
loss mechanism, i.e., incomplete propellant mass vaporization. 

7- Manifold Mass Distribution Zone Specification 

The CICM user's manual recommends that measured or calculated 
manifold mass maldistributions should be accounted for by modeling separate 
chamber flow field' zones when executing the program. Since the manifold distri- 
bution is usually calculated at the injector face plane, this technique assumes 
that the mass maldistribution will persist to the chamber throat plane. The 
manifold mass maldistribution performance loss is inevitably overpredicted with 
this technique, since turbulent mixing is ignored. Also, there is no recommended 
methodology for dividing the measured distribution into analysis zones. There- 
fore, the method of zone mass fraction assignment also becomes an user controlled 
input that affects the final performance prediction. 
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The solution of this problem is more complicated than the 
intra-element mixing problem discussed in the previous sub-section. It would 
be difficult to generalize a chamber zonal mixing model . Measurement of perfor- 
mance for thrust chamber assemblies having- negl ig.ible vaporization and intra- 
element mixing losses would seem to provide a reasonable approach for solution 
of this problem. That is, if the engine vaporization and single element mixing 
losses are small (or can be accurately calculated)- the manifold induced mal- 
distribution loss 'can be backed out from the performance data. 

C. Inconsistencies Between JANNAF Procedures and Program Operations 
1. - Background 

CICM was developed as a rigorous analytical model that 
describes the atomization, vaporization and combustion of gas/liquid coaxial 
jets in a rocket engine environment. In the context of the JANNAF series of 
performance prediction models, CICM is intended to replace the LI SP sub program 
of DER for gas/liquid coaxial elements. 

CICM-is a highly specialized program intended to be used 
for one specific injector design concept. Additionally, the model has only 
been applied, to this date, to coaxial injectors using a central liquid 0 2 
circular core surrounded by a gaseous H 2 or H 2 /0 2 combustion gas mixture 
annulus. The program input requires an extensive group of- propellant property, 
cards (644 cards for the CICM user's manual sample case) that will have to be 
generated for each propellant combination analyzed in the future. Another 
factor that currently limits the generality of CICM in the JANNAF methodology 
is that key empirical atomization rate and drop size constants, that control 
program performance predictions, have only been determined from test data 
for the LOg/GH^ propellant combination. Detailed discussion on these program 
inputs is included in Section B of this appendix. The CICM analysis documented 
in Section IV of this report was also restricted to a L0 2 /GH 2 injector. There- 
fore, there are no current plans for testing the ability of CICM to model gas/ 
liquid coaxial designs using other propellant combinations. 
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During this phase of the review task CICM was critiqued 
for its ability to function as documented in the JANNAF rigorous performance 
prediction procedure described in CPIA 246. The evaluation emphasized two 
areas; (1) test and evaluation of the CICM/STC interface procedure by running 
the CICM sample case and subsequently using the CICM input to generate an 
input deck for STC, and (2) development of a criteria 'for specifying the 
chamber axial location of the CICM/STC interface plane. 

2. Evaluation of the CICM/STC Interface Procedure 


The CICM/STC interface procedure was examined carefully 
to ensure that the JANNAF performance prediction methodology accuracy and 
utilization time is not being compromised by the currently recommended 
interface technique. It was determined that the CICM interface routine DERINI 
was incomplete and punched improperly formated cards for input to the STC 
subcritical K-Prime version. The CICM user's manual states that .DERINI punches 
input for the supercritical version of the DER program. No check was made 
to see if the punched output was compatible with the input requirements of 
that DER version. The next section of this appendix completely details the 
interface evaluation and development of a new interface procedure. This new 
DERINI version was successfully utilized during the program Task II CICM 
analysis effort.. 

3. Description of Improved CICM/STC Interface Procedure 

The JANNAF Performance Prediction Manual (CPIA 246} specifies 
that CICM will replace the LISP model for the analysis of gas/liquid coaxial 
elements. In this function, CICM must be capable of calculating spray forma- 
tion, vaporization, and gaseous combustion and of generating output which is 
consistent with STC input and operational requirements. 

The CICM/STC interface procedure was critically evaluated 
during the CICM review task of the Injection Processes Program. The CICM 
sample case documented in the user's manual was executed to determine if an 
interface with the STC program could be easily accomplished. The documented 
CICM sample case 'considers two injector zones (elements) which are each 
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divided into two intra-element mixture ratio zones through input mass fraction 
distributions. This input specification results in four separate sets (2 
interelement zones x 2 intraelement zones) of streamtube input for the STC 
subprogram of DER. The DER user's manual was then used to determine that this 
four streamtube case required an input deck consisting of eight-six separate 
cards. The CICM interface. subroutine (DERINI ) was designed to punch only the 
streamtube flowrate and drop size input cards required by STC (cards 6720, 
7010-7016, 7020-7026, 7030-7036, and 7040-7046 for this case). These cards 
comprised twenty-nine of the eighty-six cards required to correctly interface 
the CICM output and the DER input. Also, the cards punched to designate the 
streamtube and drop size group droplet flowrates (GWSPR (I, J)), velocities 
( VELD! (I, •«])), and diameters (GDIAD1 (I, J)) were improperly formated to 
be input to STC. The format error occurred because the interface subroutine 
DERINI also punched droplet temperatures for each streamtube and drop size 
group, while STC (in the subcritical DER K‘ version) does not require or allow 
for this input. 

This attempt to join the CICM sample case output with the 
subcritical STC program indicated that the interface procedure required improve 
ment. The six improvements that were grouped to result in the new interface 
procedure are detailed below. 

(1) The streamtube and droplet size group input cards 
were properly formated and labeled with their correct sequence numbers, 

as indicated in the, DER user's manual. 

(2) A number of STC inputs that have constant values 
when LISP execution does not precede STC- execution (e.g., when CICM inter- 
faces with STC) were assigned- values in the CICM interface routine DERINI 
and included in the STC input cards to be punched. 

(3) All STC inputs that are also input or internally 
calculated in CICM (e.g., combustion gas transport properties as a function 
of mixture ratio) were included in the DERINI output deck for input to STC. 

(4) STC inputs that could not be specified as constants 
and are not set by CICM input or calculation were grouped .into a new namelist 
input set for DERINI. 
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(5) Coding was included in DERINI to result in each 
STC input card having its correct sequence number, as specified in the DER 
user's manual, punched in columns 73-80. 

(6) An option was included in the new DERINI namelist 
input group to allow for writing the DERINI formulated STC input deck on a 
computer system drum file (or stratch tape) without having to punch an actual 
card deck. 


The listed improvements resulted in an interface technique, 
completely -internal to CICM, that allows generation of all required STC input. 
The new CICM/STC interface procedure is detailed *in the following three 
paragraphs that include in turn; (1) a listing of the new CICM interface routine 
DERIN and specification of required line changes and additions to generate the 
new routine from the old verison, (2)' a description of the required namelist 
input for DERINI, and (3) a description of the STC input deck generated with 
the new procedure for the CICM sample case documented in the program user's 
manual . 


A compilation of the new version of the CICM subroutine 
DERINI, designed to provide punched card or mass storage file input to the 
DER subprogram STC, is shown in Table B-II. The line modifications that were 
applied to the original' version of DERINI are detailed in Table V-III. No 
other changes are required to any CICM routine to develop the new interface 
procedure.' 


The required namelist input variables for DERINI are 
defined in Table B-IV. The DERINI namelist variable inputs must be preceded 
by a £STC specification and followed by a gEND specification (or the system 
equivalent to these Univac 1108 Exec-8 designations). The first three variables 
listed in the table designate forms by which the STC input data may be output 
from DERINI. Any combination of these three output forms may be specified. 

The remaining input variables listed are identical to descriptions given in the 
DER user's manual. The oxidizer latent heat input, DHVO, should be consistent 
with the droplet "wet bulb" temperature calculated by CICM at the interface 
plane axial chamber location. Any STC variables not listed are either input 
to, internally set within, or calculated by CICM. Liquid fuel properties- are not 
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included in the namelist because CICM requires that one propellant be gaseous 
and one liquid. The namelist input set used to check the new interface technique 
for the CICM sample case is shown in Table B-V. 

The STC input generated by DERINI for the CICM sample case 
is listed in Table B-VI. The STC program was successfully executed with the 
data set shown. 
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TABLE B-II MODIFIED CICM/STC INTERFACE SUBROUTINE (DERINI) 
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3F0R»9 OERIN1 

FORTRAN y», ISO VERSION «, 68. 02/26/76.06 l52t 12 (6») 


SUBROUTINE DERINI ENTRY POINT 002154 


STORAGE USEDI COOE(l) 002370) 0ATA{01 006755} BLANK COMKQN<2) 000000 
COMMON BLOCKS 1 

0003 TCPLF 0015J4 

0004 PROP 1* 000013 

0005 CCTA8C 000135 
000* CHCOM 000072 


EXTERNAL REFERENCES (BLOCK# NAME) 

0007 U!Cf*C 

0010 XVDHV 

0011 EQSTAT 

0012 CGPRQP 

0013 RHOGF 

0014 NRDUS 

0015 NIOtS 

0016 NI02S 

0017 SORT 

0020 NRNL* 

0021 XPRR 

0022 NkOU$ 

0023 NERR3S 


STORAGE ASSIGNMENT (BLOCK# TYPE# RELATIVE LOCATION, NaME) 


oooo 


006503 

IP 

0O01 


000003 

10L 

0001 


002224 

10576 

OOOl 


002244 

1071G 

OOOO 


006545 

l IF 

0001 


002273 

nooe 

0001 


001330 

121 

0001 


000270 

1P0L 

OOOl 


001334 

13L 

OOOl 


001340 

1«L 

0001 


000307 

1401 

• 0001 


000027 

1 45G 

0001 


001343 

15L 

OOOl 


000042 

1546 

OOOl 


•001316 

161 

0001 


000312 

160L 

0001 


000066 

1666 

OOOO 


006550 

17F 

oooo 


006551 

16F 

oooo 


006506 

2F 

OOOJ' 


0004*6 

200L. 

0001 


000523 

210L 

OOOl 


000160 

217G 

0001 


000170 

2246 

OOOl 


000227 

24JG 

0001 


000246 

247G 

0001 


000272 

261G 

oooi 


000327 

2766 

oooo 


006511 

3F 

OOOl 


000647 

300L 

0001 


000346 

307G 

0001 


000537 

340G 

OOOl 


000653 

S55G 

OOOl 


000674 

3626 

oooo 


0065 1 4 

4F 

0001 


001054 

4146 

0001 


001100 

426G 

0001. 


001273 

4726 

oooo 


006516 

5F 

OOOl 


001404 

5566 

0001 


001441 

5630 

0001 


0014*1 

576C 

oooo 


.006520 

6F 

,0001 


001067 

610L 

oooi 


001512 

620G 

0001 


001073 

62QL 

0001 


001257 

640L 

OOOl 


001574 

6606 

OOOl 


001624 

6706 

OOOl 


001650 

6775 

oooo 


006522 

7F 

0001 


001267 

700L 

OOOl 


001712 

715G 

oooi 


001721 

7 20G 

OOOl 


001751 

727G 

oaot 


002004 

7416 

0001 


002013 

744G 

OOOl 


002043 

7536 

oooo 


006532 

8f 

oooo 


006534 

9F 

0000 


0066Si 

900QF 

0001 


002310 

99991. 

0006 

R 

000046 

ACHAMC 

0006 


000014 

ACSC 

oooo 

ft 

004110 

AGO 

0005 


000002 

AMRT 

oooo 

R 

004350 

AREA! 

oooo 

R 

006424 

ARK . 

oooo 

R 

006353 

ARTOLD 

oooo 

ft 

006425 

8RK 

0006 


000012 

BSPRC 

oooo 

R 

006442 

C 

0006 


000017 

CCANGC 

„ 0006 


000015 

CLNTC 

0006 


000016 

CONRAC 

0000 

R 

006255 

CPVO 

oooo 

R 

006352 

CRTOL 

0006 


000013 

CSPRC 

•oooo 

R 

006203 

C 51 R 

0004 


000011 

CXOV 

0006 


OOOO 1 1 

OELTXC 

oooo 

ft 

006423 

DHV 

oooo 

R 

006470 

DHVF - 

oooo 

R 

006357 

DHVQ 

oooo 

ft 

005107 

DOO • 

0000 

R 

006430 

DRLDT 

oooo 

ft 

006431 

DNS 

oooo 

R 

006372 

EMRCG 

0006 


000004 

EMRGJC 

oooo 

R 

006432 

EMRI 

0000 

R 

004160 

EM** 

oooo 

R 

006377 

EMhCS 

0006 


000006 

emngjc 

000.4 

R 

000004 

EMWL 

‘ 0004 


000012 

EHfcpR 

0004 R 

000006 

EMWV 

OOOQ 

ft 

006443 

FA 

oooo 

R 

006400 

FCHA 

oooo 

R 

006362 

FCHAM 

oooo 

R 

004610 

FFMJK 

oooo 

R 

0064 1 7 

FN 

oooo 

R 

004660 

FUMIX 

oooo 

R 

004730 

FSOER 

oooo 

R 

006406 

FI 

oooo 

R 

003770 

GAM 
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0008 


000007 GAHGJt 

0000 

R 

oooooo 

GA8FU 

OOOO 

R 

001060 

GO I ADI 

OOOO 

R 

005607 

•GMR 

OOOO 

R 

006413 

G8H8 

0000 

« 

005763 

GT* 

0000 

R 

002020 

GTOOt 

oooo 

R 

002760 

GVELD1 

9000 

R 

000120 

GK0PR 

oooo 

R 

006421 

HD 

0000 

I 

006365 

t 

0000 

I 

006440 

IC 

oooo 

T 

006472 

ICARO 

oooo 

I 

006347 

ICRC 

oooo 

I 

006355 

IDRUM 

0000 

I 

006450 

IFJLE 

0000 

I 

006410 

II 

oooo 

I 

006420 

U 

oooo 

I 

006451 

IlISP 

ooco 


006715 

INJPS 

0000 

I 

006351 

Iprmst 

0000 

1 

006350 

1PRSST 

oooo 

I 

006356 

I PUNCH 

oooo 

I 

006467 

IPUN3D 

oooo 

1 

006455 

1ST 

0000 

I 

006452 

I STC 

0000 

I 

006454 

I TDK 

oooo 

I 

006453 

ITRANS 

oooo 

I 

00b35« 

IWRITE 

oooo 

1 

006405 

11 

0000 

I 

00o360 

J 

oooo 

I 

00640? 

JJ 

oooo 

I 

00o4)2 

JJJ 

oooo 

I 

006335 

JK 1 

OOOO 

I 

006336 

JSPC 

0000 

1 

00b4t 1 

K 

0000 

I 

006500 

L 

oooo 

I 

006475 

M 

0006 


000001 

M2C 

oooo 

1 

006476 

N 

0000 

I 

006447 

NASE6 

oooo 

I 

0 064 7 4 

NC 

' 0006 

I 

OOOOOO 

NCHAMC 

0006 


000002 

NL0N4C 

oooo 

1 

006361 

NDER 

0000 

I 

006366 

NOSC 

0000 

I 

006367 

Nf LEM 

oooo 

I 

006446 

NGF 

oooo 

I 

006364 

NGO 

0 00 Q 

I 

006444 

NGT 

0000 

I 

006363 

NMJXZ 

oooo 

I 

006346 

NMSTI 

oooo 

I 

006464 

NOZUN 

0009 

I 

006344 

NP 

0003 

I 

001512 

NPCP 

0000 

I 

006345 

NSSTI 

oooo 

I 

006445 

NST 

oooo 

I 

006465 

NS7PZ 

0005 

I 

OOOOOO 

NT AB 

0003 

I 

0 0 1 5 1 Q 

NTCP 

0000 

I 

006337 

NIK 

oooo 

I ' 006466 

NUG 

oooo 

R 

003720 

P 

oooo 

R 

006370 

PC 

oooo 

fi 

00644 l 

PCI 

0004 


000000 

PCR1T 

oooo 

■R 

006«36 

PHI GH 

oooo 

R 

006437 

PLOW 

oooo 

R 

00«230 

PUS 

0006 


000020 

RCBCC 

0006 


000021 

RC7C 

OOOO 

R 

006427 

RG 

oooo 

R 

006426 

RHOO 

oooo 

R 

004540 

RHOG 

0013 

R 

oooooo 

HHQGF 

0000 

R 

00646 1 

RhDLF 

oooo 

R 

00O343 

RHOtO 

oooo 

R 

006460 

rmonbf 

oooo 

R 

006342 

RHONBO 

oooo 

R 

000050 

SMR6 

0000 


006553 

STC 

0006 


000005 

STGJC 

0004 


000010 

STUCMR 

0005 


OOOOOJ 

5TT 

QOOQ 

R 

006414 

SUM1 

0000 

R 

0064J5 

SUM2 

oooo 

R 

006416 

SUMS 

oooo 

R 

006373 

SWSPR 

oooo 

R 

006456 

TBF 

oooo 

fi 

006341 

TOO 

0000 

R 

006305 

T CON VO 

0003 


000050 

TCP 

0004 

« 

000002 

TCRIT 

0005 

R 

000047 

TGAM 

oooo 

H 

004470 

TGAS 

0000 

R 

005563 

TMl 

0003 

R 

000670 

THOL 

oooo 

R 

00b376 

TLI 

0005 

R 

000003 

TMR 

0005 

R 

000071 

TMW 

0000 

R 

006457 

TNBF 

oooo 

R 

006340 

TNBO 

• oooo 

R 

004040 

TO 

oooo 

R 

005253 

TDD 

0003 

R 

000024 

TPCPL. 

0003 

R 

000000 

TTCPt 

0005 

R 

000025 

TTU 

0005 

R 

000113 

TVIS 

oooo 

R 

006225 

TVO 

oooo 

R 

006374 

VCG 

0000 

R 

00««20 

VGAS 

oooo 

R 

006433 

vise 

oooo 

R 

006375 

VLJI 

oooo 

R 

005417 

VOO 

oooo 

R 

004300 

VUS 

0000 

R 

006371 

WCG 

oooo 

R 

006403 

WFE 

0006 


000003 

WGJC 

oooo 

« 

006402 

WUXE 

oooo 

R 

006404 

WSPE 

0000 

R 

004 743 

KSPR 

oooo 

R 

006435 

NT 

oooo 

R 

006401 

WTE 

oooo 

H 

006462 

WTMLLF 

oooo 

R 

006463 

KTMLVF 

0006 

R 

000022 

XCMAMC 

0006 


000010 

XLMC 

oooo 

R 

006473 

XNAP 

oooo 

R 

006502 

XNGT 

oooo 

R 

006477 

XNMR 

0000 

R 

006501 

XNTK 

oooo 

R 

006434 

XOV 

oooo 

R 

006422 

XV 

oooo 

R 

005471 

ZST ART 






ooioi 

1* 

00103 

2* 

00103 

3* 

00103 

4* 

00103 

5* 

00104 

6* 

00105 

7* 

00106 

6* 

00107 , 

9* 

00107 

10* 

00110 

1 1* 

00110 

12* 

00111 

13* 

00111 

14* 

00112 

15* 

00112 

16* 

00112 

17* 

00113 

IB* 

00113 

19* 

00116 

20* 

00117 

21* 

00120 

22* 

00121 

23* 

00122 

24* 

00123 

25* 

00120 

26* 


SUBROUTINE DErTNI (ID£R,ACHAM,XHInD£, IRDER) 

DIMENSION «ASFL C40 SMRGC40), GWSPRU2»«0>, GDIA0U!2»«O)* 

1 GTUOl Cl2,4Q), GVELDl(t2,40) r P{40), GAHC^O J r TOC40), 

2 iGOfilO), EMhcilo), PUS(40) t VUS(flO), AREAl(40>, VGA$<40>, 

3 TGA$ (40 ) f RMOGC40) 

DIMENSION FFM[X(qO), FOMIX(ftQ), FSOfR(li) 

DIMENSION hSPR(lOO), DUO(iOO), TODUOO), VO0UOO), THttZO) 
DIMENSIUN GmR ( 6, I8),GTK(e,24) # C5TRUB)»TVO(24),CPVO(24)#TCONVOCe«) 
COMMON /TCPLF / TTCPU20,!), TPCPL(20,1), TCP{20,20»1>, 

1 THfU(20,20,l), NTCPC2), NPCP(2) 

CUMMON /PROP 1 / PCHITC2), |C«IT(2># EMWL (2 } * EMWVC2), 

1 STOCMR, CXOV» tMhpR 

CONMON/CGTAec/NTAB,$TT, AHRT,TMR<ltt>, 770(18), TGAM (IB)# THNU8)* 

* TvlS(lB) 

CUMMDN/CHCOM/NCHAHC , M2C,NCON4C, wGJC#fcMRGJC /STGJC,£MHGJC, GAHGJC, 

A XLHC»D£LTXC»8SPRC»CSPHC* ACSC,CLNTC»CDNRACf CCANGC, 

* , RCBCC,flCTC,XCHAMCC2Q),ACHAMC(20) 


QOOOOO10 

00000020 

00000030 

00000040 

00000050 

00000060 

00000070 

00000000 
00 000090 
00000100 
OOOQOitO 


DATA JKt,4SPC/0»l/ 00000120 

1 FORMAT (41l2,24X,I3) 

2 FORMAT (3£l2,6,S6X»I8) 

3 FOR MAT ( 4 E 12 , 6 , 24 X, 10 ) 

4 FORMAT { 1216# 18) 

5 FQRMAT{6£ 12,6,18) 

6 FORMAT (6J 1 2 , IB) 

7 FORMAT (6X , • STC INPUT FROM GJCM PROGRAM CASE > j»7?3»I0) 
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0012 s . 

27* 

< 6 

FORH*T(72X f 18) 



00126 

26* 

9 

FORMAT (6X» 'D2/H2 5*8 PROPERTIES FROM CIC* INPUT * * T7: . I" ) 



00127 

29* 

It 

FORMAT<2Et2,6,4BX,18) 



00130 

30* 

17 

FORMAT ( 1216) 



00131 

31* . 

18 

FORMAT C6E12, 6) 



00132 

32* 


NAMELIST /STC/TVO* TCGNVO*NfK» TN8Q* T30#RM0NBDf RHQLOfNP*NS$7'I -NM37I < 



00132 

33* 

* 

ICRC.IPRSST.IPRMST.CRTOL.ASTOlO.lWRlTE.IDRUM.IPUNCH. 



00112 

34* 

* 

CPVO»OHVO 



00133 

35* 


CSTAR(xm«, TO»GAM)*5QR T {49677, *GAM*TO/XM«/((2,/CGAM+ !#))**( (GAH*1,: 



00133 

36* 

*/(GAM»t,))))/GAM 



00134 

37* 


J * 0 

ooooci®o 


00 5 35 

36* 


NOER s 0 

00000200 


00136 

39* 


FChiAM a 0,0 

00000220 

V a v V V » 

00137 

40* 

10 

RLA0CIH0ER»17) nMJXZ.NGQ 


000003 

00141 

41 * 


HEA0C1RDER, 18) {FFMlX(I),FOMIX(n,I»l,NMIXZ) 

' 

000014 ' 

00152 

42* 


RE A 0 ( IHDER » 18 ) (FSDER { 1 ) , X* 1 .NGO ) 


000033 

00160 

43* 


READ t ISO£R# 17) NDSC.NEIEM 


00004S 

00164 

44* 


REA0(I«DEk,l8) (WSPH u ) ,000(1) ,TOOU),VaD(I),J«l, 100) 


000054 

00175 

45* 


RE AD ORDER# 18) PC , WCS. EMRCG # SWSPR , VCG, VI JI , Ti. I # EMWCG.FCHA 

» 

000074 

00210 

46* 


wTE = {SMSPR+WCG)*NELEM 

■ 00000310 

000 1 12 

0021 1 

47* 


WOXE a (SWSPH + WCG*EMRCG/'(1,+EMRCG))*NELEM 

00000320 

000116 

00212 

48* 


WFE a KTE-WOXE 

00000330 

000126 

00213 

49* 


hSPE a SRSPR*NElEM 

00000340 

000130 

00214 

50* 


CALL LOCF AC ( JK1 ,PC#TPCPL*NPCP( 1 ) . 1 1 #F1 ) 

00000350 

000133 

00215 

51* 


NP * NTCP ( 1 ) 

00000360 

000)43 

00216 

52* 


DO 20 I a 1 # NP 

00000370 

000151 

00221 

53* 

20 

THL(I) » THOUtl*I»l)+Fl*<THOL(lt+l»I#U«THOL{ll*I*l)) 

00000380 

000160 

00223 

54* 


DO 100 JJ»i,NMlXZ 

00000390 

000170 

00226 

55* 


I a J+JJ 

00000400 

000174 

00227 

56* 


GASFUC1) a FFMIX(JJ)*WFE>F0MIX(JJ)*(W0XE*M5PE) 

000004)0 

0001 77 

00230 

57* 


SMRGCI) a FOMlX(JJ)*{HOXE-WSPE)/tFFMIXCJJ)*«P£) 

00000420 

000203 

00231 

Sfl * 


GrtSPR C 1 « I ) » 0,0 

00000430 

000205 

00232 

59* 


GDI ADI { 1 . 1 ) « 0,0 

00000440 

000210 

00233 

60* 


GTOD1 ( 1 ; I) « 100,0 

00000450 

000211 

00234 ■ 

61* 

too 

GVttOl { 1 # I ) « 100,0 


000213 

00236 

62* 


IFCNGO.LT.NOSC) GO TO 160 

00000490 

000216 

00240 

63* 


DO 140 IIsl.NGO 

00000500 

000222 

00243 

64*. 


I = II+l 

00000510 

000234 

00244 

65* 


IF U1.GT.NDSC) gO TO 120 

G 0000520 

000237 

00246 

66* 


DO 110 KS1.NMJXZ 

00000530 

000246 

00251 

67* 


JJ » J+K 

00000540 

000246 

00252 

68* 


Gh/SPR ( I . JJ ) * WSPRUI)*NEUM*fQMIX(K) 

00000550 

000251 

002S3 

69* 


GDIAOKI.JJ) ■ POD (ID S 

00000560 

000256 

00254 

70* 


GTnOKI.OJ) a TODUI) 

00000570 

000260 

00255 

71* 

no 

GvELOl ( I . J J ) * vaocin 


000262 

00257 

72* 


GO TO 140 

00000600 

000266 

00260 

73* 

120 

DO 130 K * lfNMJXZ 


000272 

00263 

74* 


JJ a J *K 

00000630 

000272 

00264 

' 75* 


GWSPRU.JJ) « 0,0 

00000640 

000277 

00265 

76* 


GDIAD1U.JJ) * 0,0 

00000650 

000300 

00266 

77* 


GTUOHI.JJ) » 100,0 

00000660 

000301 

00267 

78* 

130 

GVELOIU.JJ) * 100.0 


000303' 

00271 

79* 

140 

CONTINUE 

00000690 

000310 

00273 

80* 


GO TO 300 

00000700 

000310 

00274 

81* 

160 

JJJ a 1 

00000720 

000312 

00275 

82* 


DO 230 U*i#NGO 

00000730 

000313 

00300 

63* 


i » mt 

06000740 

000332 



LZ~a 


00301 . 

84* 

00302 

85* 

00303 

66* 

00304 

67* 

00305 

SB* 

00306 

89* 

0031 1 

90* 

00312 

91* 

00313 

92* 

00314 

93* 

003 15 

94* 

00115 

95* 

003 16 

96* 

00320 

97* 

00 321 

98* 

00322 

99* 

00325 • 

100* 

O032« 

101* 

00326 

102* 

0032? 

103* 

00330 

104* 

00331 

105* 

00332 

106* 

00333 

107* 

00334 

108* 

00335 

109* 

00336 

no* 

00337 

lit* 

00342 

112* 

00343 

113* 

00344 

1 14* 

00345 

1 55* 

00346 

116* 

00347 ' 

117* 

00350 

1 18* 

OO3S0 

119* 

00351 

120* 

003S4 

121* 

00357 

122* 

00360 

123* 

00361 

124* 

00364 

125* 

00366 

126* 

00367 

1 27 * 

00370 

128* 

00370 

129* 

00370 

130* 

00371 

131* 

00372 

112* 

00373 

133* 

00374 

134* 

00375 

135* 

00376 

136* 

00377 

137* 

00401 

138* 

00402 

139* 

00403 

140* 


G6W8 • 0,0 

8UHI 8 0,Q 

80M2 > 0,0 

SUM3 * 0,0 
FN x FSDER(II) 

DU 170 IJ*JJJ,ND8C 
JJ = IJ 

CALL L(JCFAC(JKl,TOO(IJ3»T t CPL,NTCf»U)#Uffn 
HO s THLUU+FUCTHLCIUU-THLUU) 

CALL XVDHV<XV,DHV,ARK,8RK,PC,TOD(IJ),J$PC) 

CALL EQ5T A T (RHODj RG* DRLDT ,PCfTDO(IJ), XV»EMHL {J$PC ) » EHKCSi 
l ARK,6RK,JSPC) 

IFCWSPRCIJ) ,GT,FN*Shi5PR) GO TO 200 
GSwS 8 GSWSf*SPR<IJ> 

SUH1 s SUM! ♦ WSPRU J)*VQDCIJ) 

SUM2 a SUMS ♦ W§PH ( I J ) *HD 

SUMS « SUMi + W$PR{IJ)/( ¥00 ( I J ) **2*RH0D*D0D (IJ)**2> 

170 FN s FSOERCII) - GSWS/SwSPR 
GO TO 210 

200 OHS * fN*s«SPR 
GSWS s GSHS+OkS 
JJJ = JJ 

WSPH(JJ) e KSPR{JJ)*OWS 
sum s sum ♦ ows*voo(jj) 

SUM2 s SUM2 + OWS*H0 

SUMS s Su«3 ♦ OWS/(Y{JD(JJ)**2 *RhOD*OOD(JJ)**2> 

210 SUM2 r SUM2/GSW3 

DO 2 JO K s 1,NMJXZ 

JJ » JtK 

GwSPR(I,Jj) a SKSPR*FSDER(n 5*NELEM*F0MJX«) 

GvELOl C I , J J 5 = SUH1 /GSWS 

CALL LOCFAC C Jk 1 , SUMS, THl , NTCP( 1 ) , 1 1 ,F1) 

GTODlCUjj) S TTCPUU»l)+Fl*CTTCPL(IHl»i)»TTCPtCn»U} 

CALL XVOHVCXV,DHV,ARK,BRK f PC,GTO0i CI»JJ),JSPC) 

CALL EaSTAT(RHOD»RG»DRLDT,PC»GTUDl ( 1 , JJ) , XV,EMWL ( J8PC > f 
i 6MHCG,ARk,BRK, JSPC) 

230 GOl AD 1 ( 1 » JJ ) s SOKTfGSWS/CGVELOi (I,JJ)**2*RHODf»SiJM3)) 

300 DO 500 K * 1,NMIX2 
JJ * J+K 
GSWS s 0,0 
DO «00 I*l»NGO 
400 GS»S » GSKS+GhSPRCI# JJ) 

P(JJ) 3 PC 
£MR J 3 SMRG(JJ) 

CALL CGPROP(SMRG(jJ) f T0UJ)#EM«tJJ),GAH(JJ)#VISC*GASFUJJ) f 

1 O.0,Q,XOV,EKRl|PC,NGO r VUI»GT0DIU,JJ),GVEL01(l*JJ), 

2 GWSPRCl,JJ) F GSWS * TLI ) 

AREA 1 ( J J ) 8 FCHA*CFFMix(K)aWFE+FOMIXCK)aWOXE)/WTE*ACHAH 
PUS(JJ) « PC 
VUS< JJ 3 3 VCG 
VGAS(JJ) s VCG 

AGO(JJ) 3 Saf?T(32.2*GAH{JJ)*iS4S,*TO{JJ)/e«H{JJ)) 

TGAStJJ) 8 TOUJ>*a,-CGAM(JJ)m,)*0,5*<VGASUJ>/AGO(JJ))*A2) 
500 RHOG(JJ) « RHUGF CTGAS(JJ),PC;EMM(JJ)«2> 
fchah * fcham+fcha 

J « J+NMJXZ 
NOER s NDERtl 


00000750 

00033$ 

00000760 

000336 

00000770 

000337 

00000780 

000340 

00000790 

000341 

00000800 

000346 

00000810 

000346 

00000820 

000350 

00000630 

000362 

00000840 

000370 

.00000850 

000403 

00000660 

000403 

00000870 ’ 

000425 

00000880 

000432 

00000690 

000435 

00000900 

000441 

00000910 

000445 

000457 

00000940 

000464 

00000950 

&00466 

00000960 

000470 

00000970 

000472 

00000980 

000474 

00000990 

000500 

000 0 1 000 

000504 

00001010 

OOOSIO 

000523 

000525 

00001050 

000537 

00001060 

000542 

00001070 

000547 

00001080 

000551 

0000 1 U90 

000563 

00001100 

000571 

00001110 

000604 

00001120 

000604 

000626 


000647 

00001200 

000&56 

00001210 

000666 

00001220 

000674 

00001230 

000674 

00001240 

000677 

00301242 

000701 

00001250 

000703 

00001260 

000703 

00001270 

000703 

00001280 

000750 

00001290 

000761 

00001300 

0 00 JoS 

00001310 

000765 

00001320 

000766 

00001330 

OQIOOO 

001012 

00001360 

001027 

00001370 

001032 

00001380 

Q01O3S 
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00404 

141* 

00406 

142* 

00407 

1 43* 

00410 

144* 

0041 1 

145* 

00412 

' 146* 

00413 

147* 

00416 

146* 

00417 

l 49 * 

00421 

ISO* 

00422 

151* 

00423 

152* 

00424 

153* 

Oyu25 

154* 

00430 

155* 

00431 

1.56* 

00433 

157* 

00934 

156* 

00435 

159* 

00456 

160* 

00437 

161* 

00441 

162* 

00442 

163* 

00444 

164* 

00446 

165* 

00450 

166* 

00451 

167* 

00453 

1 68* 

00455 

169* 

00457 

170* 

00461 

171* 

00463 

172* 

00464 

173* 

00465 

174* 

00467 

175* 

00470 

176* 

00471 

177*. 

00474 

178* 

00976 

179* 

00477 

160* 

00500 

181* 

00501 

182* 

00502 

183* 

00503 

184* 

005(14 

165* 

00504 

186* 

00505 

187* 

oosto 

188* 

0051 1 

189* 

00513 

190* 

00515 

191* 

00517 

192* 

00520 

193* 

00521 

194* 

00522 

195* 

00523 

196* 


600 


610 

620 


630 


640 


700 

710 


C 


16 


12 


IF(NDER,U,IDER) CO TO 10 

8U«i * 0,0 

KT a 0,0 

PHIGH a 0,0 

PLOW * 0,0 

IC » 0 

OU 600 Id>J 

SUM1 s Sum*GASFl(P*PU) 
W T a. WT t OASPLdJ 
PCI a SUM1/KT 
GO TO 620 

PCI » CPLOW+PHIGH)/2 4 0 
SUM1 s o.O 
t>0 630 1 = 1 , J 


C i •l«4,*(PUS(I ) -PC X )*>3a,2/HHOG(I) 

IF(V0S(1)*VU8(I),LT,4,*C) GO TO 640 

VGASin S (VUS(IUSCiRT(VUSa)*VUS(n«4,*C})/2.0 

TG A S C I } c TU(I)*(l,»(GAH£l)»l,}*0,S*(VG6S{l)/AGO{I))**2) 

HHOGCI) = RHOGF(TGAS(n,PCI ,EMWCI},2) 

ARfcAUn s 144,*GASFL U ) / l VGAS U ) *RHOG C I ) ) 

SUH1 t sum + ARfcAHI) 

FA 5 SUHl/(FChA«*ACHAM) 
lFfABS(FA*l,J,tt, 0,001} GO TO 700 * 

IFfFA.LT. 1.0) PLOW a PCI 
IF (F A , GE , 1 , 0 ) PHIGH = PCI 
IC = iC+1 

iFUC.GT.bOJ GO TO 700 

IF { I C ,GE ,2 } GO TO 610 

IF (PHiGH, Lfc, 0.0, UK, PLOW, LE, 0,0) IC » 0 

IF (PLOW, LE, 0,0) PCI a PCI-1,0 

IF (PhiGH,LE,0,0) PCI * PCj+t.O 

GO TO 620 

PHIGH s pci 

IFCPLOW.GT.O.O) Go TO 610 

PCI s PC I -1,0 

GU TO 620 

DU 710 lsl,j 

AREA 1(1) a AREA! (D/FA 

NOT a NGO+1 

NST » J 

NGF a 1 

NASEG a J 

I WRITE a 1 

JORUMsO 

IPUnCHsO 

INPUT TO CICH THROUGH »STC 
READ (5# STC) 

CONTINUE 

iFUwKITE.EG.l) GO TO 12 

IFUDRUH.EG.I) GO TO 15 

IF (IPUNCH,EQ,1 ) CO TO 14 

GU TO 9999 

CONTINUE 

IFILE=6 

IWRITgaO 

rn tn tc 


00001390 

00001410 

00001420 

00001430 

00001440 

00001450 

00001460 

00001470 

00001500 

00001S10 

00001530 

00001570 

00001560 

00001590 

00001600 

00001610 

00001620 

00001630 

00001670 
00001680 
00001690 
00001700 
00001710 
00 172 

00001730 
0000} 740 
00001750 
00001760 
00001770 
00001780 
00001790 
00001800 
00001810 
00001630 

00001870 

00001680 

00001890 

00001900 


001040 
Q01043 
001044 
001045 
001046 
001047 
001054 
001054 
001037 
001063 
001065 
001067 
001073 
OU1073 
001100 
001106 
001114 
001130 
ooi i«i 

00115a 
001161 
001165 
001171 
001176 
001204 
001212 
001215 
00 1220 
001224 
001241 
001247 
001255 
001257 
001260 
001263 
001265 
001267 
001273 
001276 
001301 
001303 
001305 
001306 
001307 
001319 
001310 
001311 
001316 
001316 
001320 
001323 
001326 
901330 
001330 
nni'sn 
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00525 

00526 

00527 

00530 

00531 

00532 

00533 
0053# 

00535 

00536 

00537 
00560 
0o5.at 
00502' 
00503 
00544 
00505 
00546 
00507 

00550 

00551 

00552 

00553 
0,055 4 
00555 
00560 
00562 

00565 

00566 

00567 
005 70 
00571 
00573 
00573 
COSTS 
00600 
00601 
00602 

00603 

00604 

00605 

00606 
00610 
00610 
00611 

00613 

00614 
00617 
00622 
00623 
00627 
00630 
00637 
00640 
00647 
00650 
00657 
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199* 

5 

IF JLE«t 1 

001334 

199* 


tO»UH*0 

001335 

200* 


GO TU 15 

001336 

201* 

14 

CONTINUE 

00I3«0 

202* • 


IF ILE*7 

001340 

203* 


tPUNCH*0 

00 1 34 r 

204* 

IS 

CONTINUE 

001143 

205* 


ILISPsO 

001343 

206* 


ISTC*1 

001 343 

207* 


itrans*o 

001 345 

208* 


ITDK=0 

001346 

209* 


I ST= 1 

001347 

210* 


•TBF=0,O 

001350 

211* 


T N8F *0*0 

001351 

212* 


SHQN9F*0 , 0 

001’352 

213* 


RHUtFtO.O 

001353 

214* 


wtmluf=o,o 

001354 

215* 


WTMIVF=0,0 

001355 

216* 


NOZONsO 

001356 

217* 


NSTPZ=0 

001357 

218* 


NUG=0 

001360 

219* 


IPUN30=0 

005361 

220* 


OHVFiO.O 

001362 

221* 


zstart=xhinde 

001 363 

222* 


DU 31 I*t,NTAB 

001404 

223* 

31 

CStRU)=CSTAR(TM«(I)»TTO<l)»TSA«U}J 

001405 

224* . 


DO 715 J=1»NTAB 

001441 

225* 


GMH f 1 , JlsTMRtJ) 

001441 

226* 


GMRC2, J)eTTO(J) 

001442 

227* 


GMR{3,J)*TVIS(J) 

001444 

228 * 


GHR(4»J)*TGAM(J) 

001446 

229* 


GMRC5, J)=THWCJ) 

001450 

230* 


GMRC6, J)eCSTR(J) 

00 1452 

231* 

715 

continue 

001461 

232* 


DO 716 J*l,NTK 

001461 

233* 


GTK(t,J)*0,0 

001461 

234* . 


GTK(2,4)=0,0 

001461 

235* 


GTKC3»J)*0,0 . 

001462 

236* 


GTK(4,J)sTV0(J) 

001463 

237* 


GTHI5»J)sQPvO(Of) 

001465 

238* 


GTK(6,J)sTC0NV0(J) 

001467 

239* 

716 

CONTINUE 

001472 

240* 

9000 

FURMAT(lHt f ///,z3X#‘ClCH GENERATED INPUT DAT* FOR DER SUBPROGRAM S 

001472 

2« 1 * 


*tc *t/n 

001472 

242* 


WRITE (6*9000) 

001472 

243* 


jCAROsl 0 

001477 

244* 


WHITE CIFJLE, 7} jCARD 

001501 

245* 


DO 720 1*2,4 

001512 

246* 


ICARD* J *1 0 

001512 

247* 

720 

WRITEUFlLfcrB) JCARD 

001515 

248* 


ICAROs50 

001525 

249* 


WRITE tIFILE,!) |LISP»18TC» I TRANS, I TDK, ICARD 

001527 

250* 


ICARD*5010 

001541 

251* 


MRITE(|FILE,1) N0Z0N,NSTP2»NUG,IPUN30, ICARD 

00 1543 

252* 


|CARD*5020 

001555 

253* 


WRITE ( I FILE, 1 ) NP,NCHAMC,NTAB,NTK,XCAR0 

001557 

259* 


DO 725 I*1,NCHAMC 

001574 



QQbti 

255* 

7 

ACHAMCm»58RT<«,/3.!«159-i.‘CHAKCm> 

001578 

00**3 ‘ 

256* 

725 

CONTINUE 

001603 

00665 

257* 


xnap*nchamC 

001603 

00666 

258* 


NC*XNAP/3,t,9 

00160b 

00667 

259* 


DU 730 I a 1 » NC 

001620 

00672 

260* 


lC.AHD*5020 + J*10 

00162* 

00673 

261* 


M*|*3«2 

001630 

0067* 

262* 


N=m + 2 

00163* 

00675 

263* 

* 730 

WRITE (IT 1U#5) ( XCH AHC ( J ) , ACHAHC ( J ) , J*H, N ) , JC7RD 

001636 

00706 

26** 


ICARD*5100 

001660 

00707 

265* 


WRIUUFU.E,9) JCARO 

001662 

00712 

2 66* 


xnmr=ntar 

001670 

007 13 

267* 


■K!C = XNMR/fc,*,9 

001673 

0071* 

266* 


no 7*0 1=1,6 . 

001705 

00717 

269* 


DO 7*0 J= 1 , NC 

001721 

00722 

270* 


ICARD=5000*I*100+J*iO 

001721 

00723 

271* 


M=0 *6»5 

001725 

0072* 

272* 


N=6+5 

001731 

00725 

273* 

7*0 

•hRlTUlFjUrS) (GMRU,U#l»H,N),ICARD 

001736 

00736 

27** 


xntkcntk 

001762 

00737 

275* 


NC=XNTK76,+,9 

001765 

007*0 

276* 


DO 750 1=1,6 

001777 

007*3 

277* 


DU 750 J=1 , NC 

002013 

007*6 

276* 


ICARD=5600+I*100+U*10 

002013 

007*7 

279* 


H=j*b»5 

002017 

00750 

260* 


N = M+5 

002023 

00751 

281* 

750 

HKUKU 1LE,S) C5TKU,L),l*M»N>,ICARO 

.002030 

00762 

262* 


ICARD=65iO 

00205* 

00763 

203* 


WRITE (I FILE, 5) TNBF,TBF,RH0NBF,RH0LF,EHWL(21,EMWV(2),ICARD 

002056 

0077* 

26** 


1 C AROc6520 

002072 

00775 

265* 


WRJ Tt ( IFjlEjS) TNBU,TBU,RHONBU,RHOLU,EMWL(n,EHWV(l),ICARD 

00207* 

01006 

266* 


ICARD=fe530 

002110 

01007 

287* 


WRITE (IF HE, 3) TCRIT(2),TCRIT{1),OHVF,OHVO,1CARO 

002112 

01016 

286 * 


ICAHD=6S*0 

00212* 

o ro 1 7 

289* 


WHITE (I FILE, 6) lST,NSSTI ,NMST I , ICRC , IPRSST , IPRHST, I CARD 

002126 

01030 

290* 


ICARD=6550 

0021*2 

01031 

291 * 


WRITE t IF I LE, It) CrTOL,AHTOLD,ICARI) 

0021** 

0-10 36 

292* 


tCAHD=67)0 

0021b* 

01037 

293* 


WRITEUFJLE, 1J) PCI , Z9T ART, ICARO 

002156 

010** 

29** 


ICAHD66720 

002166 

010*5 

295* 


write cifile,!) ngt,nsf;nst,naseg,icard 

002170 

0105* 

296* 


XNGfaNGT ! 

002202 

01055 

297* 


NC = XN{;T/2, + ,9 

- 002205 

01056 

298* 


DO 800 JslfNST 1 

002217 

01061 

299* 


ICARD=7000*J*10 

00222* 

01062 

300* 


WR1TEUFILE,2) AREA 1 ( J) ,SA3FI ( J ) , 3MRG (J) , ICARD 

00223Q 

01070 

301* 


DU 600 1*1, NC 

0022*4 

01073 

302* 


ICARDalCARD+l 

0022** 

0107* 

303* 


H=I*2«1 

0022*7 

01075 

30** 


N=M+1 

002253 

01076 

305* 

800 

WRITE (IFIUE, 53 (GNSPR(L» J)»G¥Ei.Dl (L,J),GDIA0i (L,'*) ICARO 

002260 

oim 

30b* 


GO TO 16 

002306 

01112 

307* 

9999 

CONTINUE 

002310 

OIU3 

308* 


RETURN 

08002170 002310 

01118 

309* 


END 

ooonstRA ft a * ia * 



TABLE B-III CARD CHANGES TO CICM ROUTINE 
DERINI FOR IMPROVED STC INTERFACE 


9FOR , URS 

*7 


DERINI 


>11 


01 MENS ION GMR(6, 1 8 ) , GTK (6 , 20 ) , CSTR U8 ) , T VO (24 ) ,CPVO (24 ) , TCON VO C24 ) 
C0MM0N/CGTA6C/NTA8,STT ? AMRT,TMRi(l8),TT0(18),TGAM(18) f TMW(|8), 

* TVIS(IS) 

COMMON/CHCOM/NCHAMC,M2C,NCON4C,WGJC,EMRGJC,STGJC,.EMrtGJC,GAMGJC» 

* X LMC,DEITXC,BSPRC, CSPRC, ACSC , CLN TC , CUNRAC ,CCANGC, 

* R CBCC,RCTC,XCHAMC(203,ACHAMC(20) 

*13, 17 

1 F0RMAT(4I12,24X,I8) 

2 F0RMATUE12,6,36X,I8) 

3 FORMAT (4E12,6,24X, 18) 

4 FORMAT (1216, 18) 

5 FORMAT (6E 1 2 ( 6 , 18) 

6 FORMAT (6.1 1 2, 1 8 ) 

7 FORMAT (6X, »STC INPUT FROM CICM PROGRAM C ASE ' , T73 , I 8) 

8 FORMAT ( 72X , 18 ) 

9 F0RMAT(6X, ‘02/H2 GAS PROPERTIES FROM CICM INPUT •, T73, 18) 

11 FORMAT (2El2 f 6,48X, 18) 

17 FORMAT ( 1216 ) 

18 FORMAT (6E 12 ,6) 

NAMELIST/STC/TV0,TC0NVO,NTK,TN8O,T8G,RHON8O,RHOLO p NP,NSSTI,NMSTI P 

* ICRC,IPRSST, IPRMST ,CRTOt , ARTOLD, IWRITE, IDRUM, IPUNCH, 

* CPVQ,DHVQ 

*/CGAM ( 1*) ) )?/GAM , * SQRT (a9<S,77 ‘* GAM * T0/XMw/( C2,/(GAM + 1 ,))**{ (GAM*S 0 ) 

*20,20 

-22,27 

10 READ(IRDER, 17) 

READ ( IRDER , 18) 

READd-RDER, 18) 

REAOdRDER, 17) 

READ(IRDER,18) 

READflROER, 16) 

* 165,187 . 

iWRITEsl 
. IORUMsO 
IPUNCH=0 

C INPUT TO CICM THROUGH $STC 

REA0(5,STC) 

16 CONTINUE 

IF<IWRITE,EQ,1) GO TO 12 
IFdDRUM.EG, 1 ) GO TO 13 
IF ( I PUNCH , EG f 1 ) GO TO 14 
GO TO 9999 

12 CONTINUE 
IF ILE=6 
IwRlTEaO 
GU TO 15 

13 CONTINUE 
IFILEstl 
iDRUMaQ 
GO TO 15 

14 CONTINUE 


NMIXZ,NG0 

(FFMlX ( I ) , FOMI X ( I ) , I * 1 , NMJXZ ) 

(FSDER ( I ) , lx 1 , NGO) 

NDSc,NELEM 

( WSPR ( I ) , DOD ( I ) , TOD C I ) , VOD ( I ) I la 1 , 1 00 ) 
PC,WCG,EMRCG,SWSPR,VCG,VLJI,TLI,fcMWCG,FCHA 
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IFILEs? 

IPUMCHsO 

IS CONTINUE 
ILISPbO 
iSTCsl 
ITRANSaO 
ITOKsO 
1 8 T s 1 
TSFaO.O 
' TNBFsO.O 
RhONBFsO , 0 
RHOUFsOiO 
WTMtLFsd.O 
WTMLVFsO.O 
NOZON=0 
NSTPZsO 
NUG»0 
IPUN3D30 
OHV.FaO.'O 
ZSTARTaXMlNDE 
DO 31 I=1,N7AB 

31 CSTRU)sCSTAR(TMWU)*TTO(J),TGAM(I}) 

DO 715- J=1,NTA6 
GMR(1 # J)sTMR(JJ 
GMR(2,J)sTTO( J) 

GMR(3, J)aTVIS(J) 

GMR(4, J)=TGAM(J) 

GMR (5,J)sTMW(J) 

GHR(6, j)sCSTR(J) 

715 CONTINUE 

DO 716 J*1,NTK 
GTK(l,J)aO.O 
GTKt2,J)s0.0 
GTK(3,J)s0.0 . 

GTK(4,J)=TVO(J) 

. GTK(5#.J)aCPVO(J) 

GTK (6, J)aT.CUNVO(J) 

716 CONTINUE' 

9000 FORM A T(lHl,///, 2 Jx,fClCM GENERATED INPUT DATA FOR DER SUBPROGRAM S 
*TC ’ ///) 

WRITE (6,9000) 
iCARDsl.o 

WRITE{IFILE,7) ICARD 
DO 720 1=2/4 
ICARDsI *10 

720 WRITE(IFILE/B) ICARD 
ICARDsSO 

WRITE (IF I LE, 1) II ISP, ISTC,ITRaNS/ITdK,ICARD 
ICARDbSOIO 

WRITE (IFILE, 1) N0Z0N,NSTPZ,NUG,IPUN3D, ICARD 
ICARO»5020 

WRITE (IFILE, 1) NP,NCHAMC, NT AB,NTK. ICARD 
DO 725 I«1,NCHAMC 

ACHAMC(I) = SQRT (4 t /3 f 141'59*ACHAMC (D) 

725 CONTINUE 
XNApBNCHAMC 
NC*XNAP/3, ♦ ,9 
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OQ 730 Isl,NC 
* ICARD=5020+ I *1 0 

Mb I *3«2 
MaM+2 

730 HR ifE Cl FILE #5) (XCHAMC ( J) , ACHAMC ( J) , J=m,n) , ICARO 
ICARO=SlOO 

WRITE (IEILE,9) ICARD 
XNMRsNTAB 
NCsXNMR/6,+,9 
00 740 1=1, to 
• 00 740 J*1,NC 

rCARD=5000+I *100+J*10 
MsJ*6»S 
NsH+5 

740 WRJT.E(IFILE,5) (GMR ( I , L ) , L=M, N ) , Ic ARD 
XNTKsN.TK 
NCbXNTK/ 6. + ,'9 
00 750 1*1,6 
00 750 J=1,NC 
ICARD*5600 + IM00+J*10 
Ms J *6<*5 
N=M + 5 

750 ‘WRITE(IFILE,5) 

ICARDs'6510 
WRIT£(IFILE,5) 

I.CARD=6520 
WRIT£(IFILE,5) 

ICARD=6530 
WRI.TE ( IFILE , 3) 

ICARO=6540 
WRITEtlFILErfel 
ICARO=6550 
WRITE (If IlE, 1 1) CRTGL,ARTOLD, ICARD 
ICARDtsfe.7 1 0 

WRITEUFILE, uj pci, zstart, icard 
ICARD= b720 

WRITE (IFILE, 1) NGT,NGF,NST>NASEG, ICARD 

XNGTsNGT 

NC»XNGT/2,+.9 

DO 800 J = 1 , N S T 

ICARD*7000+J*10 

WRITE (IFILE #2) AREA 1 (J ), GASFL (J >, SMRG (J ), ICARO 

DO 800 1=1, NC 

ICARDsICARD+1 

Ms I *2® 1 

NsM+l 

800 WRITE (IFILE, 5) (GWSPR (L , J ), GVELD 1 (L , J ), GDI ADI (L , J ) ,L*M, N) , ICARD 
GO TO 16 
9999 CONTINUE 


(GTK(I.,L),L=M,N), ICARD 
TNBF,TBF,RH0NBF,RH0LF,EMWL(2),EMHV(2),ICARD 
TNBO, TBO,RHONBO, RHOLO,EMWL ( 1 ) , EMWV ( i ) , ICARD 
TCRlT<2),TCRlT(i),DHVF,DHVG, ICARD 

IST,NSSTI ,NMST-I , ICRC, IPRSST, IPRMST, ICARD 
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TABLE B-IV 

NAMELTST INPUT VARIABLES FOR 
IMPROVED CICM/STC INTERFACE ROUTINE 


VARIABLE NAME 

DEFINITION 

UNITS 

IWRITE* 

STC input data generated by CICM will be printed 
out when IWRITE - 1 

- 

I DRUM* 

t 

STC input data generated by CICM will be written 
on system drum file 11 when IDRUM = 1 


I PUNCH* 

- STC input data generated by CICM will be punched 
on cards when I PUNCH = 1 


NP 

- Total number of z-p lanes between z=ZSTART and 
nozzle throat 


NSSTT 

Maximum number of complete passes, marching from 
z = ZSTART to throat, in single tube analysis 

— 

NMSTI 

Maximum number of passes in multiple stream 
tube analysis 

_ 

ICRC 

■Number of corrector cycles calculated at 
each az interval 


IPRSST 

Number of Az intervals between single stream tube' 
printouts 

_ 

IPRMST 

Number of az intervals between multiple stream 
tube printouts 


CRTOL 

.Decimal tolerance, deviation of computed single 
stream tube throat contraction ratio from unity 

— 

ARTOLD - 

Decimal tolerance, deviation of computed 
multiple stream tube throat contraction ratio 
from unity • 

- 

NTK 

Number of temperatures at which propellant vapor 
specific heats and film thermal conductivity are 
tabulated 

~ 

TVO (20) 

Temperatures at which oxidizer CPVO and TC0NV0 
are tabulated 

°R 

CP VO (20) 

Oxidizer vapor specific heat at constant pressure 

Btu/lbm-°R 

TGONVO (20) 

Thermal conductivity of vapor/gas film surrounding 
oxidizer droplets 

Btu/ft-sec-°R 

DHVO 

Oxidizer latent heat of vaporization at chamber 
"wet bulb" temperature calculated by CICM 

Btu/lbm 

TNB.O 

Oxiidzer normal boiling point 

°R 


B-28 



VARIABLE NAME 


DEFINITION 


UNITS 


TBO 

RHONBO 

RHOLO 


NOTE: 


Oxidizer droplet saturation temperature at Pc 

Oxidizer density at normal boiling point 

Oxidizer density at saturation temperature 
corresponding to Pc 


°R 

Ibm/in 3 

lbm/in 3 


All parameters except those asterisked are identical to descriptions 
given in the DER users manual (Ref. 2 ). 


B-29 



TABLE B-V NAMELIST INPUT FOR MODIFIED 
C1CM/STC INTERFACE SUBROUTINE 


SSTC 

CRTOL»0,01, ARTt)LD = 0,01, IWRITEsl, IORUMsO, JPUNCHcO, 

NPb 50, NSSTI=3, NMSTla3 # K*RC B 1# IPRSST*25, IPRM$T»25, 
TNB0*162,, TB0=265,, RHoN6os, 04 1 3# RHOtOs,0271, DHV0 B 45., 
NTK«20, - “ 

T VO (1)«200,-, 265., 275,, 285,, 300, ,340., 400,, 600,, 1200 ,,1800,, 

2400.. 3000.. 3400.. 3800.. 4200.. 4600.. 5000.. 5600., 

6000. . 6400., 

CP VO ( 1 ) s • 94 , , 94 , ,55, , 43 , ,356, ,286, 0 257, ,226, , 245 , ,260, ,269, 
,276,. 280, .284, ,288, ,292,, 2955, ,301, ,304, ,307, 
TCOMVOd )=, 000009 17, ,00000917, .0000108, ,000011057,00001130, 
,0 000! 167, ,00001389, ,00001806, ,00002778, ,000035, 
,00004028, ,00004444, ,00004583, ,00004681 , ,00004722, 
,00004639, ,000045, ,000 03806, ,0 000191 7, ;0, 

SEND 
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TABLE B-VI Cl CM SAMPLE CASE GENERATES INPUT ELEMENT FOR STC 


9ELT+L STC 


03 

I 

CO 


UNIVAC 

1. 

2. 

5, 

4. 

5. 

6. 

7, 

0. 

1100 SERIES ELT PROCESSOR 
STC INPUT PROM CICM 

0 i‘ 

0 0 

50 2 

,000000+00 ,963300+01 

LEVEL MB AT 8 J 52 U 0 AM 
PROGRAM CASE 

0 0 

0 0 

16 20 

,500000+01 .255300+01 

ON THURSDAY , 

FEBRUARY 26+ 

1976 

10 

20 

30 

40 

50 

5010 

5020 

5030 


Q 2/ H 2 I 

5*3 PROPERTIES FROM CICM 

INPUT 



5100 

10. 

,000000 

,100000+00 

,500000+00 

,100000+01 

,150000+01 

, 200000 * 0 : 

5110 

n. 

,250000+01 

,300000+01 

,350000+01 

,400000+01 

,450000+01 

,500000+01 

5120 

12. 

,550000+01 

,600000+0! 

,700000+01 

,800000+01 

,000000 

,000000 

5130 

1 3 * 

,590000+03 

,721514+03 

,143978+04 

,226954+04 

.300 B 16+04 

,366881+04 

5210 

i«. 

,925666+09 

,477364+04 

,520902+04 

.556641+04 

,585309+04 

,607691+04 

5220 

15. 

,629526+09 

.636535+04 

,690942+04 

,650606+04 

,000000 

,000000 

5230 

15. 

,599022-05 

.751820-05 

, l 33 b 2 7-04 

,209151-04 

,201415-04 

,349234-04 

5310 

1 7 . 

,911210-09 

,466013-04 

,514907-04 

,555499-04 

,509147-04 

,616405-04 

5320 

is. 

,630106-09 

,654000-04 

,675904-04 

,666049-04 

,000000 

,000000 

5330 

19. 

, 1 90 S 00 + 01 

,139700+01 

,138000+01 

,133800+01 

.130100+01 

,127400+01 

541 0 

20. 

, 125900 + 01 

,124000+01 

,122900+01 

.122000+01 

,121400+01 

,120900+01 

5420 

21. 

,120500+01 

,120200+01 

, l 19900 + 01 

, 1 19600 + 01 

,000000 

,000000 

5430 

22. 

,201600+01 

,221800+01 

,302400+01 

,403200+01 

,504000+01 

,604700+01 

5510 

23. 

,705000+01 

.004000+01 

,900700+01 

,994100+01 

,108370+02 

, 116920+02 

5520 

29, 

,125020+02 

,132660+02 

,146520+02 

,150620+02 

,000000 

,000000 

5530 

25, 

,532077+09 

, S 0833 S +04 

,713601+04 

,704513+04 

,015031+04 

,828650+04 

5610 

20, 

,031991+09 

,027652+04 

,019444+04 

,600436+04 

,795305+04 

,701410+04 

5620 

27, 

,766900+09 

,752364+04 

, 723406 + 04 

,696488+04 

.000000 

,000000 

5630 

26, 

,000000 

,000000 

, 000000 

,000000 

,000000 

,000000 

5710 

29. 

,000000 

, OOOQOO 

,000000 

,000000 

,000000 

,000000 

5720 

30, 

.000000 

,000000 

,000000 

,000000 

,000000 

,000000 

5730 

31, 

,000000 

,000000 

,000000 

,000000 

,000000 

,000000 

5740 

32, 

,000000 

,000000 

,000000 

,000000 

,000000 

,000000 

50(0 

33, 

,000000 

,000000 

.000000 

,000000 

,000000 

,000000 

5820 

39, 

,000000 

,000000 

,000000 

,000000 

,000000 

,000000 

5850 

35, 

,000000 

,000000 

,000000 

,000000 

,000000 

,000000 

5840 

30, 

,000000 

,000000 

,000000 

,000000 

,000000 

,000000 

5910 

37, 

,000000 

,000000 

,000000 

,000000 

,000000 

,000000 

5920 

30, 

,000000 

,000000 

,000000 

,000000 

,000000 

,000000 

5930 

39, 

,000000 

,000000 

,000000 

,000000 

,000000 

,000000 

5940 

90, 

,200000+03 

,265000+03 

,276000+03 

,285000+03 

,300000+03 

,340000+03 

6010 

91, 

,400000+03 

,600000+03 

, 120000+04 

,180000+04 

,240000+04 

,300000+04 

6020 

92 , 
93, 

,340000+04 

.000000+04 

,300000+04 

,640000+04 

,420000+04 

,460000+04 

,500000+04 

,560000+04 

6030 

6040 

6110 

99. 

.940000+00 

,940000+00 

,550000+00 

,430000+00 

,356000+00 

,266000+00 

95, 

,257000+00 

,226000+00 

,245000+00 

,260000+00 

,269000+00 

,276000+00 

6120 

90, 

97. 

,280000+00 

,304000+00 

,204000+00 

,307000+00 

,208000+00 

,292000+00 

,295500+00 

,301000+00 

6130 

6140 

90, 

,917000-05 

,917000-05 

, 108000-04 

,110500-04 

,113000-04 

,116700-04 

6210 

99, 

,138900-04 

,180600-04 

,277800-04 

,350000-04 

,402000-04 

,444400-04 

6220 

50 , 

51, 

,458300-04 

,191700-04 

,468100*04 

,000000 

,472200+04 

,463900-04 

,450000-04 

,380600-04 

6230 

6240 

52, 

,000000 

,000000 

.000000 

,000000 ■ 

.201600+01 

,201600+01 

6510 

53, 

59, 

,162000+03 

,590000+02 

,265000+03 

,270600+03 

,413000-01 

,000000 

,271 000-01 
,450000+02 

.,320000+02 

,320000+02 

6520 

6530 


(CYCLE 



55. 

1 

3 

3 

56, 

.tooooo-ot 

,100000-01 


■57, 

,701007m 

,150000901 


58. 

12 

1 

0 

59. 

,289392+01 

,180291+01 

,235010+01 

60, 

,000000 

,100000+03 

,000000 

61. 

,167655+00 

,338165+03 

,797571-02 

62, 

, 167655 + 00 

.2951OJ+03 

,108685-01 

63. 

, t67655 + 00 

.250306+03 

,108369-01 

60. 

,167655+00 

,205803+03 

,195075-01 

65, 

,836276-01 

,180150+03 

,153935-01 

66, 

,330055+01 

,213035+01 

,287724+01 

67, 

.000000 

, 100000+03 

‘,000000 • 

68, 

,200912+00 

,338165+03 

,797571-02 

69. 

,200912+00 

, 2951*01+03 

,108685-01 

70, 

,200912+00 

,250306+03 

,106369-01 

71. 

,200912+00 

,205801+03 

,195075-01 

72. 

,102056+00 

,180150+03 

,153935-01 

73. 

,332126+01 

,220519+01 

,318061+01 

70, 

,000000 

, 100000+03 

, 000000 

75. 

.228233+00 

,330501+03 

,797777-02 

76. 

,228233+00 

,291000+03 

,106120-01 

77. 

,228233+00 

,206802+03 

,106856-01 

78. 

,228233+00 

.201528+03 

,190719-01 

79. 

.110116+00 

,175890+03 

,158903-01 

80. 

,382225+01 

,206868+01 

,212001+01 

81. 

,000000 

, 100000+03 

,000000 

62, 

,228233+00 

,330501+03 

,797777-02 

03; 

,228233+00 

,291000+03 

,108120-01 

80, 

,228233+00 

,206802+03 

,106856-01 

85. 

,228233+00 

,201528+03 

,190719-01 

86, 

,110116+00 

,175800+03 

,158903-01 


END ElT, TIMEl 0,3866 SECONDS , 


8FIN 
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,1 

25 

25 

6540 

t 



6550 

6710 

6720 

7010 

,167655+00 

.atot^+os 

.S53597-02 

7011 

,167655+00 

,317580+03 

, 921000-02 

7012 

,167655+00 

,271780+03 

,128331-01 

7013 

,167655+00 

,227906+03 

,171472*01 

7014 

,167655+00 

.169710+03 

.193210-01 

7015 

,838276*01 

,162630+03 

,126132-0! 

7016 

,200912+00 

,369139+03 

,553597-02 

7020 

7021 

,200912+00 

,317580+03 

,921000-02 

7022 

,200912+00 

,271780+03 

,128331-01 

7023 

,200912+00 

,227906+03 

,171072-01 

7024 

.200912+00 

,189710+03 

,193210-01 

7025 

,102056+00 

,162630+03 

,126132-01 

7026 

,228233+00 

,367607+03 

,518475-02 

7030 

7031 

,228233+00 

,310769+03 

,911975-02 

7032 

,228233+00 

,268513+03 

,126896-01 

7033 

,228233+00 

,223920+-03 

,170021-01 

7034 

,220233+00 

,183100+03 

,197896-01 

7035 

1 1 0 1 1 o+ 00 

,152905+03 

,126287-01 

7036 

,228233+00 

,367607+03 

,518075-02 

7000 

7001 

,228233+00 

,310769+03 

,911975-02 

7002 

,228233+00 

,268513+03 

,128696-01 

7003 

.228233+00 

,223920+03 

,170021-0! 

7004 

,228233+00 

,183100+03 

, 197896-01 

7045 

.1101J6+00 

,152905+03 

,126267-01 

7046 



Criteria for Specifying the CICM/STC Inter face 
Plane Location 


In Section 2.1.2 of CPIA 246 it is recommended that for 
CICM/STC analyses the CICM program should be executed to the axial plane at 
which the liquid jet has disappeared for all flow zones. At this point CICM 
output is transferred into STC input. There are two problems with specifying 
the interface plane in this manner. 

(1) The CICM program contains an advanced droplet heat- 
up and vaporization model.- The subcritical K-Prime STC version assumes' a 
constant "wet bulb" propellant temperature. If CICM execution is limited to 
the point of liquid jet dissipation a significant percentage of the liquid 
droplets will not have yet heated to the "wet bulb" temperature. It is phy- 
sically incorrect to ignore this effect and to characterize all the liquid 
droplets with a constant temperature and latent heat of 'vaporization in the 
STC input. 


(2) CICM performance predictions are controlled, in 
large part, by two empirically correlated input constants, C A and B^. C ft is 
an atomization {jet stripping) rate constant and is a drop siz e constant. 

The recommended input values for these coefficients were backed out from CICM 
by correlating hot test data. In these instances, CICM was allowed to compute 
to the chamber throat plane. Thus, the technique used to derive the constant 
input values is inconsistent with the recommended procedure of joining the CICM 
. and STC analyses at an intermediate chamber axial plane. 

■ CICM improves the JANNAF methodology for subcritical 
propellants because it allows for the droplet temperature transient. However, 
it is economically unrealistic to use CICM to the chamber throat plane 
because of the high computation time this technique requires. Also, using CICM 
to the throat plane would cause the coaxial injector analysis technique to be 
inconsistent with the JANNAF conventional liquid/liquid injector methodology, 
which utilizes STC. The most physically realistic technique is to have CICM 
execute until all the calculated oxidizer drop size groups have been heated 
to the chamber "wet bulb" temperature. As previously cited, for oxygen the 
unsteady state typically comprises only approximately 10 percent of the total 
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time required to vaporize 99 percent of the propellant. Thus, STC would 
still be responsible for calculating the majority of the liquid mass transfer 
to the gaseous phase. Importantly, the STC assumption of constant liquid drop 
temperature is verified when CICM calculates the complete unsteady state time 
period . 


5. Specification of Intra-Element and Manifold Zone Mass 

Distributions 

There is one additional technical problem in interfacing 
the CICM and STC programs. CICM does not contain formulations for calculating 
intra or inter-element mixing. The subjects were previously discussed in’ 
Sections B.6-7 of this appendix. CPIA 246 and the CICM user's manual recommend 
the following two solutions. 

(1) Manifold mass maldistributions should be accounted 
for by modeling separate chamber flow field zones. 

(2) Intra-element mass maldistributions are modeled by 
using empirical -single element cold flow data to input distinct radial mass 
distribution sub-zones to CICM for each chamber flow field zone designated 
as described in (1) above. 

There are at least the following four limitations to these suggested solution 
techniques. 

0) The JANNAF programs can not allow for the dissipation, 
due to diffusion mixing, of the face plane measured manifold distributions. 

(2) The JANNAF methodology does not recommend where the 
single element mass and mixture ratio distribution should be specified. If 
the distribution is measured at the face plane the solution will be in error 
because coaxial elements rely on shear (gas/liquid a V) mixing to produce nearly 
uniform mass distribution at the chamber throat plane. 

(3) The CICM and DER literature list only one example of 
application of the .recommended coaxial mass distribution specification technique 
(the J-2S sample case that is included in CPIA 246). The method that was used 
to specify the given flow distribution is not described. However, it was stated 



that the given distribution was known to result in low performance predictions. 
This would be expected if the given distribution did not account for shear 
mixing to the chamber throat, plane. 

(4) As previously cited in Section B.6 of this appendix, 
the test cases used to back out the recommended atomization and drop size inputs 
to DER assumed that the thrust chambers in question had uniform throat plane 
mixture ratio distributions. For most real coaxial injectors there will be a 
finite mixing loss because the coaxial element 'is a relatively slow mixing 
element. It is apparent that the correct values for the C A and E$ A coefficients 
will be directly dependent on the assumed single element mixture ratio dis- 
tribution. Unless a standard method for measuring or calculating single 
element and chamber mixture ratio distributions is developed it is extremely 
doubtful that universal values for the and constants can be verified. 
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APPENDIX C 


JANNAF SIMPLIFIED PREDICTION 
PROCEDURE FOR CICM ANALYSIS 




The M-l sea-level, pressure fed facility for ablative chamber 
testing is shown in Figure C-l . The corresponding instrumentation code sheet 
follows Figure 1. Figure C-2 specifies chamber pressure tap axial and cir- 
cumferential locations. 

The M-l injector design layout is shown in Figure C-3. The 
injector contained 3248 coaxial elements with gaseous hydrogen being injected 
annularly around the oxidizer. A row of orifices, drilled through the porous 
face, was located around the injector periphery and provided the chamber wall 
film cooling. Approximately 3.7% of the total fuel flow rate was used for 
chamber wall film cooling. Total fuel element flow rate was 89.8% of the thrust 
chamber fuel flow rate with a baffle fuel film cooling flow rate of 3.9%. The 
remaining 2.6% of the fuel flowed through the rigimesh injector face. The 
element consisted of two basic components which were threaded together. An 
oxidizer tube was recessed within the fuel sleeve producing a fuel annulus 
between the two parts. The fuel annulus was fed by four holes having an area 
four times that of fuel annulus. The oxidizer tube was flared at a fifteen 
degree included angle and was recessed 0.231 inches from the injector face. 
Elements were arrayed in 33 concentric rows. 

The low area ratio combustion chamber used for t esti ng with the 
M-l injector is comprised of an outer steel shell and an inner ablative liner. 
The assembled combustion chamber (Figure C-4) consists of an upper fuel torous 
and a lower conical combustion chamber. The thrust chamber design parameters, 
as related to the ODK input parameters, are identified in Figure C-5. 

The test 009 nominal computer input decks for the JANNAF programs 
utilized during the M-l analysis are shown in Figures C-6 through C-9. 
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M-l INSTRUMENTATION TAP LOCATIONS 


Measurement 

Oxid. 

Tap 

Loc. 

Fuel 

Tap 

Loc. 

Tank Pressure (POT, PFT) 

01 

FI 

Flow Meter Pressure (POFM, PFFM) 

- ■ 

F2 

Flow Meter Temperature (TOFM, TFFM) 

| 

F4 

GH 2 Mixer Pressure (PFMIX-2) 


F5 

Thrust Chamber Pressure-1 (POTCV-1, PFTCV-1) 

06 

F6 

Thrust Chamber Pressure-2 (POTCV-2, POTCV-2) 

07 

F7 

Thrust Chamber Temperature (TOTCV-2, TFTCV-2) 

08 

F8 

Thrust Chamber Injector Pressure (POJ, PFJ) 

09 

F9 

Thrust Chamber Injector Temperature (TOJ, TFJ) 

010 

F10 

Thrust Chamber Pressure (Pc4B-l & 2) 

11 

11 


FIGURE T,-I .(cont. ) INSTRUMENTATION CODE 
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Row 

No. 

Number of Orifices' 
Req'd per Row 

,1 

10 

2 

16 

3 

24 

<1 

30 

5 

36 

6 

42 

7 

48 

8 

48 

9 

48 

10 

66 

11 

72 

12 

76 

13 

84 

14 

90 

15 

84 

16 

84 

17 

96 

18 

96 

19 

108 

20 

108 

21 

120 

22 

120 

23 

132 

24 

132 

25 

144 

26 

144 

27 

156 

28 

156 

29 

168 

30 

168 

31 

180 

32 

180 

33 

180 


FIGURE C-3. M-l INJECTOR DESIGN 
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Injector 

Face 


38.31 Dla. 


23.96 


30.00 Dia. 
Throat 


43.1 Dia. 


40.62 

Dia. 


5.79 



FIGURE C-4. M-l ABLATIVE CHAMBER FUEL TORUS ASSEMBLY 
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PARAMETER 

ODK INPUT 
NAME 

DESIGN VALUE 

Chamber Radius, R 

c 

- 

20 . 31 in. 

Throat Radius, R^ 

RSTAR 

15.0 in.. 

Contraction Ratio 

ECRAT 

1.833 - 

Inlet Angle, au 

THETAI 

11.3° 

Cylindrical Length, L c 

- 

0.0 in. 

Chamber Length, L-j- 

- 

29.75 in, 

Normalized Inlet Radius 

rwtu/r t < 

RWTU 

2.132 

Normalized Outlet Radius 

rwtd/r t 

RWTD 

.213 

Expansion Angle, aD 

THETA 

29.9° 

Exit Angle, aE 

THE 

29.9° 

Exit Radius 

• - 

21.55 

Expansion Ratio 

EPS 

2.06:1 


Figure c-5 M-l Thrust Chamber Design Parameters 
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5.6 
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12,721 

1,142 

6.658E-05 
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6190, 
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1,136 

7.040E-05 

7.0 

6323. 

14,726 

1.130 

7.337E-0S 
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6344, 

15,937 

1,128 

7.464E-05 

10,0 
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17,930 

1.129 

7 , 462E-05 

12.0 

6053, 

19.514 

1,132 

7.327E-05 


1 

1 0 

10 
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23,874 

.0 

97. 

2.016 

1.4 

.05 

120, 

,06 

4,05 


1295,9 
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1,833 

11.29 

.0 


INJECTION PROCESSES PROGRAM TASK IIB M*l 
TE8T OOt XMINDE64 , 05 


.0 31.96 

TE8T OAT* CORRELATION 


0 

1 

,087614 
,071221 ,1 

,428879 169 , 


.2207 


3248 0 

0 3 

,231 ,57265 -i 

,0 ,021564 , 

169, -.049087 6000, 

.0 97. 2, 

26,09 0,5 0, 

.0 600, 


.0 

97, 

3,0553 

1.4 

0,010 


.0 

.0 

,037854 
4, OS 
1. 


NOTE: FIRST 624 CARDS IDENTICAL TO 
SAMPLE CASE IN REF. 6 


FIGURE C-6. Cl CM INPUT DECK 
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An Experimental-Analytical Method to Study Steady Spray 

Combustion 

Frediano V. Bracco* 

Guggenheim f.aborutoi ies, Princeton University, Princeton, N. J. 


Theme 

A N cxpeiimental-miulytical method is presented by which 
the local combustion gas pat a meters and flux ofliquid fuel 
drops resulting from the steady burning of a fuel spray in a 
gaseous oxidizer can be determined. The method does not 
require any knowledge of the droplet distribution function, 
ding and vaporization equations. Instead, it requites local 
static pressure measurements. Results from the application 
of this method to a liquid oxygen-ethanol rocket combustor 
:uc given. They relate Mostly to the axial uniformity of the 
vaporization tale and of the combustion gas variables. 

Confent 

This method is essentially a technique to obtain maximum 
information out of a set of static pressure measurements based 
on substituting the measured static pressure values into the 
conservation equations and in solving them for other un- 
known quantities. For clarity, the technique is here illus- 
trated using a simplified set of equations which require more 
assumptions than arc necessary. The necessary assumptions 
arc listed after the technique lias been introduced. Details 
about the technique and its extensive application to various 
configui ations of a liquid oxygen-ethanol rocket motor can be 
found m the thesis icfeiied to in the footnote. 

Considci a constant cross-sectional area combustor in which 
a liquid fuel and a liquid oxidizer arc injected. Assume that 
(he oxidiz.er vapoiizes much faster than the fuel and consider 
that part of the combustor where only gaseous oxidizer, com- 
bustion pi oducts, and liquid fuel drops exist. Further assume 
that the combustion is steady, that at the station of interest 
the flow is one dimensional (uniform through the cross sec- 
tion) and that Ihcrc is no iccirculalion. Neglect heat transfer,, 
and viscosity cifects. Temporarily assume also that all 1 fuel 
drops have initially the same velocity and radius and that there 
aic no collisions, bicakups, or nucleations. The follow-ing 
equations, relating pi o perl ies at the injector end to properties 
at any downstieam station, can then be wiitten 

/>«■-( W r Wo>) I- M'o* (1) 

pip Po --p ■ • (IP/tlr - I Fo, //,,,) -b IVojl/jr^ (2) 

pn[h H (u 1 ! 2)]= — [lF,(A f 1 iP, 12 -|- h° t ) - W 0f 

(A 1 « ! o,/2 j-/, 0 ,) ,p 0 J2 t /,„„)] (3) 
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T~ T{p, p, X t ) (4) 

h — h{p, p, Xt) AT) 

Flip, T, X,' 2 . W o r , W 0 >, W F ) - 0 /=-- 1, 2 ... / 

W F jr 3 -= W 0t /r 0 3 (6) 

Hi {(iurldx) ------ JC D p\u - u F \ (n — u F )/rp L (7) 

Uytdrjdx) ~ - A-[j 1- g RSVSr (8) 


Where p , u, p, T, h arc the combustion gas density, velocity, 
pressure, temperature, and latent enthalpy respcclix'cly ( p a is 
the value of p at the injector end). W t (WJ is the local liquid 
fuel (oxidizer) flux and W 0r (W a4 ) is its value at a - 0 (in- 
jector end). Hi (/q>) is the liquid fuel (oxidizer) drop velocity, 
itoiiu'ot) is the injection velocity and u XF iu r ^) is its component 
in the a direction. A,(A«) and /i°r(/i < .°) are the vaporization 
cneigy and tire cnlhalpy of foimation respectively. X, are 
the number of moles of product i per mole of burned fuel, 
r is the local drop radius. p L is the specific gravity of the 
liquid fuel. C D and R e are the drag coefficient and the Rey- 
nolds number respectively and k, x, g, q me properly selected 
constants. Equations (1-3) express mass, momentum, and 
energy consci vation, respectively. Equation (4) is the thermal 
equation of stale of the combustion products, and Eq. (5) is 
the caloric equation of state. F, stands for a set of /equations 
which arc necessary to relate the amount of vaporized propel- 
lants to the variables of the gas (they arc as many as the 
chemical species of winch the gas is assumed to be made up.) 
Equation (6) states the conservation of the drop number. 
Equations (7) and (8) aie possible forms of the diag and 
vaporizalion equations for individual diops. If the con- 
ditions at tjic injector end and basic thci modynamic data arc 
knovvn, these S | / equations contain the following 8 f- / 
unknowns* p, u, p, T, h, W F , tir, r, X,. 2 .. 

It is then observed that the first 5 -|- / equations could be 
solved if any two of the 7 T J unknowns appearing in them 
were given, in which case the last three equations could be 
dropped. Notice that the knowledge of two parameters al- 
lows the elimination of three equations since Eq. (fi) contains 
the drop radius which appears in the last two equations but 
not in (he first 5 /. 

Actually, the measurement of just one parameter is sufficient 
to obtain useful solutions of the first 5 1-/ equations. Indeed 
the terms containing the liquid drop velocity (tti) in Eqs. (2) 
and (3) are small so that the solution of the system is not very 
sensitive lo its value. Accordingly, the i alio 0 *" iii/u <1 1 was 
selected as one of the two parameters. 

The selection of the parameter which is to lie measured is 
dictated by the criterion that it must be easy to measinc and 
the solution of the first 5 1- / equations must ise sensitive to 
its value. The static pressure (actually the loss of viatic pres- 
sure between the injector and any axial location), whkh meets 
those i equipments, cun be selected. 

In conclusion, the first 5 -i - 1 equations can be solved, at 
various axial locations, for selected values of it, ju and using 
static pressure measurements. All the gas variables and the 
liquid fuel flux are thus determined without any assumption 
about the droplet drag, and vaporization processes. 
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j. spacecraft 


Typical results arc presented in Figs. I and 2. They were 
obtained with an oxygen-ethanol rocket motor of constant 
cross sectional area (7.62 cm ID). The injector was made up 
of 16 impinging like-on-like doublets with a distance between 
injector units of 1.5 cm. Chamber pressure, nozzle entrance 
Mach number and injection mixture ratio ( O/F ) were 20 atm, 
0. 1 5 and 2.33, respectively. The static pressure difference be- 
tween the injector and various downstream stations was 
measured accurately and repeatedly by water manometers and 
is given in Fig. 1 (this technique is feasible only for low cham- 
ber pressutes and/or small nozzle entrance Mach numbers), 
The first 5 + / equations were then solved for p, u, T, h, W t 
and the local concentrations of O, H, 0 2 , H 2 , OH, CO, H 2 0. 
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Fig. 2. Calculated dimensionless gas velocity, gas 
temperature, and liquid fuel llu\ \s distance 
from the injector. 


and C0 2 (/ - 8). 1 he calculated local gas velocity , gas temp- 
erature, and flux of liquid fuel are given in Fig 2 where u t and 
T/ are the complete combustion values and ir„ r is the injection 
value. The validity.of the approach and of the assumptions 
embodied in the first 5 -}• /equations wete fut ther checked by 
measuring the gas velocity (by streak photography) at two 
stations (vertical bars in Fig 2 ) Of particular interest is the 
calculated dimensionless liquid fuel flux; a parameter which is 
important for both efficient and stable rocket chambei design. 

A re-examination of the assumptions actually needed for 
the application of this method, leads to the conclusion that 
only two assumptions influence the results markedly. The 
first assumption is that of no recirculation at the station of 
interest (not everywhere between the injector and the station 
of interest). Near the injector, within distances of the order 
of the distance between injector units or of the jet break-up 
length, which ever is longer, recirculation can be expected to 
be active and the first 5 equations do not apply. Indeed in 
Fig. 2 no results at e given foi .v <T 3 cm. The second assump- 
tion is that which must be made to write out explicitly the 
F, equations. In this study, instantaneous mixing and reaction 
of the vapori/ed fuel to equilibrium reaction products was 
assumed. Notice that the assumption that all fuel drops have 
initially the same velocity and radius and that there arc no 
collisions, break-ups, or nucleations ate not necessary as it is 
indicated by the relative insensitivity of the solution to the 
value of the parametci «,/«. 

The objective of this Synoptic has been the explanation of 
the method. However, the results which wete obtained by its 
extensive application to the liquid oxygen-ethanol system are 
also of piactical importance. They are discussed in the foot- 
note and concern both optimal and steady combustion cham- 
ber design and reseatch. The conclusions, which should be 
valid for liquid oxygen-hydrocarbon systems of practical 
interest, include: the assumption of chemical equihbi ium of 
the reaction products appears to be a valid one (it simplifies 
considerably computations), the liquid fuel vaporizes and 
burns uniformly in the axial direction (sec Fig 2, for example) 
rather than actively near the injeetoi and very slowly fat Horn 
it (it limits the usefulness of the eoncentialeJ combustion 
models sometimes used in stability studies); the gas pata- 
nietets are not axially unifotm (as exemplified in 1 ig. 2, it 
explains why the obsetved longitudinal instability shock wave 
frequency is found to be close to the complete combustion 
acoustic chambet frequency); the energy soutce is not pro- 
portional to the mass sou ice (a consequence of the chemical 
equilibtiuiu of the reaction products, it complicates consider- 
ably stability studies); the initial momenta of the liquids ate 
important in steady state computations (they account, foi 
example, foi the observed increase of static pressure neat the 
injector, as shown in Fig 1). 
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NOMENCLATURE LIST 


A 


-B 



A 


D 

.D 

J A 

K 

L' 

M 


ODE 

O/F 

Pc 

Rg D 

T 


U 

V 

V 

W i. 

ZOM 

Z 

A 

% C* 


Area 

Drop size constant 
Drop size constant 
Constant 

Atomization rate constant 
Drag coefficient 
Diameter 

Mass median drop diameter 
Drop secondary breakup constant 

Mass transfer coefficient 
Chamber length 
Flowrate 

One dimensional equilibrium 

Mixture ratio 

Chamber pressure 

Reynold's number 

Temperature 

Velocity 

Velocity 

Weber number 

Local mass flux- 

Cold flow collection plane' distance 

Axial distance 

Difference 

Characteristic exhaust velocity efficiency 
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Nomenclature List (cont.) 


n^* Same as % C* 

a Diffusion correction factor 

a Surface tension, standard deviation 

y . Viscosity 

p Density 

e Chamber contraction ratio 

a Diffusion correction factor 

Subscripts 


a,s 

BL 

C 

d 

D 

eff 

f 

g 

HL 

HW 

j 

Kin 

1 

MIX 

0 

OPP 

PRED 

P 


Vapor pressure at droplet surface 

Boundary layer 

Chamber 

drop 

Drop, diameter 

Effective throat stagnation pressure 

Fuel 

Gas 

Heat loss 
Hot wa.x 
Jet 

Ki net-i cs 
Liquid 
Mixing 
Oxidizer 
Opposite orifice 
Predicted ' 

Propellant 
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Nomenclature List (cont.) 


Subscripts (cont.) 

r Relative 

s Surface, static pressure 

T Throat, total 

TD Two-Dimensional 

Test Test Value 

V Vapor 

VAP Vaporization 
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